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ABSTRACT 


Large  scale  structure  and  mixing  processes  are  investigated  in  chemically 
reacting  wakes  and  shear  layers  to  which  a  periodic  disturbance  is  applied.  The 
experiments  employ  a  diffusion-limited  acid-base  reaction  to  directly  measure 
the  extent  of  mixing.  Optical  diagnostics  used  include  laser  absorption  and 
laser  induced  fluorescence.  Absorption  of  laser  light  by  reacted  product  pro¬ 
vides  a  measure  of  cross-stream  average  product.  Fluorescence  was  measured 
by  a  self-scanning  linear  photodiode  array  using  high  speed  computer  data 
acquisition  to  obtain  the  product  distribution  across  the  layer. 

Previous  results  showing  that  forcing  alters  the  structure  and  growth  rate  of 
shear  layers  are  confirmed.  Forcing  artificially  extends  the  lifetime  of  vortices 
whose  size  is  consistent  with  the  disturbance  wavelength.  Amalgamation  of 
smaller  vortices  is  enhanced  over  that  in  the  natural  layer  until  the  frequency 
locked  scale  is  achieved.  At  high  Reynolds  number  product  measurements  show 
reduction  of  product  with  forcing.  At  moderate  Reynolds  numbers,  on  the  other 
hand,  there  is  an  increase  in  product  when  forced.  In  one  case  a  five  fold 
increase  in  product  was  observed.  The  differences  are  related  to  the  different 
effects  of  forcing  on  entrainment,  composition  ratio  and  secondary  structure.  • 

A  dramatic,  order  of  magnitude  increase  in  mixing  was  discovered  for  certain 
forced  wake  flows.  This  effect  is  strongly  associated  with  an  interaction  between 
the  spanwise  organized  wake  vortices  and  the  test-section  side  walls. 
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Chapter  1 


Introduction 


1.1.  Historical  Background 

The  experiments  reported  in  this  thesis  are  the  latest  in  a  series  of  investiga¬ 
tions  of  the  mixing  layer  which  began  with  the  discovery  by  Brown  and  Roshko 
(1971,  1974)  that  organized  vortical  structures  play  a  prominent  role  in  those 
shear  flows,  even  when  they  are  turbulent  at  high  Reynolds  numbers.  Another 
important  step  occurred  when  Oster  and  Wygnanski  (1982)  found  that  those 
large  structures  could  be  controlled  by  introducing  relatively  small  periodic  per¬ 
turbations  into  the  flow  and  that  the  development  of  the  shear  layer  could 
thereby  be  drastically  altered.  In  the  meantime,  Konrad  (1978)  and  Breidenthal 
(1978),  (1981)  had  been  investigating  the  relation  of  those  structures  to  mixing 
at  the  smallest  scales,  i.e.,  molecular  mixing.  It  had  become  evident  that  the 
primary  structures  play  a  prominent  role  in  bringing  fluids  from  the  two  sides 
mto  the  mixing  layer,  the  so-called  "entrainment”  process,  but  it  was  not  clear 
how  the  ultimate,  ‘‘chemical”  mixing  is  accomplished  and  how  it  depends  on 
Reynolds  number.  Their  principal  finding  was  that,  in  flows  which  develop  from 
laminar  boundary  layers,  the  initial  Kelvin-Helmholtz  instability  and  the  large, 
primary  vortices  which  result  from  it  are  followed,  at  some  distance  down¬ 
stream,  by  a  “mixing  transition''  where  the  amount  of  mixed  (or  chemically 
reacted)  fluid  in  the  shear  layer  dramatically  increases  although  the  interac¬ 
tions  and  mechanics  of  the  primary  vortices  do  not  appear  to  be  greatly  or  qual¬ 
itatively  altered.  The  enhanced  mixing  was  attributed  to  increase  of  interface 
area  between  fluids  from  the  two  sides,  created  by  secondary  and  smaller  scale 
motions.  Konrad  found  that  after  the  mixing  transition  the  rate  of  production 
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of  mixed  fluid  reached  what  seemed  to  be  a  plateau,  remaining  independent  of 
increasing  downstream  distance,  and  thus  of  increasing  local  Reynolds  number, 
within  the  accuracy  of  the  measurements. 

The  next  obvious  question  following  from  the  measurements  of  Konrad  and  of 
Breidenthal  and  the  experiments  of  Oster  and  Wygnanski  is  the  following:  what  is 
the  effect  of  periodic  forcing  on  the  chemical  mixing?  Is  it  similar  to  the  effect 
on  growth  rate,  which  is  connected  with  momentum  mixing?  Certainly  the 
development  of  negative  Reynolds  stress  will  not  have  a  counterpart  in  species 
unmiying!  The  main  aim  of  the  present  experiments  was  to  study  these  effects. 

1.2.  Recent  Developments 

Experimental  investigations  more  recent  than  the  pioneering  ones  discussed 
above  are  also  relevant  to  the  present  work. 

Some  insights  into  the  structure  of  the  secondary,  stream  wise  vortices,  which 
precede  and  contribute  to  the  mixing  transition,  was  obtained  by  Bernal  (19B1) 
who.  in  particular,  obtained  unexpected  views  of  the  cross-sections  of  those  vor¬ 
tices.  using  laser  induced  fluorescence  (L.I.F.)  techniques  developed  in  this 
laboratory  by  Dimotakis,  et  al.  (1983). 

Konrad's  measurements  through  the  mixing  transition  were  made  in  gas 
flows  while  Breidenthal's  were  in  water.  The  post-mixing-transition  plateau  in 
water  was  found  to  be  lower  by  a  factor  of  two  than  in  air  (Figure  1.1).  Some  of 
this  could  be  attributed  to  a  systematic  experimental  effect  (see  Breidenthal, 
1978)  but  also  to  the  effects  of  Schmidt  number  (Broadwell  and  Breidenthal. 
1981)  which  is  about  0.7  for  air  and  600  for  water.  More  recent  measurements 
by  Mungal  (1983)  in  gas  flow,  using  a  chemically  reacting  shear  layer,  defined  a 
somewhat  lower  post-mixing-transition  plateau  level  than  Konrad's:  Kooches- 
fahani,  working  in  water  with  a  technique  different  from  Breidenthal’s, 
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measured  a  lower  post-mixing-transition  plateau  Level  than  Breidenthal’s.  The 
ratio  between  the  gaseous  and  aqueous  plateau  valves  is  approximately  2  for  the 
Konrad-Breidenthal  set  (Figure  1.1)  and  about  the  same  for  the  Mungal- 
Koochesfahani  set.  but  the  latter  are  about  20%  lower  than  the  former.  Thus 
there  is  some  uncertainty  about  the  correct  values  for  mixing  product.  What  is 
not  controversial  and  relevant  for  our  investigations  is  that  there  is  a  mixing 
transition. 

Returning  to  the  aspects  of  the  problem  concerning  forcing  effects,  an  impor¬ 
tant  contribution  was  made  by  Ho  and  Huang  (1982).  They  forced  shear  layers 
by  superimposing  periodic  perturbations  on  the  free  stream;  Oster  and  Wygnan- 
ski  had  utilized  an  oscillating  flap  at  the  trailing  edge  of  the  splitter  plate. 
Whereas  Oster  and  Wygnanski's  experiments  were  at  high  Reynolds  number,  and 
they  observed  post-mixing-transition  flow,  Ho  and  Huang's  were  mainly  in  pre¬ 
mixing-transition  flows. 

1.3.  Variables  and  Parameters  of  the  Problem 

In  the  unforced  layer,  the  mixing  effects  are  usually  expressed  in  terms  of  the 
growth  of  a  thickness  cS(x);  this  may  be  any  measure  of  the  thickness  based  on 
velocity  profile  or  on  visual  appearance  or,  especially  in  this  investigation,  the 
reaction  or  mixing  product  thickness  6p  (Konrad  and  Breidenthal  called  it  P). 

The  rate  of  growth,  or  of  production,  44  is  then  a  measure  of  the  mixing 

dx 

rates.  For  a  self-similar  shear  layer,  growing  linearly  with  x,  these  rates  are 
constant.  Of  interest  is  how  they  depend  on  parameters  of  the  flow  and  under 
what  circumstances  are  they  not  constant  in  i,  an  important  point  in  question 
being  the  mixing  transition. 

For  unforced  layers  in  incompressible,  homogeneous  flows,  the  effect  of 
.  UZ 

different  velocities  in  the  two  streams  -rr~  =  r  can  usually  be  normalized  by  the 
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ratio 


Ux-Uz 

Ui  +  Uz 


to  which  the  post-mixing-transition  growth  rates  are  proportional.  This  leaves 
the  possible  effects  of  viscosity,  which  are  sometimes  expressed  as  a  local  Rey- 
( U,  -  U»)x 

nolds  number  - - - to  non-dimensionalize  the  development  distance  x  ; 

alternatively,  the  downstream  distance  can  be  made  nondimensional,  — ,  with 

an  initial  length  such  as  the  momentum  thickness  9X,  of  the  splitter-plate  boun¬ 
dary  layers.  With  this  the  role  of  Reynolds  number  is  expressed  somewhat 
indirectly,  through  its  role  in  establishing  the  value  of  in  a  particular 
geometry,  the  direct  effect  of  Reynolds  number  should  then,  in  principle,  require 

Ux6\ 

another  parameter  such  as  — - — .  Of  course,  if  the  Reynolds  number  is  “high 
enough"  it  is  often  assumed  that  there  is  no  direct  effect  of  Reynolds  number. 

When  periodic  perturbation  or  forcing  is  applied  to  the  shear  layer,  additional 
parameters  are  introduced  into  the  flow:  the  perturbation  amplitude,  say  U\ 
which  could  be  normalized  by  Ux  or  Ux  —  U2  :  and  the  forcing  frequency  F. 
This  latter  was  established  by  Oster  and  Wygnanski  as  a  very  prominent  parame¬ 
ter  for  forced  layers.  In  presenting  the  effects  of  F  on  growth  rates  <5(x).  they 
found  that  the  effects  of  F  on  a  given  shear  layer  could  be  correlated  by  using 

Ue  , 

the  normalizing  length  —  where  Ue  ■  jr(C/x  +  U2)  is  taken  to  be  the  convec- 

r  2 

tion  velocity  of  the  large  organized  vortices.  Thus  rr—  is  a  function  of  ~r— .  To 

Ue  Ue 

( £7i  -  Ug) 

correlate  the  effects  of  U2/U\,  the  growth  parameter  c  a  -rr - — =r-  is  incor- 

\U\  +  Ug) 


porated  so  that, 


=  /n(xwo) 


(1.1) 
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where  we  have  introduced  what  we  call  the  "Wygnanski-Oster  parameter”  X wr). 

Oster  and  Wygnanski  did  not  actually  show  such  plots;  that  was  done  by 
Browand  and  Ho  (1983),  who  incorporated  the  Oster-Wygnanski  and  the  Ho- 
Huang  results  into  a  plot  which  is  reproduced  in  Figure  1.2.  The  dependent 
thickness  here  is 


/ 


{Uj-UHU-  Uz) 
(U1  -  Uz)z 


dy 


This  figure  displays  the  dramatic  effects  of  forcing  on  the  growth  rate  and  exhi¬ 
bits  the  three  different  response  regions  of  the  shear  layer,  which  had  already 
been  introduced  by  Oster  and  Wygnanski.  The  characteristics  of  these  regions 
are,  briefly,  as  follows: 

I.  Xwo  <  1  ;  the  growth  rate  is  enhanced  by  up  to  a  factor  of  2 
(reached  at  a  saturating  value  of  forcing  amplitude) 

II.  1  ss  Xwo  £  2;  in  this  "frequency  locked”  region  growth  rate  is 
inhibited,  even  reduced  to  zero;  Reynolds  stresses  are  reversed. 

III.  Xwo  >  2;  relaxation  to  unforced  growth  rate. 

The  significance  of  Region  II  is  that  the  passage  frequency  of  the  vortices 
(which  are  all  equally  spaced  in  that  region)  is  equal  to  the  forcing  frequency. 
The  spacing  or  scale  of  the  vortices  entering  that  region  has  been  detern  med  by 
whatever  amalgamation  processes  are  needed  to  accomplish  this  in  region  1. 
The  forcing  wave  length  is  large  compared  to  the  scales  of  region  I  and  this  has 


the  effect  of  enhancing  growth  in  that  region. 

The  Browand-Ho  plot  in  Figure  1.2  gives,  ideally,  a  single  correlation  curve 
and  it  is  rather  impressive  that  both  pre-mixing-transition  and  post-mixing- 
transition  data  fall  quite  well  onto  it.  Nevertheless,  one  should  expect  to  see,  for 
a  wide  range  of  conditions,  some  dependence  on  other  (dimensionless)  parame¬ 
ters.  Indeed,  a  later  plot  prepared  by  Ho  and  Huerre  (1964)  have  on  it  other 
data  which  do  not  all  fall  onto  the  '‘universal”  curve  so  well. 

On  the  basis  of  our  experiments  to  be  described,  we  have  chosen  as  an  addi¬ 
tional  parameter  Fq.  the  frequency  of  the  initial  Kelvin-Helmholtz  structure  in 
the  shear  layer.  This  is  in  lieu  of  the  initial  pare  meter  6X  which  we  discussed 
above  in  connection  with  unforced  layers;  in  fact,  F0  is  proportional  to  Uc/ex. 
Finally,  the  parameter  characterizing  direct  Reynolds  number  effects  would  be 
Re#i.  so  that  equation  1.1  generalizes  to 

f  = /n <*„;£.  Re. ,) 

The  plot  of  Browand  and  Ho,  which  is  influenced  heavily  by  the  high  Reynolds- 
number  results  of  Oster  and  Wygnanski,  should  be  appropriate  for  the  asymp- 

F 

totic  conditions  Re®,  »  1  and  <<  1-  The  latter  may  also  be  expressed  as 

■~»  1,  i.e.,  the  forcing  wave  length  is  large  compared  to  initial  scales. 

Ao 

In  the  experiments  to  be  described  we  attempted  to  define  a  corresponding 
plot  for  the  distribution  of  reaction  product,  ip(x).  along  the  layer  at  the 
asymptotic  conditions  above,  but  were  somewhat  Limited  by  the  range  of  condi¬ 
tions  available  in  the  flow  facility  we  used;  in  fact,  a  rather  large  range  of  condi- 

F 

tions  extending  to  low  Reynolds  number  and  values  of  -  1  were  covered. 

f  o 
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1.4.  Outline 

A  brief  description  of  the  flow  facility  and  the  measurement  techniques  is 
given  in  Chapter  2. 

The  major  part  of  the  experiments  was  devoted  to  mixing  layers,  the  results 
for  which  are  presented  in  Chapter  3,  the  main  body  of  this  report.  Here  the 
effects  of  forcing  on  the  large  vortex  structure,  on  growth  rate,  and  on  mixing 
product  are  described,  together  with  their  relation  to  the  mixing  transition.  - 

Chapter  4  is  a  briefer  description  of  the  effects  of  forcing  on  the  wake  of  the 
splitter  plate,  i.e.,  the  flow  to  which  the  mixing  layer  degenerates  when  (7a  =  Ux, 
which  had  also  been  investigated  by  Breidenthal. 

A  short  discussion  in  Chapter  5  summarizes  the  observations,  measurements 
and  some  of  their  ramifications. 

Appendix  A  is  a  short  note  on  some  aspects  of  the  laser  technique  for 
measuring  product  thickness.  Appendix  B  was  motivated  by  the  visual  pictures 
which  were  obtained.  It  is  a  model  computation  of  a  mixing  layer  utilizing  point 
vortices  to  represent  the  spanwise  organized,  large  structures.  The  main  pur¬ 
pose  was  to  compute  the  patterns  produced  by  the  interface  between  the  two 
sides  of  the  layer  as  it  is  rolled  up  and  distorted  by  the  action  of  the  vortices,  for 
comparison  with  those  obtained  in  the  laboratory. 
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Chapter  2 

Experimental  Facilities  and  Instrumentation 

2.1.  Breidenthai's  Mixing-Layer  Apparatus 

The  flow  facility  constructed  by  Breidenthal  (1978)  was  used  for  all  the  exper¬ 
iments.  Some  minor  improvements  were  made  in  the  honeycomb  and  screen 
configuration  in  an  effort  to  reduce  the  freestream  turbulence  level. 
Breidenthal  reported  a  turbulence  intensity  of  0.5%  and  our  measurements 
confirmed  this.  Prof.  Hassan  M.  Nagib,  with  his  considerable  experience  in  such 
matters,  suggested  some  changes  that  reduced  the  turbulence  intensity  to  0.3% 
for  the  high  Reynolds  number  flows  we  studied. 

The  screens  were  replaced  with  ones  made  from  smaller  diameter  wire  spaced 
on  a  coarser  mesh  yielding  a  65%  open  area  ratio  thereby  reducing  the 
anomalous  effects  common  to  screens  possessing  60%  or  less  open  area  ratio. 
The  six-inch-long  honeycomb  section  was  replaced  with  a  pair  of  two-inch  honey¬ 
combs  spaced  two  inches  apart.  Each  honeycomb  had  a  chopper  screen  on  the 
downstream  end.  The  reader  is  refered  to  Loehrke  and  Nagib  (1976)  for  the  rea¬ 
soning  behind  these  changes. 

A  reference  scale  for  the  photographs  throughout  this  work  is  provided  by 
the  test  section  dimensions;  7  cm  high  x  11  cm  span  x  45  cm  length. 

2.2.  The  Forcing  Mechanism 

A  considerable  part  of  the  effort  in  the  preliminary  phase  of  the  work  was 
aimed  at  producing  a  controllable  periodic  disturbance.  Wygnanski  and  Oster 
(1980)  forced  the  flow  with  an  oscillating  flap  that  was  hinged  to  the  trailing 
edge  of  the  splitter  plate  while  Ho  and  Huang  (1982)  used  rotating  valves  in  the 


supply  lines  to  each  stream  to  provide  periodic  disturbances  to  the  free 
streams.  For  the  present  experiments,  it  was  important  that  the  forcing 
mechanism  should  not  limit  optical  access  or  add  complexity  to  test  section 
maintenance.  For  this  reason  a  flapping  splitter  plate  tip  was  not  considered 
viable.  Instead,  a  periodic  pressure  drop  device  seemed  a  better  choice.  The 
forcing  mechanism  provided  a  periodic  loss  in  head  which  produced  a  periodic 
component  in  the  mass  flow  through  the  test  section. 

The  periodic  pressure  drop  was  produced  by  a  variable  area  orifice  plate.  The 
drop  in  pressure  across  an  orifice  plate  depends  on  the  ratio  of  area  of  the 
orifice  to  area  of  the  pipe  in  which  it  is  placed;  the  smaller  the  ratio  the  greater 
the  pressure  drop.  The  forcing  mechanism  emulated  an  orifice  plate  whose  area 
ratio  was  a  periodic  function  of  time.  This  was  accomplished  by  a  fixed  stator 
which  was  a  thin  plate  with  an  off-center  circular  hole  and  a  rotor  which  was  a 
disc  with  a  pie-shaped  wedge  cut  out  of  it.  Vs  the  rotor  turned  at  constant 
speed,  the  flow  saw  an  open  area  that  was  a  periodic  function  of  time.  Figure  2.1 
shows  the  rotor  and  stator  viewed  from  the  downstream  side  looking  upstream. 
Figure  2.1a  shows  the  largest  open  area  and  Figure  2.1b  demonstrates  the  smal¬ 
lest  area  for  the  particular  rotor  and  stator  shown  in  the  photographs.  The 
rotor  was  powered  by  a  DC  electric  motor  and  the  speed  of  rotation  was  moni¬ 
tored  with  an  optical  pick-up. 

Figure  2.2  shows  the  complete  forcing  mechanism  with  drive  motor  and  opti¬ 
cal  pick-up  as  it  appears  before  installation.  The  device  was  inserted  from  above 
into  the  drain  pipe  leading  to  the  high  speed  settling  chamber  and  supported 
with  angle  iron  braces  as  shown  in  Figures  2.3  and  2.4.  The  motor  was  powered 
by  a  DC  power  supply  and  the  pulse  train  from  the  optical  pick-up  was  moni¬ 
tored  with  a  counter  from  which  frequency  was  inferred.  The  frequency  was 
very  easy  to  change  as  it  was  only  necessary  to  literally  turn  a  knob  on  the 
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power  supply  to  change  rotation  speed  and  therefore  frequency.  The  amplitude 
and  shape  of  the  waveform  were  changed  by  either  changing  stators  (different 
circular  cut-outs  with  different  off-center  dimensions)  or  rotors  (different  size 
pie-shaped  cut-outs).  Trial  and  error  produced  a  reasonable  disturbance  with  at 
least  90%  of  the  total  energy,  as  measured  with  a  spectrum  analyzer,  in  the  pri¬ 
mary  purely  sinusoidal  component  and  amplitudes  less  than  2 %  of  the  frees- 
tream  high  speed  velocity.  The  forcing  was  always  applied  to  the  high  speed 
stream  only. 

2.3.  Laser  Doppler  Velocimeter 

A  Laser  Doppler  Velocimeter  (LDV)  and  phase-locked-loop  processor  was  con¬ 
structed.  complete  with  three-dimensional  traverse  and  processing  electronics 
to  allow  accurate,  non-intrusive,  velocity  measurements.  The  LDV  and  phase- 
locked  loop  were  copies  of  successful  designs  perfected  by  Prof.  Paul  Dimotakis 
and  his  students  in  this  laboratory. 

A  photograph  of  the  LDV,  traverse,  and  test  section  is  included  in  Figure  2.5. 

The  short-duration,  blow-down  nature  of  the  facility  did  not  allow  detailed 
mean  velocity  measurements  on  the  high-speed  flows.  Measurements  of  this 
type  were  made  only  on  the  very  low  speed  laminar  shear  layer.  The  LDV  was 
used  mainly  for  setting  the  speeds  prior  to  a  run  and  for  measurement  of  the 
amplitude  of  the  forcing  disturbance. 

2.4.  Measurements  of  Structure  and  Growth  Rates  from  Photographs 

Flow  visualization  was  used  extensively  in  these  experiments.  As  found  by 
past  investigators,  it  is  an  attractive  technique,  providing  a  relatively  simple, 
almost  unique  method  for  seeing  and  recording  detailed  spatial  structure  in  the 
flow  We  used  passive  dye  marking  as  well  as  reaction-product  marking  tech¬ 
niques  to  obtain  either  short  exposure  photographs  to  "freeze"  instantaneous 
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spatial  structure  or  long  time  exposures  to  give  measures  of  mean  growth  rates. 
Two  different  visualization  techniques  were  employed.  The  first  was  the  phenol- 
phthalein  method  of  Breidenthal  (197B)  and  the  second  was  the  laser  induced 
fluorescence  method  of  Dimotakis,  Miake-Lye  and  Papantomou  (1962).  Both 
allow  one  to  visualize  the  product  of  a  chemical  reaction  (i.e.,  mixed  fluid)  or  to 
visualize  concentration  of  a  particular  reactant  (i.e.,  as  a  passive  contaminant). 

Simple  photographic  methods  were  employed  to  record  the  photographs. 
Three  different  35mm  cameras  were  used;  a  Nikon  Model  F,  a  Pentax  ME  Super, 
and  a  Konica  FS-1.  All  were  essentially  equivalent  as  regards  performance. 
Ektachrome  160  and  Kodachrome  40  color  slide  films  were  used  to  photograph 
the  phenolphthalein  experiments.  The  laser  induced  fluorescence  experiments 
required  "faster”  film  so  Tri-X  black-and-white  film  gave  better  results  than 
color  slide  film.  Time  exposures  with  high  contrast  were  difficult  to  obtain  with 
the  available  lighting.  A  very  dark  green  filter  (Kodak  Number  54)  and 
Panatomic-X  black-and-white  film  were  used  when  taking  time  exposures  of  the 
red  reacted  phenolphthalein.  The  Alter  passed  green  and  blue-green  light  while 
stopping  red  light  almost  entirely.  Acceptable  time  exposures  of  as  long  as  ten 
seconds  were  obtained  by  this  method. 

Because  snap  shots  require  only  a  fraction  of  a  second,  flows  where  measure¬ 
ments  could  be  made  only  tediously  or  not  at  all  could  at  least  be  photographed. 
The  many  photographs  obtained  during  the  course  of  this  research  were  care¬ 
fully  cataloged  for  future  reference.  The  phenomena  documented  ranged  from 
truly  laminar  flows  to  turbulent  post-mixing  transition  shear  layers.  Wake  flows 
(U i  =  Uz)  were  also  studied  for  a  wide  range  of  Reynolds  numbers  by  photo¬ 
graphic  means. 

Some  quantitative  information  was  obtained  from  time  exposures.  The  edge 
of  the  mixed  region  in  a  time  exposure  was  selected  and  digitized  with  a  Hen  .ett- 
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Packard  9874A  digitizer  at  several  points  along  the  top  and  bottom  of  the  side 
view  wedge  of  reacted  fluid.  The  resulting  data  were  curve  fitted  to  an  eight 
oarameter  function  using  a  least  squares  method  to  produce  a  continuous  curve 
on  top  and  bottom  that  represented  the  long  time  average  edges  of  the  mixed 
region.  The  perpendicular  distance  between  edges  was  measured  and  found  to  be 
nearly  equal  to  the  visual  thickness, 

2.5.  Measurements  of  Reaction  Product  by  Light  Absorption 

The  absorption  technique  of  Breidenthal  (1978)  was  used  to  obtain  quantita¬ 
tive  information  about  the  mixing.  For  this  technique  and  the  laser  induced 
fluorescence  technique  to  be  described  later,  the  product  of  a  chemical  reaction 
was  measured.  A  pH  sensitive  dye  was  placed  in  one  stream  and  mixed  with  an 
acid  such  that  the  dye  was  in  the  "ofH  3tate,  i.e.,  was  transparent.  A  base  was 
added  to  the  other  stream.  Hence,  the  pure  fluids  are  transparent  but  when 
mixed  at  a  molecular  level,  the  base  causes  the  pH  indicator  to  turn  "on"  thus 
marking  the  product  of  the  reaction. 

For  the  absorption  technique,  the  pH  indicator  was  phenolphthalein  which  is 
transparent  in  the  "ofT1  state  and  red  in  the  "on"  state.  Breidenthal  used  a 
mercury-argon  discharge  lamp  as  a  light  source  and  passed  a  beam  through  the 
test  section  perpendicular  to  the  plane  of  the  layer,  thus  measuring  the 
integrated  concentration  of  product  along  the  line  of  sight  of  the  beam  as  a 
function  of  time.  The  technique  was  modified  for  the  present  work  in  the  follow¬ 
ing  way:  an  argon-ion  laser  was  used  as  a  source  and  the  beam  was  spread  mto 
a  sheet  so  that  a  spanwise  average  was  measured. 

The  exponential  relationship  between  incident  and  transmitted  light  was  the 
basis  for  calculating  the  product  thickness.  For  a  sheet,  it  was  necessary  to 
approximate  the  exponential  function  if  the  measurement  was  to  be  interpreted 
as  spanwise  average  product.  The  details  of  how  this  was  done  and  other 
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relevant  information  on  the  absorption  method  are  discussed  in  Appendix  A. 

During  the  absorption  runs,  four  channels  of  data  were  recorded;  the  output 
of  the  photodiode,  the  input  laser  power,  the  free  stream  velocity,  and  a  refer¬ 
ence  pulse  train,  were  all  digitized  at  a  500  hertz  rate  and  multiplexed.  Thirty 
seconds  of  data  were  taken  at  sixteen  downstream  locations  for  eight  different 
flows.  The  data  acquisition  computer  was  designed  and  built  by  D.  Nosenchuck. 
The  reader  is  referred  to  his  thesis,  Nosenchuck  (1962),  for  details. 

2.6.  Measurements  of  Structure  and  Reaction  Product  by  L.I.F. 

Koochesfahani,  Dimotakis  and  Broadwell  (1963)  have  shown  that  fluorescence 
intensity  can  be  used  to  observe  product  concentration  in  a  turbulent  shear 
layer.  This  technique  uses  a  pH  sensitive  dye.  fluorescein,  which  fluoresces 
efficiently  in  the  presence  of  green  laser  light  (the  5145  angstrom  line  of  a  Lexel 
Model  95  Argon-ion  laser  was  used).  The  intensity  of  fluorescence  is  linearly  pro¬ 
portional  to  concentration  of  dye  in  the  measuring  volume.  The  laser  beam  was 
passed  through  the  shear  layer  perpendicular  to  its  plane  and  imaged  with  a 
camera  lens  onto  a  self-scanning  linear  photodiode  array  called  a  "Reticon 
Array*'.  The  end  result  is  that  it  was  possible  to  measure  product  concentration 
as  a  function  of  time  at  1024  points  simultaneously  in  an  entirely  non-intrusive 
way.  These  points  were  equally  spaced  along  the  height  of  the  test  section. 

The  present  experiments  using  this  technique  were  conducted  in  collabora¬ 
tion  with  M.  Koochesfahani  who  developed  the  hardware  and  data-aquisition 
software  for  this  application.  The  reader  is  referred  to  his  thesis.  Koochesfahani 
(1984),  for  detailed  discussion  of  the  limitations  of  this  technique. 

One  limitation  worth  mentioning  is  the  resolution  problem.  Each  photodiode 
in  the  array  provided  a  measurement  which  was  an  average  in  two  ways.  First,  it 
was  an  average  over  space.  The  optics  collected  light  from  a  small  volume  rf 
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fluid  and  focused  it  onto  the  detector  which  simply  added  up  the  total.  For  this 
experiment  a  length  scale  for  the  volume  was  about  0.1  millimeters.  Second,  it 
was  an  average  in  time.  The  current  from  the  diode  charged  a  capacitor  which, 
m  a  sense,  was  a  bucket  that  collected  and  added  up  the  current  flow  for  a  short 
period  of  time  and  then  dumped  it  out  when  the  measurement  was  made.  The 
time  scale  was  the  period  of  the  sampling  rate  which  in  this  case  was  2  mil¬ 
liseconds.  The  smallest  length  scale  in  the  flow,  it  turns  out,  was  always  smaller 
than  the  measuring  volume  for  conditions  in  this  experiment.  For  pre-mixing- 
transition  flows  the,  smallest  scale  was  the  interface  thickness  (i.e.,  molecular 
diffusion  length)  which  is  small  for  a  liquid  (estimated  to  be  0.3^im  in  these 
flows).  For  post-mixing-transition  flows  the  length  associated  with  the  fine  scale 
turbulence,  namely  the  Kolmogorov  scale,  is  small  and  the  interfacial  diffusion 
scale  in  even  smaller.  Both  were  estimated  to  be  smaller  than  the  measuring 
volume  for  the  Reynolds  numbers  studied  in  this  work.  This  resolution  difficulty 
must  be  kept  in  mind  when  interpreting  the  fluorescence  data. 

All  of  the  laser  induced  fluorescence  data  were  taken  at  a  location  25  centim¬ 
eters  downstream  of  the  splitter  plate  tip.  This  represented  about  2/3  of  the 
test  section  length.  A  total  of  39  runs  were  made.  Approximately  20  seconds  of 
data  were  taken  for  each  run,  meaning  that  10,000  data  points  were  recorded 
for  each  diode  (pixel)  on  the  array  Only  every  fourth  pixel  was  sampled,  so  out¬ 
put  from  256  pixels  was  recorded,  implying  that  over  two  million  data  points 
were  recorded  for  each  20  second  run.  The  data  were  stored,  temporarily,  on 
hard  disc  and  later  transfered  to  mag  tape  for  permanent  storage.  Processing 
was  conducted  on  our  computer,  a  PDP  11/44  based  system. 
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2.7.  Measures  of  Shear  Layer  Thickness 

Throughout  this  irork,  reference  is  made  to  many  different  measures  of  thick¬ 
ness.  All  of  them  have  been  introduced  by  previous  investigators  and  we  have 
used  the  conventional  symbol  for  each  whenever  possible.  The  following  is  a  list 
of  the  thicknesses  along  with  a  description  of  each. 


Symbol  Description 

6,  6U  Vorticity  thickness  of  the  mean  velocity  profile.  It  is  computed 

from  u{y )  as  follows: 


V  i  ~  Uz 

du 

,  ,  max 


Within  the  assumption  of  the  boundary  layer  approximation.  5  is 
thus  physically  the  same  as  the  maximum  slope  thickness. 


Integra  thickness  of  the  mean-velocity  profile.  It  is  computed 


from  iZ(y)  as  follows: 


a-  f 


g(y)-cfr  u  i-g(y) 

Vx-Vz  U I  -  u2 


9  is  often  incorrectly  refered  to  as  the  momentum  thickness 
because  it  appears  to  be  a  generalization  of  the  formula  for 
momentum  thickness  of  a  boundary  layer  appropriate  for  a  shear 
layer.  For  U2  =  0  it  is  true  that  6  is  momentum  thickness.  How¬ 
ever,  for  U2*  0  it  can  be  shown  that  6  is  unrelated  to  momentum. 


Momentum  thickness  of  the  high  speed  boundary  layer.  This  is  the 
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<5, 


conventional  thickness  used  to  represent  initial  conditions.  It  can 
be  computed  from  the  measured  boundary  layer  mean-velocity 
profile  in  the  usual  way.  In  the  case  of  a  laminar  boundary  layer, 
Thwaite's  method  can  be  applied  to  get  an  estimate  of 
Breidenthal  computed  by  this  method.  Throughout  the  present 
study,  was  computed  using 


Si 


l,xet  t 

Ur 


j. 

2 


(2.3) 


where  xett  is  an  effective  development  length  equal  to  5  centime¬ 
ters  for  our  facility. 

1 %  thickness.  The  1%  here  refers  to  1%  of  the  maximum  mean 
product  concentration.  <5i  is  the  distance  from  the  1%  level  on  the 
low  speed  side  to  the  1%  level  on  the  high  speed  side.  Mungal 
(1983)  and  Koochesfahani  (1984)  have  found  it  is  very  nearly  equal 
to  the  visual  thickness. 


<5,*,  Visual  thickness.  Brown  and  Roshko  (1974)  introduced  this  thick¬ 

ness  as  a  means  of  extracting  quantitative  information  from  pho¬ 
tographs.  It  is  the  thickness  obtained  by  drawing  straight  lines 
tangent  to  the  top  and  bottom  edges  of  the  structures  in  a  short 
exposure  photograph.  In  the  present  work,  the  edge  in  a  long  time 
exposure,  taken  with  the  method  described  earlier,  was  nearly 
equal  to  the  visual  thickness  and  is  thus  refered  to  as  . 

Product  thickness  referenced  to  free  stream  conditions  on  the 
high  speed  side.  Mungal  (1983)  and  Koochesfahani  (1984)  present 
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their  product  measurements  in  terms  of  this  thickness  and  6pg. 
Mungal  used  the  symbols  Px  and  Pz  for  and  <5pg  respectively. 
For  <5pt  the  reference  concentration  is  the  high  speed  free  stream 
value.  It  is  the  width  of  a  top-hat  profile  whose  height  is  the  free- 
stream  concentration,  Cl0,  and  whose  area  is  the  same  as  the 
mean  product  concentration  profile. 


7r- /<7(y)  *y 

u  lo 


(2.4) 


BE 

L  ■ 

l  - 
[ 


&pm  <  Pe  .  P  Product  thickness  referenced  to  free  stream  conditions  on  the  low 
speed  side.  This  is  the  same  as  6Pl  except  the  reference  conditions 
are  on  the  low  speed  side. 

Spt-  ’^■fc^(y)dy  (2.5) 

dpa  is  the  same  as  P  in  Breidenthal's  work  and  Pz  in  Mungal's  work. 
In  the  present  work,  the  marker  fluid  in  all  the  shear  layer  experi¬ 
ments  was  carried  by  the  low  speed  fluid  so  when  we  speak  of  pro¬ 
duct  thickness  we  mean  6Pa. 
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Chapter  3 

Effects  of  Forcing  on  Shear  layer  Structure  and  Mixing 

In  this  chapter  several  experiments  on  the  mixing  layer  are  presented  and 
discussed.  Emphasis  is  on  the  effects  of  forcing.  Changes  m  the  large  scale 
structure  and  in  streamwise  development  are  discussed  for  a  large  range  of  Rey¬ 
nolds  numbers.  The  effects  of  forcing  on  the  mixing  characteristics  of  shear 
layers  are  discussed  in  detail. 

Before  discussing  the  effects  of  forcing,  the  basic  structure  and  mixing 
phenomena  in  the  unforced  shear  layers  studied  in  this  investigation  are 
described.  This  provides  a  necessary  reference  for  the  discussion  of  the  results 
obtained  with  forcing 

3.1.  Reynolds  Number  Dependence  of  the  Unforced  Layer 

A  wide  range  of  flow  speeds  was  studied  during  the  course  of  this  work.  For  a 
fixed  velocity  ratio,  the  initial  Reynolds  number  of  the  high  speed  boundary 
layer  characterizes  changes  in  flow  conditions  for  changing  speeds.  This  Rey¬ 
nolds  number  is  defined  as  follows: 

UX9X 

Re*t  =  — ~  O-D 

where  Ux  is  the  high  speed  velocity  and  ex  is  the  momentum  thickness  of  the 
high  speed  boundary  layer.  This  Reynolds  number  is  obviously  connected  with 
initial  conditions.  Also  useful  is  the  unit  Reynolds  number  defined  as  follows: 

Re'  =  — 


(3  2) 


Re'  has  dimensions  of  inverse  length  and  when  multiplied  by  the  vorticity  thick¬ 
ness  becomes  the  large-structure  Reynolds  number.  The  range  of  unit  Reynolds 
number,  Re',  was  from  approximately  100  cm'1  to  4340  cm'1.  Correspondingly, 
the  range  of  initial  Reynolds  number  which  could  be  achieved  varied  over  a 
range  from  30  to  190.  Photographs  of  shear  layers  at  three  different  Reynolds 
numbers  within  this  range  are  shown  in  Figure  3.1  (the  y  dimension  of  the  test 
section  in  these  photographs  is  7  cm).  They  provide  examples  of  the  three  flow 
regimes  through  which  a  shear  layer  developes  with  downstream  distance,  start¬ 
ing  from  an  initially  laminar  state:  (i)  steady  laminar  flow;  (ii)  pre-mixing  transi¬ 
tion  flow;  and  (iii)  post-mixing  transition  flow, 

In  Figure  3.1a,  the  unit  Reynolds  number  is  low  enough  for  the  laminar  part 
of  the  layer  to  extend  over  the  whole  field  of  view.  The  mixing  occurs  on  a  thin 
interface  separating  high  speed  and  low  speed  fluid.  The  mean  velocity  profile 
for  the  flow  in  Figure  3.1a,  measured  with  the  LDV  at  a  location  20  cm  down¬ 
stream.  is  shown  in  Figure  3.2.  At  this  location  the  maximum  slope  thickness  of 
the  profile  is  0.6  mm  and  the  edge-to-edge  thickness  is  approximately  12  mm, 
while  the  diffusion  interface  thickness  in  Figure  3.1a  is  less  than  1  mm. 

Of  course,  the  diffusion  layer  will  be  much  thinner  than  the  vorticity  layer 
because  the  former  begins  at  the  trailing  edge  of  the  splitter  plate  while  the 
latter  has  developed  from  the  boundary  layers  leaving  the  plate  (estimated  to 
have  an  edge-to-edge  thickness  of  4mm).  If  the  observed  thickness  of  12mm.  at 
x  =  20  cm  were  to  correspond  to  a  shear  layer  beginning  with  zero  thickness  at 
the  same  origin  as  the  diffusion  layer,  then  the  available  similarity  solution 

1 2  771771 

implies  a  diffusion  layer  thickness  approximately  equal  to  -  j— -,  where 

Sfc* 

Sc  =  is  the  Schmidt  number;  for  water,  the  diffusion  layer  thickness  would 


be  approximately  0.5  mm. 
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With  increasing  unit  Reynolds  number,  instability  appears  at  some  point 
downstream.  In  Figure  3.1b  that  point  is  at  the  far  left  of  the  picture  (a  similar 
example  may  be  seen  in  Figure  3.6a).  Figure  3.1b  displays  another  feature,  the 
mixing  transition,  at  about  the  midpoint  of  the  picture.  Thus  this  picture 
displays  all  three  flow  regimes  defined  above. 

In  Figure  3.1c  the  mixing  transition  (not  well  defined  in  this  picture)  is  com¬ 
plete  by  about  x  -  9 cm  (as  estimated  from  results  discussed  in  Section  3  3), 
while  the  transition  from  the  steady  regime,  if  any,  has  been  squeezed  into  a 
small  region  near  the  splitter  plate. 

The  above  examples  show  that  the  unforced  shear  layer  goes  through  two 
transitions.  The  first  is  the  transition  from  steady,  laminar  flow  to  a  two- 
dimensional  flow  with  growing  instabilities  which  eventually  form  vortices  (as  in 
Figure  3  6a).  Above  the  Reynolds  number  where  vortices  form,  no  further 
change  in  primary  structure  takes  place.  The  first  transition  is,  therefore,  a 
transition  in  structure.  The  second  transition  is  a  transition  in  mixing.  The 
exact  details  of  how  this  transition  takes  place  are  not  well  understood  but  it  is 
associated  with  development  of  three-dimensional  structure  It  is  known  that 
an  order  of  magnitude  increase  in  mixing  occurs  in  the  high  Schmidt  number 
case,  e  g.,  Breidenthal  (1978). 

The  fact  that  the  two  transitions  occur  at  different  values  of  the  Reynolds 
number  means  that  three  flow  regimes  exist.  The  first  is  the  Laminar  shear 
layer  This  occurs  when  the  Reynolds  number  is  smaller  than  both  transition 
values.  For  the  laminar  shear  layer,  simple  diffusion  of  vorticity  and  mass  in  the 
shear  layer  are  the  dominant  phenomena  In  the  second  regime  the  flow  is  a 
turbulent  shear  layer  with  Reynolds  number  smaller  than  that  required  for  the 
mixing  transition,  i.e. ,  the  Reynolds  number  is  larger  than  the  structure- 
transition  Reynolds  number  but  is  smaller  than  the  mixing-transition  Reynolds 


number. 


We  shall  refer  to  this  as  the  "pre-mixing- transition  turbulent ”  regime. 
Growth  rate  and  momentum  mixing  by  Reynolds  stress  are  all  well  developed. 
Mixing  is  limited  by  diffusion  over  an  interface  in  this  case.  The  final  case  is  the 
one  where  both  transition  Reynolds  numbers  are  exceeded.  This  is  the  "posf- 
mucing- transition  turbulent"  regime.  Entrainment  by  inviscid,  large-scale 
structures  is  the  limiting  process  which  determines  mixing  in  this  case. 

3.2.  Forced  Shear  Layer  Structure  —  Visual  Results 

In  this  section,  results  obtained  from  photographs  of  reacting,  forced  shear 
layers  are  presented.  Photographs  provide  a  way  to  discern  changes  in  large- 
scale  structure,  and.  to  a  limited  degree,  changes  in  mixing  that  take  place  due 
to  forcing. 

Prior  to  discussion  of  the  photographs,  some  terminology  is  needed  to  sim¬ 
plify  the  discussion.  When  the  Reynolds  number  is  larger  than  that  required  for 
the  transition  in  structure,  discussed  above,  we  refer  to  the  flow  as  “turbulent”. 
When  the  Reynolds  number  is  smaller  than  the  structure  transition  value  the 
flow  is  "laminar".  To  identify  regions  of  the  flow  with  respect  to  the  mixing  tran¬ 
sition,  we  will  use  the  terms  "pre-”  and  "post-mixing-transition”.  Figure  3.1a  is 
thus  a  laminar  shear  layer.  Figure  3.1b  is  primarily  a  turbulent,  pre-mixing- 
transition  shear  layer  though  the  downstream  half  of  the  picture  is  post-mixing 
transitional.  Figure  3.1c  is  a  turbulent,  post-mixing-transition  shear  layer. 

Each  shear  layer,  which  is  initially  laminar  when  it  leaves  the  splitter  plate 
has  an  initial  Instability  and  roll-up  frequency  which  will  be  called  the  "natural 
frequency”.  It  is  the  passage  frequency  of  the  first,  and  smallest,  organized 
structures  in  the  mixing  layer.  Freymuth  (1966)  and  Ho  and  Huang  (1982) 
found  that  this  frequency  is  predicted  well  by  simple  parallel-flow,  linear 
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stability  theory  applied  to  the  initial  mean  velocity  profile;  c.f.,  Michalke  (1964) 
for  the  classical  treatment  of  a  hyperbolic  tangent  profile.  Michalke’s  temporal 
theory  implies  that  the  most  unstable  wavelength,  A0,  is  28 .5  times  the  momen¬ 
tum  thickness,  9 l 


A0  =  28.5  0!  (3.3) 

The  natural  frequency,  Fa,  is  equal  to  the  convection  velocity  divided  by  the 
most  unstable  wavelength, 


Throughout  the  present  work,  the  natural  frequency  was  estimated  using  these 

U !  +  U2 

relationships  together  with  an  estimate  of  9,  from  Eq.  (2.3)  and  Uc  =  - - - . 

For  the  later  discussion  of  the  effects  of  forcing,  it  is  useful  to  recall  some 
characteristics  of  the  spectra  of  passage  frequencies  of  the  large  vortical  struc¬ 
tures  (Winant  and  Browand,  1974)  or,  equivalently,  the  probability  distributions 
of  their  spacings  (Brown  and  Roshko,  1974)  Near  the  splitter  plate,  in  the  early 
stages  of  a  shear  layer  developing  from  laminar  boundary  layers,  the  (nearly) 
fixed,  natural  frequency  defines  a  spectral  peak  (or  delta  function  in  the  wave¬ 
length  or  vortex-spacing  distributions)  (Ho  and  Huang,  1982).  Eventually  far 
downstream,  in  the  ”well  developed”  turbulent  shear  layer,  a  broad  spectrum 
develops.  Its  maximum  corresponds  to  the  most  probable  passage  frequency 
(or  vortex  spacing);  it  is  close  to  the  mean  value  (Brown  and  Roshko.  1974) 
Consistent  with  a  linear  growth  rate  for  the  shear  layer,  the  local  passage  fre- 

i 

quency  decreases  downstream  as 

Th<--  forced  shear  layer  of  Wygnanski  and  Oster,  described  in  the  Introduction, 
was  in  the  latter,  late  stages  of  development  (post-mixing-transition)  while  the 


forcing  experiments  of  Ho  and  Huang  were  mainly  on  a  shear  layer  in  Its  early 
stages  of  development,  i.e..  pre-mixing-transition. 

F 

The  work  of  Ho  and  Huang  (1982)  focused  mainly  on  — -  —  1  while  Oster  and 

F0 

F 

Wygnanski  (1982)  studied  «  1.  The  present  research  covered  the  entire 

.To 

range  and  the  observed  effects  of  forcing  on  structure  agree  well  with  the 
findings  of  the  above  investigations.  Much  of  what  follows,  therefore,  may  be 
found  in  the  above  works.  The  aim  in  this  discussion  is  to  unify  the  results,  to 
present  the  necessary  background  for  those  unfamiliar  with  the  previous 
research  and.  most  important,  to  provide  a  reference  for  the  effects  of  forcing 
on  the  large-scale  structure  and  growth  rate  of  the  shear  layers  for  which  pro¬ 
duct  measurements  were  made. 

First  consider  the  case  when  the  forcing  frequency  equals  the  natural  fre- 
jr 

quency,  i.e.,  •=—  =  1.  In  this  case,  vortices  are  shed  from  the  splitter  plate  at  the 
•f*o 

natural  frequency  and  move  downstream  in  an  orderly  fashion,  without  amalga¬ 
mation,  for  an  unnaturally  long  distance.  Eventually,  the  growth  of  subharmon¬ 
ics  of  the  input  frequency  causes  pairings,  triplings,  etc.,  like  those  observed  in 
the  unforced  layer.  Figure  3.3  shows  am  example  of  a  forced,  pre-mixing- 
transition  shear  layer.  Laser  induced  fluorescence  was  used  to  label  pure  fluid 
from  the  low  speed  side  for  this  photograph,  The  forcing  wavelength 
corresponds,  in  Figure  3.3,  to  the  spacing  between  the  initial  structures  on  the 
left,  i.e.,  the  upstream  side.  The  onset  of  a  pairing  is  evident  on  the  downstream 
end.  Even  though  energy  is  being  put  into  the  flow  at  the  initial  roll-up  fre¬ 
quency,  the  growth  of  subharmonics,  available  in  the  background,  eventually 
causes  a  pairing. 


-  24  - 


Forcing  at  frequencies  larger  than  the  natural  frequency,  i.e  .  —  >  1,  causes 

*  0 

no  change  in  the  behavior  outlined  above.  Vortices  form  at  the  forcing  fre¬ 
quency  and  move  downstream  without  amalgamation  until  subharmonics  grow 
sufficiently.  The  flow  adjusts  itself  to  the  forcing  in  this  case  by  a  simple  reduc¬ 
tion  of  the  flow  scale  to  accommodate  the  input  wavelength.  There  must  exist 
an  upper  limit  to  the  forcing  frequency  at  which  this  is  true.  Linear  stability 
theory  predicts  that  the  neutral-point  frequency  is  roughly  twice  the  natural 
frequency  High  frequencies  were  difficult  to  achieve  in  the  present  experiments 
so  this  topic  was  not  pursued. 


Forcing  at  frequencies  smaller  than  the  natural  frequency,  —  <  1,  on  the 

c  0 

other  hand,  causes  a  variety  of  more  interesting  flow  phenomena.  In  a  sense, 
this  is  the  more  important  case  from  an  engineering  standpoint  since  most 
shear  layers  of  technical  interest  are  large  Reynolds  number  flows  for  which  the 
natural  frequency  is  very  high  and  background,  excitation  spectra  are  more 
likely  to  contain  lower  frequencies;  the  initial  roll-up  frequency  {F0)  can  be  so 


high  in  fact  that  forcing  it  may  not  be  practical.  Behavior  of  the  structures  for 
tr 

—  <  1  is  best  understood  by  first  considering  forcing  at  integer  subharmonics 
F0 


of  the  natural  frequency,  i.e., 


F  1 

—  =  —  where  n  is  an  integer  The  first  subhar- 
F  o  n 


monic  corresponds  to  n  =2,  the  second  subharmonic  to  n  =  3,  etc. 


If  the  flow  is  forced  at  exactly  one-half  of  the  natural  frequency  .  — —  =  a 

r  o  2 

controlled  pairing  takes  place.  Vortices  are  shed  from  the  splitter  plate  at  the 
natural  frequency  and  immediately  undergo  a  pairing  to  form  a  larger  vortex. 
The  larger  vortex  moves  downstream  without  interaction  along  with  other  vor¬ 
tices  of  the  same  type  In  our  facility,  it  was  never  obvious  that  these  larger  vor¬ 
tices  would  themselves  undergo  a  pairing  before  they  reached  the  end  of  the 
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test  section.  Figures  3.4a  and  3.4b  show  the  difference  between  forcing  the 
natural  frequency  and  the  first  subharmonic. 

Amalgamations  of  more  than  two  vortices  can  be  caused  by  forcing  at  higher 
integer  subharmonics  of  the  natural  frequency.  Figure  3.4c  and  3.4d  show 
examples  of  controlled  tripling  and  controlled  quadrupling.  The  forcing  fre¬ 
quency  is  one-third  the  natural  frequency  in  Figure  3.4e  and  one-quarter  the 
natural  frequency  in  Figure  3.4d. 

At  this  point,  the  effect  of  forcing  on  structure  has  been  outlined  for  forcing 
at  frequencies  larger  than  or  equal  to  the  natural  frequency  and  for  forcing  at 
integer  subharmonics.  Now  consider  what  happens  at  Intermediate  frequencies. 
No  new  types  of  flow  structure  occur.  Instead,  rescaled  versions  of  the  flow  pat¬ 
terns  discussed  above  are  the  result. 

This  is  best  explained  by  example.  Consider  the  case  when  the  input  fre¬ 
quency  is  less  than  the  natural  frequency  but  is  greater  than  the  first  subhar¬ 
monic.  Recall  that  forcing  at  the  natural  frequency  causes  no  amalgamations 
and  forcing  at  the  first  subharmonic  causes  pairing.  For  forcing  frequencies 
between  the  two,  a  rescaled  version  of  the  flow  forced  at  the  natural  frequency 
occurs.  This  was  first  observed  by  Ho  and  Huang  (1982)  and  is  confirmed  by  the 
present  research.  The  "response  frequency1’  of  the  shear  layer,  a  term  intro¬ 
duced  by  Ho  and  Huang  (1902),  is  the  frequency  at  which  vortices  form  at  the 
splitter  plate.  When  the  forcing  frequency  is  between  the  natural  frequency  and 
the  first  subharmonic,  the  "response  frequency”  equals  the  forcing  frequency. 
All  flows  forced  between  the  natural  frequency  and  the  first  subharmonic,  i.e., 
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—  <  <  1.  thus  have  a  flow  structure  that  is  a  rescaled  version  of  the  flow 
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resulting  from  forcing  at  the  natural  frequency.  The  rescaling  is  such  that  the 
spacing  between  the  vortices  shed  from  the  splitter  plate  is  equal  to  the 


wavelength  of  the  forcing.  Hence,  the  scale  of  the  flow  is  increased  for 


In  a  similar  manner,  forcing  between  the  first  subharmonic  and  second 
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subharmonic,  —  <  — -  <  —  causes  a  rescaled  pairing  flow  pattern.  In  this 
3  Fa  2 

case,  the  response  frequency  is  twice  the  forcing  frequency  The  rescaling  is 
such  that  spacing  between  the  vortices  shed  from  the  splitter  plate  is  equal  to 
one-half  the  forcing  wavelength.  Similarly,  forcing  between  the  second  subhar¬ 
monic  and  third  subharmonic  causes  a  rescaled  tripling  flow  pattern  and  so  on. 

From  the  foregoing  discussion,  one  may  be  tempted  to  extend  the  findings  to 
infinite  subharmonics  and,  in  so  doing,  develop  a  recipe  for  the  flow  resulting 
from  any  arbitrary  frequency  3maLler  than  the  natural  frequency.  For  example, 
the  recipe  alluded  to  would  imply  that  forcing  at  the  19th  subharmonic  causes  a 
controlled  amalgamation  of  20  small  vortices  into  a  single  large  vortex.  How¬ 
ever,  for  large  numbers  of  vortices,  the  amalgamation  does  not  take  place  in  a 
single  event  as  it  does  when  smaller  numbers  join.  Depending  on  the  availability 
of  disturbances  in  the  background,  a  series  of  pairing,  tripling,  etc.,  interactions 
may  take  place  until  a  row  of  large  vortices,  in  tune  with  the  forcing  frequency, 
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is  set  up.  This  behavior  is  characteristic  of  the  -r—  <<  1  regime  and  was  the 
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basic  result  of  Wygnanski  and  Qster  (1980).  Ho  and  Huang  (1982)  refer  to  the 
case  when  large  numbers  of  vortices  come  together  in  a  single  event  as  a  "col¬ 
lective  interaction"  and  state  that  large  amplitudes  are  required  to  produce  it. 
In  this  work,  collective  interactions  were  difficult  to  generate  because  the  top 
and  bottom  of  the  test  section  limit  the  flow. 

If  consideration  were  limited  to  small  amplitudes,  which  are  more  important 
from  a  practical  standpoint,  then  the  suggested  recipe  would  be  misleading 


when  the  forcing  frequency  is  very  small  compared  to  the  natural  frequency, 
—  «  i.  For  this  case,  the  results  of  Wygnanski  and  Oster  (1980)  show  that  a 
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series  of  amalgamations  occurs  until  vortices  in  tune  with  the  input  frequency 
are  developed.  This  occurs  at  the  downstream  location  where  the  mean  local 
passage  frequency  of  the  large-scale  structures  in  the  unforced  flow  is,  roughly, 
equal  to  the  forcing  frequency.  Browand  and  Ho  (1903)  find  that  the  growth 
rate  upstream,  in  the  interaction  region,  may  be  as  much  as  twice  that  of  the 
natural  shear  layer,  depending  on  the  amplitude  of  the  forcing  (see  Figure  1.2). 
In  this  sense,  one  can  say  that  the  forcing  causes  the  early  amalgamations  to 
occur  at  an  unnaturally  high  rate.  The  upstream  region  of  enhanced  interac¬ 
tion  is  called  Region  I  by  Wygnanski  and  Oster  (1980).  The  subsequent  region, 
where  the  amalgamated  vortices  are  in  tune  with  the  forcing,  is  referred  to  as 
Region  II.  The  layer  possesses  basically  no  growth  in  Region  II.  Eventually,  again 
depending  on  availability  of  disturbances  in  the  background,  the  frequency- 
locked  vortices  begin  to  interact.  Region  III,  as  Wygnanski  and  Oster  (1980)  cadi 
it,  is  this  downstream  region  where  a  relaxation  to  the  unforced  state  takes 
place.  Eventually,  at  far  enough  distances  downstream,  the  original  growth  rate 
is  achieved  (see  Browand  and  Ho  (1983)  or  Figure  1.2). 

To  summarize,  the  effects  of  forcing  on  structure  described  above  can  be 

F 

separated  into  two  categories.  The  first  case  occurs  when  is  near  unity  but 
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not  very  small  compared  to  unity.  In  this  case,  controlled  amalgamations 
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result.  The  second  case  occurs  for  <<  1.  In  this  case,  the  flow  can  be 

F o 

divided  into  three  regions,  as  shown  in  Figure  1.2.  In  the  first  region  (I),  the 
forcing  causes  an  unnaturally  high  growth  rate  due  to  enhanced  interactions  of 
the  vortices.  In  the  second  region  (II),  large  scale  vortices  are  locked  to  the 
forcing  frequency.  No  amalgamations  and  virtually  no  growth  characterize  this 
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region.  Finally,  in  the  third  region  (III),  the  vortices  from  Region  II  begin  to 
interact  and  the  growth  rate  returns  to  the  unforced  value. 


Examples  of  Forced  Shear  Layers  for  Various  Reynolds  Numbers.  With  the 
basic  effects  of  forcing  on  structure  having  been  outlined,  the  remainder  of  this 
section  will  be  devoted  to  examples  of  the  various  flow  regimes  outlined  above. 
The  examples  are  taken  from  the  many  photographs  of  the  flow  recorded  during 
the  course  of  this  research.  The  range  of  available  forcing  frequencies  was  lim¬ 
ited  to  approximately  1  Hz  to  20  Hz.  Low  Reynolds-number  flows  have  low 
natural  frequencies  and  high  Reynolds  numbers  have  high  natural  frequencies 
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so,  for  the  most  part,  low  Reynolds-number  flows  fall  into  the  — — -  l  category 
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while  high  Reynolds-number  flows  have  -=—  <<  1.  The  examples  will  begin  with 
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the  smallest  Reynolds  number  studied  and  end  with  the  largest  value 


The  first  example  is  the  steady  laminar  shear  layer  of  Figure  3.1a.  Figure  3.5 
shows  photographs  of  the  same  flow  when  forced  at  F  -  0.5 He.  Although  the 
shear  layer  in  Figure  3.1a  is  not  visibly  oscillating,  possibly  because  oscillations 
have  not  yet  sufficiently  amplified,  an  estimate  of  the  initially  most  unstable  or 
“natural"  frequency  can  be  made  from  equations  2.3.  3.3  and  3  4  Although  the 
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Reynolds  number  is  low  ( — - —  *  60)  the  result  from  inviscid  stability  theory  is 
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good  enough  for  our  purposes.  This  gives  F$  =  0.35  Hz .  thus  -=—  =  1.4.  Alterna¬ 
te) 
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tively,  using  conditions  further  downstream,  at  x  =  20  cm,  — —  =  2.0.  Figure  3.5a 
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shows  simultaneous  plan  and  side  views  of  the  flow  while  Figure  3.5b  shows  the 
view  from  above  the  plane  of  the  layer  looking  upstream.  The  photographs  show 
that  forcing  produces  a  wave-like  distortion  of  the  diffusion  interface.  The  wave 
shape  is  approximately  two-dimensional  in  the  central  part  of  the  test  section 
with  wave  crests,  i.e.,  the  dark  sections  in  the  plan  view,  parallel  to  the  splitter 


plate.  Near  the  side  walls,  the  crests  split  into  symmetric  sets  of  upstream  and 
downstream  tilted  crests.  The  side  wall  interaction  is  due  to  the  relatively  thick 
side-wall  boundary  layers  present  in  this  low  speed  flow.  A  simple  estimate  of 
the  side-wall  boundary  layer  thickness  is  consistent  with  measurements  taken 
from  the  photograph. 

Figure  3.5  demonstrates  that  forcing  alters  the  flow  and  might  even  enhance 
mixing.  The  fact  that  the  crests  are  dark  in  the  photograph  indicates  that  pos¬ 
sibly  more  product  is  present  there.  Whether  or  not  forcing  increases  the  mean 
product  is  speculative  since  product  measurements  were  not  made  in  this  flow. 

The  next  example  is  the  previously  cited  Figure  3.3.  The  flow  speeds  for  this 
case  are  roughly  four  times  those  for  the  laminar  flow  discussed  above.  The 
forcing  frequency  is  greater  than  the  first  sub-harmonic  so  vortices  form  from 
the  splitter  plate  at  the  forcing  frequency  and  move  downstream  in  an  orderly 
fashion.  The  visualization  technique  is  a  non-reacting  one  in  which  pure  high 
speed  fluid  is  black  and  pure  low  speed  fluid  is  white.  The  nearly  vertical  black 
lines  in  the  white  area  are  shadows  caused  by  bubbles  on  the  upper  surface  of 
the  test  section.  Figure  3.3  shows  qui';e  clearly  that  the  cores  of  the  structures 
contain  pure  high  speed  fluid.  Appendix  B  presents  a  model  of  this  flow  which 
reproduces  many  of  the  features  of  the  photograph.  The  model  is  an  extension 
of  a  model  proposed  by  Jimenez  (1980)  in  which  the  vorticity  is  modeled  with 
point  vortices  and  the  interface  is  passive,  i.e.,  contains  no  vorticity.  Jimenez 
considered  the  case  when  the  interface  intersected  the  centers  of  the  vortices 
and  found  that  the  resulting  roll-up  was  symmetric.  The  model  presented  in 
Appendix  B  considers  an  off-center  interface.  The  roll-up  is  non-symmetric  for 
this  ease  and  the  interface  shapes  are  remarkably  similar  to  Figure  3.3.  Appen¬ 
dix  B  provides  a  complete  presentation  of  the  solution  and  more,  in  depth  dis¬ 
cussion  of  the  results. 
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The  next  example  is  a  turbulent,  pre-mixing-transition  shear  layer.  Figure 
3,6  shows  photographs  of  a  reacting  shear  layer  for  several  different  forcing  fre¬ 
quencies.  The  high  speed  velocity  is  11.6 cm/sec  and  the  low  speed  velocity  is 
5  9cm /sec  for  all  entries  in  Figure  3.6.  From  these,  estimates  of  the  other 
relevant  parameters  give  r  =  0.51,  Uc  -  8.75  cm.  /sec,  e,  =  0.66  mm,  =  4.7  Hz, 

=  76,  and  ^-^-  =  570 cm-1.  Each  photograph  contains  simultaneous  plan 

and  side  views  of  the  flow.  Only  two-thirds  of  the  plan  view,  including  one  wall,  is 
shown.  An  inert  dye  was  added  to  the  lower  free-stream  so  that  pure  high  speed 
fluid  appears  grey  in  the  photographs.  Reacted  fluid  appears  black  in  both  side 
and  plan  views. 

The  changes  in  primary  structure  caused  by  forcing  are  evident  in  the  photo- 
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graphs  and  are  consistent  with  the  description  given  earlier  for  the  — -  ~  1. 

Figures  3.6a,  3.6b,  and  3.6d  are  all  rescaled-natural-frequency  type  flows,  i.e., 
vorticies  form  at  the  splitter  plate  at  the  same  frequency  as  the  forcing  fre¬ 
quency.  This  is  consistent  with  the  earlier  statement  that  the  response  fre- 
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quency  equals  the  forcing  frequency  for  —  <  —  <  1.  Figure  3.6b  shows  that 
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the  frequency-locked  vortices  eventually  begin  to  interact;  a  change  in  scale, 
caused  by  pairing  of  previously  frequency-locked  vortices,  is  evident  in  the  mid¬ 
dle  of  the  side  view.  Figures  3.6e,  3.6f,  and  3.6g  are  examples  of  controlled  amal¬ 
gamations.  It  is  clear  in  these  cases  that  the  vortices  which  participate  in  the 
amalgamation  are  not  of  equal  size.  Figure  3.6f,  for  example,  shows  pairing  of 
large  vortices  with  smaller  upstream  vortices.  In  Figure  3  6g,  the  smaller  vortex 
is  more  clearly  visible  while  in  Figure  3  6e,  the  smaller  vortex  is  a  mere  kink  in 
the  braid  connecting  the  larger  vortices. 

The  development  of  three-dimensional  instabilities  is  quite  interesting  in  the 
series  of  pictures  in  Figure  3.6  The  unforced  case,  Figure  3.8a,  has  three- 


dimensional  structure  but  it  is  hard  to  say  whether  this  is  caused  by  interaction 
with  the  thick  side-wall  boundary  layers  or  if  it  represents  the  onset  of  the  ins¬ 
tability  that  leads  to  the  streamwise  vortices.  The  forced  flows  show  a  variety  of 
interesting  three-dimensional  structure.  In  Figure  3.6b,  the  primary  structures 
at  the  upstream  end  are  seen  to  have  a  certain  waviness  even  at  the  start. 
Notice  also  In  this  figure  how  the  change  in  9cale  in  the  middle  of  the  picture  is 
accompanied  by  development  of  streamwise  structures.  Figures  3.6c  and  3.6d 
show  an  interesting  side  wall  effect  which  occured  frequently  for  the  low  Rey¬ 
nolds  number  forced  flows  encountered  during  this  research.  The  interaction  of 
the  primary  vortices  with  the  side-wall  boundary  layer  produces  three- 
dimensional  motion  which  separates  the  flow  into  a  wall  region  and  a  central 
region.  In  the  central  region,  the  primary  structure  is  spanwise  coherent  while 
in  the  wall  region,  the  vorticity  gets  tipped  by  the  velocity  defect  of  the  boun¬ 
dary  layer.  The  boundary  between  the  two  regions  consists  of  a  complicated 
three-dimensional  motion  which  apparently  mixes  the  two  reactants  efficiently. 
Another  example  of  this  type  of  flow,  taken  at  similar  conditions,  is  shown  in  Fig¬ 
ure  3.7.  Continuing  with  the  plan  views  of  Figure  3.8,  it  is  seen  that  the  side-wall 
influence  is  not  as  strong  in  3.6e,  3.6f,  and  3.8g.  It  appears,  in  3.6g,  that  the 
forced  pairing  causes  the  development  of  streamwise  vortices  at  the  down¬ 
stream  end  of  the  plan  view.  This  may  Indicate  that  pairing  of  the  primary 
structures  triggers  the  development  of  the  secondary  structures. 


In  the  previous  example,  the  mixing  transition  is  not  an  issue,  i.e.,  the  Rey¬ 
nolds  number  was  small  enough  that  the  mixing  transition  did  not  occur  in  the 
test  section.  The  next  example,  shown  in  Figure  3.8.  is  at  a  large  enough  Rey¬ 
nolds  number  for  the  mixing  transition  to  occur  near  the  downstream  end  of 
the  test  section.  The  speeds  for  this  flow  are  C/,  =  17  4  cm /sec  and 
C/2  =  0.3  cm/sec  and  implied  by  them  are  the  following:  r  =  0.48,  =  0.54mm, 


Fq  a  8.3  Hz , 


U  1*1 


=  94,  and 
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905  cm-1.  Black  areas  in  the  photographs 
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indicate  regions  where  mixing  occurs.  However,  the  wedge-shaped  dark  area  at 
the  downstream  end  of  the  plan  views  is  a  shadow  from  an  optical  source.  Quan¬ 
titative  measurements,  presented  later,  do  not  extend  that  far  downstream. 
Unlike  the  previous  example,  the  optics  were  arranged  to  give  the  full  wall-to- 
wall  plan  view  of  the  test  section.  Changes  in  the  primary  structure  for  this  case 
are  very  similar  to  those  discussed  earlier  for  Figure  3.6.  Figures  3.8a,  3.8b,  and 
3.8c  are  rescaled  natural  frequency  flows.  Figures  3.8d  and  3.8e  show  controlled 
pairing.  Figures  3.3f  and  3.8g  show  controlled  tripling,  and  so  on  for  the  remain¬ 
ing  entries  in  Figure  3.8.  A  rough  feeling  for  the  way  forcing  affects  growth  rate 
of  the  layer  can  be  obtained  from  these  photographs.  The  unforced  layer 
demonstrates  the  expected  linear  growth  rate  while  the  forced  layers  each  con¬ 
tain  a  region  where  there  is  virtually  no  growth.  In  each  case,  this  region 
corresponds  to  the  location  where  the  passage  frequency  of  the  primary  vortices 
equals  the  forcing  frequency.  As  the  frequency  decreases  for  succesive  photo¬ 
graphs  in  Figure  3.8,  the  frequency-locked  region  moves  further  downstream. 

As  mentioned  above,  the  mixing  transition  begins  to  occur  in  the  unforced 
case  at  the  downstream  end  of  the  photograph.  Figure  3.8a.  A  picture  alone  is 
not  enough  to  make  such  an  assessment  but  product  measurements,  to  be 
presented  shortly,  verify  that  this  is  the  case.  The  clear  presence  of  the  secon¬ 
dary,  streamwise  vortices  is  a  clue,  however,  that  the  mixing  transition  is 
imminent.  Figure  3.9  is  an  enlargement  of  the  photograph  in  Figure  3.8a  which 
allows  the  three-dimensional  structure  of  the  unforced  case  to  be  more  care¬ 
fully  scrutinized.  Crude  measurements  of  the  streamwise  structure  spacing 
were  obtained  from  this  photograph  and  others  like  it.  The  result  is  that  the 
ratio  of  the  streamwise  to  spanwise  structure  spacing  varied  from  0.3  to  1.0 
which  is  consistent  with  the  range  of  values  measured  by  Bernal(  1981).  It  is 
clear,  particularly  in  Figure  3.8g  and  3.8h,  that  forcing  can  enhance  the  secon¬ 
dary  vortices.  Notice  how  clearly  defined  the  streamwise  structures  are  in  these 
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two  cases.  An  enlargement  of  Figure  3.8g  is  shown  in  Figure  3.10. 

A  fined  note  with  regard  to  the  plan  views  of  Figure  3.8  concerns  the  reduction 
in  the  side-wail  influence  compared  to  the  previous  example  shown  in  Figure  3  0 
The  reason  for  this  is  not  entirely  clear.  The  speeds  are  roughly  a  factor  of  two 
larger  for  Figure  3.8  which  means  the  side  wall  boundary  layers  are  only  30% 
smaller  than  those  in  Figure  3.6.  But  the  flow  near  the  side-walls  in  the  plan 
views  is  drastically  different  in  the  two  cases.  Perhaps  the  more  fully  developed 
three-dimensional  structure  in  Figure  3.8  prevents  the  side-wall  interaction 
from  developing  as  strongly. 

As  a  final  example,  a  high  Reynold-number  shear  layer  will  be  discussed.  This 
was  the  largest  Reynolds  number  at  which  data  were  acquired.  Figure  3.11 
shows  four  side  view  photographs  of  fluorescent  product;  one  unforced  and 
three  forced.  Two  cylindrical  lenses  were  used  to  spread  the  laser  beam;  one  in 
the  spanwise  direction  and  one  in  the  streamwise  direction.  This  produces  a 
laser  beam,  normal  to  the  shear  layer,  with  dimensions  5  cm.  spanwise  and 
40  cm  streamwise.  The  spanwise  width  was  the  same  width  used  in  the  absorp¬ 
tion  measurements.  One  observes,  in  each  picture,  an  instantaneous  view  of  the 
spanwise  averaged  product  in  the  center  5  cm  of  the  span.  The  actual  height  of 
the  test  section  in  each  picture  is  7  cm  and  the  flow  is  from  left  to  right  with  the 
high  speed  side  on  top. 


Figure  3.11a  shows  the  unforced  flow.  Three  forced  shear  layers  and  the 
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value  of  for  each  are  shown  in  Figure  3.11.  Figure  3.11b  shows  the  three 

regions  of  growth  discussed  earlier.  At  the  left  of  the  picture,  two  small  vortices 
are  pairing;  they  will  later  form  a  larger  vortex  in  the  frequency-locked  region. 
Notice  the  apparent  lack  of  structure  in  the  recovery  region  (right  side  of  the 
picture).  It  is  expected  that  large  scale  structures  would  eventually  emerge. 
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The  time  scale  necessary  for  this  to  take  place  depends  on  the  availability  of 
background  disturbances  at  long  wavelengths  and  on  the  growth  rate  of  instabil¬ 
ities  associated  with  those  wavelengths.  It  is  not  unreasonable  for  this  time 
scale  to  be  so  long  that  we  were  unable  to  visualize  the  emergence  of  larger 
scales  in  the  length  of  the  test  section. 

Conditions  in  the  frequency  locked  region  correspond  to  the  neutral  point  on 
the  stability  curve  (Ho  and  Huerre,  1984).  Figure  3  11b  provides  the  necessary 
information  to  check  this  assertion.  The  thickness  was  measured  in  the  locked- 
in  region  of  Figure  3.11b  by  the  method  described  in  Chapter  2.  This  thickness 
along  with  the  known  forcing  frequency  and  convection  velocity, 
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Ue  -  0.5  (U1  +  U2y  imply  a  measured  value  of  dimensionless  wavelength.  — 

approximately  equal  to  0.62.  This  value  is  nearly  equal  to  the  value  at  the  neu¬ 
tral  stability  point  and  is  consistent  with  the  results  presented  by  Ho  and 
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Huerre'  (1984)  if  one  assumes  — — »  8.  When  the  forcing  frequency  is  at  the 

neutral  point,  the  first  subharmonic  is  thus  very  near  the  most  unstable  fre¬ 
quency  and  is  almost  certain  to  emerge  and  cause  amalgamation  of  frequency 
locked  vortices.  Such  events  occur  at  the  begining  of  Region  III. 

Figure  3.11c  shows  the  locked-in  region  most  clearly.  This  picture  is  very 
similar  to  photographs  by  Wygnanski  and  Oster  (1980)  taken  at  Reynolds 
numbers  ten  times  as  high.  Notice  the  penetration  of  unmixed  fluid  all  the  way 
across  the  layer  in  the  third  vortex  from  the  left.  This  suggests  that  the  flow  is 
very  nearly  two-dimensional  at  this  point;  remember  the  picture  is  a  spanwise 
average.  The  spiral  arms  of  the  structures  are  much  longer  and  wrap  around 
more  of  the  core  than  in  the  unforced  case. 
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Figure  3.1  Id  gives  a  good  example  of  the  collective  interaction  of  Ho  and 
Huang  (1982)  which  is  typical  for  small  forcing  frequency  and  large  enough 
amplitude.  Six  vortices,  visible  at  the  left  of  the  picture,  are  wrapped-up  into  a 
single  large  vortex  by  the  time  they  reach  the  downstream  end.  It  is  also  evident 
in  this  picture  that  the  bottom  wall  exerts  an  influence  on  the  flow  as  the  vor¬ 
tices  actually  run  into  it  at  the  downstream  end.  Product  measurements  were 
restricted  to  the  upstream  2/3  of  these  pictures. 

Just  below  each  of  the  forced  flow  pictures,  a  line  is  drawn  to  represent  the 
forcing  wavelength,  X  .  for  the  conditions  under  which  the  picture  was  taken. 
Wygnanski  and  Oster’s  scaling  argument  (see  Browand  and  Ho  (1983)  or  Figure 
12)  is  qualitatively  apparent.  If  the  pictures  were  enlarged  so  as  to  make  the 
wavelengths  equal,  the  three  forced  flows  would  appear  very  similar  in  the 
frequency-locked  region.  This  supports  the  hypothesis  that  the  thickness  of  the 
forced  shear  layer  scales  with  the  forcing  wavelength.  One  feature  which  would 
not  appear  similar  in  such  a  scaling  is  the  number  of  vortices  which  amalgamate 
to  make  up  the  frequency  locked  structure.  Initial  conditions  must  determine 

F 

that  number,  i.e.,  the  number  depends  on  . 
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One-second  duration  time  exposures  of  the  phenolphthalein  reaction  are 
shown  in  Figure  3.12  for  flow  conditions  identical  to  those  in  Figure  3.11.  Figure 
3.12a,  the  unforced  flow,  shows  a  more  or  less  linear  growth  as  expected  while 
Figure  3.12b  shows  very  clearly  the  three  regions  of  growth  characteristic  of  the 
forced  shear  layer. 

The  edge  of  the  mixed  fluid  is  fairly  easy  to  identify  in  these  photographs. 
Several  points  along  these  boundaries  were  selected  and  digitized  as  described 
m  Chapter  2.  The  points  selected  for  each  picture  in  Figure  3.12  are  shown  in 
Figure  3.13.  The  horizontal  and  vertical  scales  are  x  and  y  in  centimeters.  The 


line  connecting  the  points  on  each  edge  is  the  function  resulting  from  the  least 
squares  fit  described  in  Chapter  2.  The  y  distance  between  the  fitting  functions 
was  measured  and  found  to  be  roughly  the  same  as  the  visual  thickness,  6 ^  ,  of 
Brown  and  Roshko  (1974).  The  thickness  was  measured  at  sixteen  x  locations 
spaced  1  5  cm  apart. 

The  thickness  was  scaled  with  forcing  wavelength  and  plotted  against  Wygnan- 
ski  and  Oster's  non-dimensional  downstream  distance  parameter  The  result  is 
shown  in  Figure  3  14.  The  case  at  the  highest  frequency,  F  -  12  Hz ,  agrees  well 
with  Figure  1.2.  The  F  -  8  Hz  case  has  acceptable  agreement  though  it  is  not  as 
good  as  the  higher  frequency  case.  The  F  =  4  Hz  data  have  a  smaller  growth 
rate  for  small  Xwo  than  Figure  1.2  would  predict.  If  the  growth  rate  were  as 
large  as  it  “should"  be,  the  structures  would  exceed  the  vertical  extent  of  the 
test  section.  Thus,  the  test  section  restricts  the  growth  in  the  Xwo  <  1  region.  It 
is  possible  that  initial  conditions  are  responsible  for  the  sma’l  growth.  However, 
the  similarity  growth  hypothesis  for  the  forced  shear  layer  suggests  that  the 
beginning  of  the  frequency-locked  region  is  independent  of  initial  conditions. 
Both  Wygnanski's  and  Ho's  data  are  consistent  with  this  observation  even 
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though  Wygnanski  and  Oster  used  small  values  of  —  while  Ho  and  Huang 

worked  with  values  as  high  as  1.  A  recent  review  paper  by  Ho  and  Huerre  (1904) 
presents  forced  shear-layer  data  plotted  in  the  same  dimensionless  coordinates. 
There  is  a  considerable  amount  of  scatter  in  the  results.  The  data  in  Figure  3.14 
fall  within  the  range  of  values  in  Ho  and  Huerre’s  plot.  Another  reason  for  the 
variation  in  the  present  data  could  be  the  use  of  visual  thickness  as  a  measure 
of  width,  In  the  unforced  shear  layer,  the  various  measures  of  thickness  are  all 
proportional.  The  use  of  visual  thickness,  <5^,  in  the  present  case  instead  of 
integral  thickness,  9,  as  Ho  and  Huerre,  assumes,  without  justification,  that  this 
is  also  true  in  the  forced  layer. 


Figure  3.14  also  shows  that  flow  at  the  end  of  the  frequency-locked  region  for 
F  -  12  Hz  corresponds  to  conditions  at  the  neutral  point  in  a  linear  stability 
analysis  as  pointed  out  previously. 

3.3.  Product  Thickness  Measurements  by  Absorption 

In  this  section  we  present  the  basic  measurements  of  reaction  product  made 
by  the  absorption  method  described  in  Chapter  2.  The  first  part  describes  the 
format  in  which  the  measurements  will  be  presented  and  the  remainder 
describes  the  results.  While  much  of  the  material  in  the  previous  section  on 
structure  was  discovered  by  previous  investigators,  the  contents  of  the  present 
section  and  those  that  follow  represent,  to  the  best  of  our  knowledge,  the  first 
measurements  of  product  in  forced  shear  layers.  Prior  to  these,  it  was  not  clear 
whether  there  was  any  correlation  between  the  effects  of  forcing  on  growth  rate 
and  the  effects  on  product  production.  The  basic  measurement  presented  in 
this  section  is  the  streamwise  distribution  of  product  thickness 
r  tf<59 

<5p(x)  =  Jecp{y)dy  and  its  derivative,  ,  which  will  be  called  "the  product 

growth  rate".  These  give  one  measure  of  the  mixing  processes.  Another  and 
possibly  better  measure  would  be  the  product  flux  f  cp  {y  )u  (y )  dy  ;  its  deriva¬ 
tive  m  the  flow  direction  would  give  the  local  rate  of  chemical  production.  Since 
simultaneous  measurements  were  not  made  of  cp(y,t)  and  u{y,t),  it  was  not 
possible  to  determine  this  quantity,  which  can  be  written  as  the  sum  of  the  two 
terms  f^udy  +  fa  cp'u'  dy .  If  the  second  term  is  neglected  and  the  first  one 
is  written  as  Uyfip.  where  Up  is  a  suitable  mean  velocity,  then  our  product 

dip 

growth  rate,  . r  .  multiplied  by  u -,  would  be  a  measure  of  chemical  produc- 
ax 

tion  rate,  provided  Up  does  not  vary  with  x.  While  the  quantity  is  not  as 

precise  a  measure  of  chemical  production  rate,  we  believe  it  i9  a  useful  indicator 
of  local  mixing  rates. 
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Eight  data  sets  were  recorded  and  analyzed.  The  results  of  each  set  are 
presented  in  Figures  315  to  3.22  in  a  standard  format  so  that  the  reader  may 
compare  one  set  to  another  or  to  their  own  data.  The  main  result  of  the 
analysis  of  each  data  set  is  the  mean  product  thickness  as  a  function  of  down¬ 
stream  distance.  This  plot  appears  in  the  center  on  the  left  side  in  the  standard 
format.  Figure  3.15  may  be  referred  to  at  this  time  to  better  understand  this 
description.  The  circles  in  the  product  thickness  plot  are  the  actual  measured 
values  while  the  solid  line  represents  what  is  judged  to  be  a  reasonable  smooth 
fit  indicative  of  the  trend  of  the  measurements.  Above  the  product  thickness 
plot,  the  derivative  of  the  solid  line  “fit”  is  shown.  This  represents  the  rate  of 
change  of  product  thickness  with  downstream  distance  or  more  simply  stated, 
the  product  growth  rate.  Below  the  product  thickness  plot,  is  shown  a  photo¬ 
graph  of  the  flow  for  which  the  measurements  were  taken.  The  photograph  was 
enlarged  in  such  a  way  that  the  scales  agree  with  those  on  the  plots,  eg. 
x  =25  cm  in  the  plots  is  the  same  as  *  =  25  cm  in  the  photograph- and  the  ori¬ 
gin  in  the  plot  is  the  tip  of  the  splitter  plate  in  the  picture.  In  this  way.  one  can 
compare  the  product  thickness,  its  growth  rate,  and  the  visual  information 
along  a  single  vertical  line  in  this  standard  format. 

Finally,  in  the  upper  right,  the  various  values  of  important  parameters  for 
the  flow  in  question  are  presented.  At  the  top  of  the  list  are  the  three  basic 
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parameters  of  the  flow;  Ux,  the  velocity  of  the  high  speed  side;  r  =  — - ,  the  velo- 

u  i 

city  ratio;  and  F,  the  forcing  frequency.  The  remaining  entries  in  the  list  are 
derived  quantities.  e:  is  the  initial  momentum  thickness  of  the  high  speed 

boundary  layer  in  millimeters.  This  allows,  for  example,  the  value  of  —  to  be 
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quickly  estimated  for  a  given  x  location.  F0  is  the  most  unstable  frequency  from 
linear  stability  theory,  which  can  be  compared  to  the  forcing  frequency.  The 


next  two  entries  in  the  list  are  the  initial  and  unit  Reynolds  numbers  respec¬ 
tively.  The  unit  Reynolds  number,  when  multiplied  by  some  measure  of  the  local 
thickness  of  the  layer,  provides  an  estimate  of  the  large-structure  Reynolds 
number.  The  last  entry  provides  an  estimate  of  the  downstream  location  where 
the  frequency-locked  region  begins.  Recall,  from  Figure  1.2,  that  the  frequency- 
locked  region  ends  at  roughly  twice  this  value. 

Data  were  recorded  at  two  different  Reynolds  numbers.  These  will  be  referred 
to  in  the  remainder  of  this  section  as  simply  "high”  and  "low"  Reynolds  number 
cases.  In  the  high  Reynolds  number  case,  which  will  be  considered  first,  the 
unforced  flow  and  two  forced  flows  will  be  presented.  For  the  low  Reynolds 
number  case,  the  unforced  flow  and  four  forced  flows  will  be  considered.  The 
photographs  for  the  high  Reynolds  number  flows  have  the  high  speed  side  on  the 
bottom  and  an  inert  dye  added  to  the  high  speed  side.  The  low  Reynolds  number 
photographs  have  no  inert  dye  and  the  high  speed  side  is  on  top  in  the  side  view. 

The  high  Reynolds  number,  unforced  flow,  Figure  3.15,  represents  the  founda¬ 
tion  of  the  present  data  since  it  was  as  close  to  the  limiting  case  of  an  entirely 
post-mixing-transition  flow  as  the  facility  and  data  acquisition  methods  would 
allow.  The  growth  rate  of  product  thickness  is  small  for  small  x  and  increases 
monotonically  with  increasing  x.  The  derivative  plot  shows  that  the  slope  of  the 
product  curve  is  constant  beyond  z  *  9  cm  with  the  value  of  the  slope  equal  to 
approximately  0.04.  The  measurements  are  consistent  with  the  information  in 
the  photograph.  Notice  that  the  beginning  of  the  constant  slope  section 
corresponds  well  with  increased  amounts  of  mixed  fluid  in  the  photograph. 

The  beginning  of  the  constant  slope  region  at  about  x  =  9  cm  in  Figure  3.15 
apparently  represents  the  beginning  of  the  post-mixing  transition  region  for  this 
flow.  The  only  data  available  for  comparison  are  the  data  of  Breidenthal  (1978) 
which  were  taken  in  the  same  facility  using  the  same  method.  The  only 


difference  is  that  Breidenthal  measured  product  at  a  single  spanwise  location 
while  the  present  measurements  gave  a  spanwise  average.  Breidenthal's  meas¬ 
urements,  shown  in  Figure  1.1,  indicate  the  mixing  transition  is  complete  at  this 
velocity  ratio  atr»  20  cm.  Hence,  a  discrepancy  by  a  factor  of  two  or  more 
exists  between  Breidenthal's  and  the  present  measurements.  The  spanwise 
averaging  in  the  present  case  might  account  for  some  of  the  difference  but  it  is 
hard  to  believe,  looking  at  the  photograph  in  Figure  3.15,  that  it  makes  a 
difference  of  a  factor  of  two.  In  the  present  work,  it  is  the  effect  of  forcing  on 
mixing  that  is  the  main  interest  and  from  this  stand  point  it  makes  sense  to 
compare  forced  flows  with  unforced  flows  measured  in  the  same  way  For  this 
reason,  completion  of  the  mixing  transition  at  x  -  9  cm  for  this  unforced  flow 
will  be  used  as  the  reference  value  for  comparison  with  forced  flows  in  what  fol¬ 
lows. 

Figure  3.16  shows  the  result  of  forcing  the  above  flow  at  F-3.9HZ.  The 
frequency-locked  region  begins  at  x  -  27.6  cm,  which  is  consistent  with  the 
visual  information  in  the  photograph.  The  product  distribution,  in  this  case,  is 
not  very  different  from  that  in  the  unforced  case  except  at  large  x,  where  a 
noticeable  increase  in  mean  product  and  production  rate  are  observed.  For 
small  values  of  x  in  this  flow,  the  forcing  causes  a  simple  "flopping’’  of  the  layer 
from  side  to  side,  which  does  not  increase  mixing.  However,  the  “flopping”  even¬ 
tually  leads  to  a  strong  roll-up  at  approximately  x  =  20  cm  which  causes  twice 
the  asymptotic  product  growth  rate  of  the  unforced  case. 

Figure  3.17  shows  one  of  the  most  interesting  data  sets  obtained  during  this 
research.  The  forcing  frequency  in  this  case  is  8  Hz  which  means  the  frequency- 
locked  region  begins  at  13.4cm  and  ends  at  26.8cm.  The  flow  goes  through  the 
mixing  transition  in  much  the  same  way  as  the  unforced  case  with  the  asymp¬ 
totic  value  of  production  rate  achieved  by  x  =7  cm  However,  the  product 
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thickness  stops  growing  at  the  begining  of  the  frequency-locked  region.  Approxi¬ 
mately  half  way  into  the  frequency-locked  region,  the  product  begins  to  grow 
again  at  a  rate  slightly  larger  than  the  asymptotic  unforced  rate.  The 
frequency-locking  phenomenon  thus  causes  a  temporary  zero  product  growth 
rate.  The  structures  apparently  convect  and  rotate  without  doing  any  mixing. 
Comparison  with  the  unforced  flow  shows  that  forcing  creates  more  product 
upstream  of  x  =  15  cm  and  less  product  downstream  of  that  point. 

This  concludes  the  presentation  of  the  high  Reynolds  number  data.  The  low 
Reynolds  number  case  will  now  be  addressed.  The  flow  speeds  for  these  data 
were  roughly  a  factor  of  four  smaller  than  in  the  previous  flows  so  the  initial 
Reynolds  number  is  half  the  value  for  the  above  flows  and  the  unit  Reynolds 
number  is  one-quarter  the  previous  value. 

Figure  3.10  shows  the  unforced  low  Reynolds-number  data.  In  this  case,  there 
is  much  less  product  than  that  measured  for  the  high-Reynolds  number  case, 
which  is  consistent  with  pre-mixing  transition  conditions.  At  approximately 
x  =  20  cm  the  production  rate  begins  to  rise  sharply,  indicating  the  beginning 
of  the  mixing  transition. 

The  same  flow,  when  forced  at  F  =  1.8 Hz  is  shown  in  Figure  3.19.  There  is 
very  little  mixing  until  the  beginning  of  the  frequency-locked  region.  This  is  an 
interesting  case  since  the  beginning  of  the  frequency-locked  region  here  coin¬ 
cides  with  the  location  at  which  the  mixing  transition  would  begin  in  the 
unforced  case.  The  product  growth  rate  and  the  product  thickness  increase 
sharply  at  that  point  for  the  forced  flow.  The  product  growth  rate  actually 
exceeds  the  post-mixing  transition  value  measured  for  the  high  Reynolds 
number  unforced  case.  This  implies  that  the  product  growth  rate  overshoots 
since  it  is  expected  that  the  unforced  asymptotic  value  will  be  approached  for 
sufficiently  large  x  in  Figure  3.19. 
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Figure  3.20  shows  the  low  Reynolds  number  flow  for  a  forcing  frequency  of 
2  6 Hz.  Here,  as  in  the  previous  forced  case,  the  product  and  the  product 
growth  rate  both  increase  rapidly  at  the  begining  of  the  frequency-locked 
region.  Unlike  the  high  Reynolds-number  forced  flows,  the  product  growth  rate 
is  Largest  in  the  frequency-locked  region  for  these  low  Reynolds-number  flows. 

Figure  3.21  for  a  forcing  frequency  of  4.5  Hz  once  again  shows  the  same 
behavior  as  seen  in  the  previous  two  examples.  The  product  and  product  growth 
rates  are  both  small  until  the  begining  of  the  frequency-locked  region.  In  this 
case,  however,  the  product  begins  to  level  off  and  the  product  growth  rate 
reduces  at  about  z  =20 cm.  This  is  slightly  downstream  of  the  end  of  the 
frequency-locked  region  which  occurs  at  z  =  16.2  cm  for  this  case 

Figure  3.22  shows  the  final  Low  Reynolds  number  data  set.  Unfortunately,  a 
photograph  for  this  forcing  frequency  of  5.9  Hz  was  not  obtained.  The  trends 
seen  in  the  previous  forced,  Low  Reynolds-number  flows  are  also  true  for  this 
case.  The  product  thickness  and  product  growth  rate  increase  at  the  begining 
of  the  frequency-locked  region.  After  this  the  product  thickness  levels-off  with 
corresponding  Low  growth  rate.  This  flow  is  particularly  interesting  since  the 
actual  product  thickness  at  z  =  12  cm  is  almost  5  times  the  unforced  value  at 
the  same  point.  With  regard  to  the  missing  photograph,  pictures  were  obtained 
at  forcing  frequencies  both  higher  and  lower  than  5,9  Hz  and  there  was  nothing 
unusual  about  them.  The  flow  picture  would  look  like  a  rescaled  version  of  the 
photograph  in  Figure  3.21. 

3.4.  Product  Thickness  Fluctuations  from  Absorption  Measurements 

In  this  section,  the  fluctuating  product  thickness  time  series  will  be  dis¬ 
cussed,  The  mean  product  thickness,  presented  in  the  previous  section,  is  the 
long-time  average  of  a  fluctuating  signal  which  represents  the  product  thickness 
as  a  function  of  time  at  a  fixed  z  location.  The  product  thickness  fluctuates,  of 


course,  due  to  the  passage  of  the  large-scale  structures  (as  well  as  smaller 
ones).  The  shape  of  the  waveforms  give  a  indication,  therefore,  of  the  distribu¬ 
tion  of  product  carried  by  the  structures. 

Figure  3.23  shows  three  product-thickness  time  histories  at  z  locations  of 
6 cm,  15cm,  and  24cm  for  the  high  Reynolds  number,  unforced  case  previously 
discussed  in  Figure  3.15.  In  the  absorption  experiments,  32  seconds  of  data 
were  recorded.  The  long-time  mean  values  were  computed  by  taking  the  mean 
for  the  entire  run.  The  time  series  shown  for  each  z  location  represent  only  a  2 
second  portion  of  the  data  recorded  at  that  location.  The  vertical  scale  in  each 
plot  is  the  product  thickness  in  centimeters.  In  addition  to  the  downstream 
location,  z,  and  mean  product  thickness,  <5/>a,  the  root  mean  square  of  product 
thickness  fluctuations,  d*.  is  included  in  each  plot.  The  signals  in  Figure  3.23 
show  that  the  presence  of  the  large-scale  structures  causes  a  large  fluctuation 
in  the  product  thickness.  Spectral  analysis  of  these  signals  showed  a  rather 
broad  peak  with  the  frequency  corresponding  to  the  maximum  decreasing  like 

as  expected.  This  change  in  frequency,  or  period,  of  the  fluctuations  is 

apparent  in  the  plots  of  Figure  3.23.  Note  also  that  the  ratio  of  rms  product 
thickness  to  mean  product  thickness  is  approximately  equal  to  0.3  in  the  post- 
mixing-transition  region. 

Figures  3.24  and  3.25  show  the  time  traces  corresponding  to  the  data  sets 
presented  in  Figures  3.16  and  3.17  respectively.  The  3 .9 Hz  forced  case,  Figure 
3.24,  shows  little  change  from  the  unforced  case  which  is  consistent  with  the 
similarity  between  the  mean  product-thickness  distributions  for  the  two  cases. 
Figure  3.25,  on  the  other  hand,  shows  a  strong  response  to  the  forcing  frequency 
of  8 Hz.  The  period  of  the  forcing  corresponds  well  with  time  between  peaks  in 
the  three  traces.  Recall  that  the  high  Reynolds-number,  8  Hz  flow  had  a  zero 
production  rate  in  the  first  part  of  the  frequency-locked  region,  (c.f.  Figure 


-  44  - 


3  17).  The  three  traces  in  Figure  3.25  correspond  to  locations  upstream,  in  the 
middle,  and  downstream  of  the  zero  growth  region  The  signals  change  very  lit¬ 
tle  through  this  region. 

Figures  3.26  through  3.30  show  examples  of  time  traces  corresponding  to 
data  sets  presented  in  Figures  3.18  through  3  22  respectively.  These  are  the  low- 
Reynolds-number  data  sets.  Note  the  change  in  the  vertical  scale  from  the  high- 
Reynolds-number  plots,  Figure  3.26,  for  the  unforced  case,  shows  trends  in 
mean  product,  rms  product,  and  passage  frequency  similar  to  those  discussed 
for  the  high  Reynolds  number  unforced  case.  The  actual  values  of  mean  and 
rms  product  are  of  course  smaller 

Figures  3.27  and  3.28,  corresponding  to  forcing  frequencies  of  1.8  Hz  and 
2.6  Hz  respectively,  show  that  large  amplitude  fluctuations  are  caused  by  forcing 
the  flow  in  this  frequency  range.  Notice  that  the  ratio  of  rms  to  mean  product  is 
almost  a  factor  of  two  larger,  at  the  x  =  2  km  location,  than  it  is  for  the  post¬ 
mixing  transition  region  of  the  high  Reynolds  number  unforced  flow,  c.f.  Figure 
3  23.  This  is  due  to  the  controlled  amalgamation  which  takes  place  in  these 
flows.  When  the  vortices  merge  at  the  beginning  of  the  frequency-locked  region, 
their  mixing  is  enhanced  (see  Figures  3.19  and  3.20)  resulting  in  rather  large 
regions  of  mixed  fluid  separated  by  large  regions  of  unmixed  fluid,  c.f.  photo¬ 
graphs  of  the  flows  in  Figures  3.19  and  3.20.  This  causes  the  large  excursions  in 
the  product  thickness  seen  in  Figures  3.27  and  3.28. 

Figure  3.29  shows  time  traces  for  the  low  Reynolds-number  flow  with  a  forcing 
frequency  of  4,6 Hz.  The  x  =  6cm  trace  demonstrates  an  interesting  signature 
of  the  large  scale  structures.  This  being  a  pre-mixing-transition  case,  the  pro¬ 
duct  resides  in  a  thin  interface  which  is  wrapped  up  in  the  large  scale  vortices. 
When  the  interface  is  tangent  to  the  laser  sheet  used  for  the  measurement,  the 
instantaneous  product  thickness  is  large  and  thus  produces  a  spike  in  the  time 


trace.  The  number  of  spikes  for  a  given  structure  is  determined  by  the  number 
of  times  the  interface  has  wrapped  around  the  vortex  at  an  given  x  location. 
Another  interesting  point  in  Figure  3.29  is  the  small  rms  value  of  product  thick¬ 
ness  at  the  15cm  location.  This  is  the  end  of  the  frequency-locked  region  for 
this  flow  and.  as  was  demonstrated  in  Figure  3.21,  a  location  where  the  rate  of 
change  of  product  is  high.  The  last  trace  in  Figure  3.29  3hows  that  the  rms  pro¬ 
duct  thickness  increases  in  the  recovery  region  (region  III)  for  this  flow. 

Figure  3.30  shows  the  low  Reynolds-number,  5.9  Hz  forced  time  traces.  The 
spikes  in  product  thickness,  mentioned  above,  occur  for  the  smallest  x  location 
and  the  other  two  traces  show  the  increase  in  rms  product  and  change  in  scale 
of  the  structures  characteristic  of  the  recovery  region. 

3.5.  Product  Measurements  by  Laser  Induced  Fluorescence 

Results  of  product  measurement  by  the  method  of  laser  induced  fluorescence 
are  presented  in  this  section.  The  absorption  measurements,  presented  above, 
are  of  integrated  values  across  the  layer.  The  product  thickness  is.  in  fact,  an 
average  of  the  y  distribution  of  product  concentration.  The  laser  induced 
fluorescence  technique,  employed  in  the  present  experiments,  allows  the  details 
of  this  distribution  in  y  to  be  measured.  In  this  sense,  the  present  section 
presents  more  detailed  information  on  the  flows  discussed  previously. 
Specifically,  it  is  the  variation  of  product  concentration  in  the  y  direction  that 
is  of  interest. 

Measurements  of  laser  induced  fluorescence  were  conducted  on  four  high- 
Reynolds-number  and  four  low-Reynolds-number  shear  layers  at  the  conditions 
described  earlier.  For  all  these  measurements,  fluorescence  of  reacted  fluid  was 
measured,  i.e.,  it  was  the  product  of  a  chemical  reaction  that  was  measured.  All 
of  the  measurements  were  made  for  a  single  downstream  location  at  x  =  25 cm.. 
For  the  high-Reynolds-number  flows,  the  unforced  flow  and  forced  flows  with 
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F  =  8  Hz  ,  12  Hz  ,  and  20  Hz  were  studied.  For  the  low-Reynolds-number  flows, 
the  unforced  flow  and  forced  flows  at  F  -  1.5  Hz  ,  2.5  Hz  ,  and  5,0  Hz  were  stu¬ 
died.  Methods  described  in  Chapter  2  were  used  to  compute  the  product  con¬ 
centration  at  256  equally  spaced  points  across  the  layer.  A  single  one  of  these 
points  is  referred  to  as  a  pixel.  The  concentration  time  histories  were  averaged 
for  the  length  of  the  run  (approximately  20  seconds  of  data)  to  obtain  the  mean 
product  concentration  profiles. 

Product  Concentration  Profiles.  The  product  concentration  profiles  for  the 
high-Reynolds-number  cases  are  shown  in  Figure  3.31.  Only  every  fifth  measure¬ 
ment  is  shown  in  this  plot,  y  is  expressed  in  this  plot  in  pixels;  one  pixel  is 
0  37  mm  so  1  cm  =a  25  pixels  in  this  plot.  The  high  speed  side  is  on  the  left  and 
the  low  speed  side  is  on  the  right. 

The  width  at  1%  of  the  maximum  concentration,  called  6j,  is  nearly  equal  to 
the  visual  thickness,  <5^,  as  defined  by  Brown  and  Roshko  (1974).  Values  of  5,. 
for  each  profile  in  Figure  3.31,  are  presented  in  Table  3  1  below  along  with 
corresponding  values  of  X TO. 


TABLE  3. 1 

Frequency 

exF 

Uc 

<5i 

6Pa 

6Pe 

<5. 

unforced 

5.2  cm 

7.2  mm 

0.14 

5  Hz 

1.9 

4.2  cm 

4.8  mm 

0.12 

12  Hz 

2.0 

4.1  cm 

5.4  mm 

0.13 

20  Hz 

4.6 

4.4  cm 

5.6  mm 

0.13 
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Returning  to  the  profiles  in  Figure  3.31,  it  is  interesting  that  the  product- 
concentration  profiles  appear  similar  in  shape  even  though  the  local  structure 
is  quite  different  for  the  different  forcing  frequencies.  The  greatest  effect  on  the 
profiles  occurs  on  the  low  speed  side  (high  pixel  number).  Recall  that  the  forc¬ 
ing  is  applied  only  to  the  high  speed  stream. 

The  product  thickness,  <5Pg  in  Table  3.1,  is  the  area  under  the  product  concen¬ 
tration  profile.  It  is  apparent  that  forcing  reduces  the  amount  of  product  at 
x  -  35  cm  in  this  shear  layer,  for  the  frequencies  measured,  which  is  consistent 
with  results  of  absorption  measurements  on  the  same  flows. 

The  mean  concentration  profiles  for  the  low-Reynolds-number  runs  are 
shown  in  Figure  3.32.  The  profiles  show  that  a  considerable  change  in  the  shape 
results  from  forcing,  in  contrast  to  the  situation  at  high  Reynolds  number  (Fig. 
3.31)  where  all  profiles  are  of  roughly  similar  shape.  The  mean  quantities, 
implied  by  the  profiles,  are  collected  in  Table  3.2  below. 

TABLE  3.2 


Frequency 

s  xF 
Uc 

Si 

Sp* 

SP2 

Si 

Unforced 

4.2  cm 

2.6  mm 

0.06 

5.0  Hz 

3.1 

3.0  cm 

3.5  mm 

0.12 

2.5  Hz 

1.6 

4.4  cm 

3.7  mm 

0.00 

1.5  Hz 

0.9 

4.9  cm 

3.7  mm 

0.06 

In  contrast  to  the  high-Reynolds-number  results,  all  of  the  low-Reynolds- 
number  forced  shear  layers  had  more  product  than  the  natural  flow.  The  5  Hz 
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forcing  produced  a  29%  decrease  in  visual  thickness  but  a  35%  increase  in  pro¬ 
duct  thickness.  At  2.5  Hz  forcing,  the  product  is  larger  by  42%  and  the  layer  is 
thicker  by  5%.  Finally  at  1.5  Hz  forcing  the  layer  is  17%  thicker  and  contains 
42%  more  product. 

Time  Traces.  The  signals,  whose  average  values  produce  the  mean  concentra¬ 
tion  profiles,  will  now  be  examined.  For  the  high-Reynolds-number  flow,  product 
concentration  as  a  fuction  of  time,  at  five  selected  pixels  on  the  array,  is  shown 
m  Figure  3.33  for  the  unforced  flow  and  in  Figure  3.34  for  the  8  Hz  forced  case. 
The  five  pixels  selected  are  spaced  20  pixels  apart  m  each  figure,  i.e.,  approxi¬ 
mately  8  mm  The  center  trace  in  each  figure  is  at  the  location  of  the  maximum 
in  the  mean  concentration  profile,  which  was  arbitrarily  choosen  to  be  the  origin 
in  y  The  value  of  y for  each  of  the  other  traces  in  both  figures  is  indicated. 
These  traces  display  information  equivalent  to  the  temperature  time  series 
obtained  by  Mungal  (1983)  in  a  gas  flow  and  are  precisely  the  same  as  the  pro¬ 
duct  time  histories  obtained  by  Koochesfahani  (1984)  in  a  liquid  flow 

The  boundaries  of  structures  are  identified  in  these  plots  by  product  concen¬ 
tration  jumps  from  relatively  low  values  (outside  the  structure)  to  a  higher 
value  (inside).  It  is  seen,  in  the  traces,  that  the  product  concentration  within 
the  structure  is  independent  of  y,  i.e.,  the  duration  of  passage  of  the  mixed 
areas  decreases  near  the  edges  of  the  layer  but  the  concentration  level  within 
the  mixed  regions  does  not  change  significantly  with  y  location.  This  is  true  for 
forced  as  well  as  unforced  flow.  Mungal  (1982)  and  Koochesfahani  (1984)  drew 
similar  conclusions  from  their  data  on  unforced  flows.  Comparing  the  unforced 
time  histories,  Figure  3.33,  to  the  8  Hz  forced  case,  Figure  3.34,  shows  that  the 
product  concentration  inside  the  structures  is  lower  for  the  forced  flow  indicat¬ 
ing  that  the  composition  of  mixed  fluid  is  altered  by  forcing.  This  point  will  be 
discussed  in  more  detail  later. 
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In  the  low-Reynolds-number  case,  time  traces  of  the  measured  product  con¬ 
centration  on  the  centerline  for  three  of  the  data  sets  presented  in  Table  3.3  are 
shown  in  Figure  3.35.  The  time  series  for  unforced  flow  demonstrates  that  mix¬ 
ing  occurs  primarily  on  large-structure  interfaces  and  that  there  is  very  little 
small-scale  structure  in  this  shear  layer.  Notice  how  the  signal  contains  many 
spikes,  each  corresponding  to  an  interface  passing  through  the  measuring  sta¬ 
tion.  Comparison  of  the  unforced  trace  with  the  traces  in  the  high-Reynolds- 
number  layer,  Figure  3.33,  shows  that  there  are  more  of  the  spikes  in  the  low- 
Reynolds-number  flow  than  in  the  high-Reynolds-number  flow  or,  more  accu¬ 
rately.  that  the  low-Reynolds-number  interfaces  are  more  resolvable  than  are 
the  high-Reynolds-number  interfaces.  The  time  trace  of  the  F  =  5  Hz  data  has 
rather  large  sections  of  time  when  no  product  is  present  on  the  centerline.  At 
first,  it  was  thought  that  these  blank  spots  occur  between  structures  and  that 
the  mixed  regions  correspond  to  the  cores  of  the  structures.  It  was  found  that 
the  opposite  is  true,  i.e.,  the  unmixed  regions  correspond  to  the  cores  of  the 
structures  and  the  mixed  regions  occur  in  the  braids.  This  fact  will  be  obvious 
later  when  the  "flow  images"  are  introduced.  The  F  =  3.5  Hz  data  in  Figure  3.35 
contain  the  same  type  of  spikes  as  those  observed  in  the  natural  case.  The 
interfaces  in  this  case  appear  thicker  than  those  in  the  unforced  case. 

Flow  fmages.  Perhaps  the  best  way  to  display  time-series  data  obtained  in 
the  experiments  is  with  the  so-called  "flow  image".  This  is  essentially  an  inver¬ 
sion  of  the  process  that  was  used  to  take  the  data.  The  data  were  obtained  from 
intensity  information  along  a  line  in  the  flow.  The  process  can  be  inverted  by 
translating  the  data  into  intensity  information  along  a  line  on  a  video  screen. 
Moreover,  since  successive  scans  were  recorded,  the  corresponding  intensity 
lines  on  the  video  screen  can  be  placed  in  proper  sequence  to  obtain  a  "picture" 
called  the  "flow  image”.  It  is  not  an  instantaneous  spatial  picture  as  one  would 
obtain  with  a  camera  in  the  laboratory  but  rather  it  is  a  y -t  diagram  of  the 


product  concentration  along  a  fixed  line  in  the  flow  as  ti  e  vortices  in  the  shear 
layer  pass  by  that  line.  The  difference  between  the  flow  image  and  a  photograph 
is  best  realized  by  considering  the  conditions  under  which  the  two  might  be 
identical.  If  the  Taylor  hypothesis  is  valid,  namely  if  the  vortices  pass  by  the 
measuring  station  in  a  “frozen"  state,  then  the  flow  image  is  the  same  as  a  spa¬ 
tial  photograph.  The  structures  are,  of  course,  rotating  and  growing  as  well  as 
translating  so  strictly  speaking  the  Taylor  hypothesis  does  not  hold  For  more 
discussion  on  the  flow  image  the  reader  is  referred  to  Koochesfaham  11984). 

Once  the  data  are  translated  into  intensity  on  a  video  screen  in  the  form  of 
the  flow  image,  a  variety  of  image  processing  tricks  can  be  performed.  One  very 
illuminating  technique  is  pseudo-color  assignment.  The  idea  here  is  to  assign 
different  colors  to  different  intensity  levels  such  that  relative  concentration  lev¬ 
els  are  easy  to  indentify.  One  could,  for  example,  assign  blue  to  the  low  intensity 
levels  and  identify  blue  areas  in  the  resulting  image  as  regions  of  low  product 
concentration. 

Flow  images  for  the  four  high-Reynolds-number  runs,  discussed  earlier,  are 
shown  in  Figure  3.36.  Side-view  photographs  of  the  same  flows  are  shown  in  Fig¬ 
ure  3.5.  Time  in  these  flow  images  increases  from  right  ta  left  and  each  image 
contains  512  consecutive  scans,  Thus,  approximately  a  one  second  portion  of 
the  data  for  each  run  is  shown.  The  appearance  of  the  large-scale  structure  in 
the  unforced  case,  Figure  3.36a,  is  very  similar  to  the  structure  at  x  =  25  cm  in 
Figure  3.5a;  i.e.,  two-thirds  of  the  distance  downstream  in  Figure  3.5a.  Figure 
3  36b  shows  the  F  =  8  Hz  flow  image  and  it  is  apparent  that  the  forcing  organ¬ 
izes  the  vortices  in  this  case.  Comparing  Figure  3.36a  and  3  36b,  it  is  seen  that 
the  unforced  flow  possesses  more  small  scaie  turbulence  than  the  F  -  8  Hz  case. 
The  braids  in  Figure  3.36b  are  thicker  than  those  typically  found  in  the  natural 


Both  flow  images  in.  the  recovery  region  (Figures  3.36c  and  3.36d)  look  very 
much  like  the  unforced  case.  The  F  -  12  Hz  case  (Figure  3.36c)  for  which 
Xm  -  2.8,  appears  to  be  at  an  intermediate  mixed  condition  between  the 
frequency-locked  structure  of  Fig.  3.36b  and  the  new,  larger-scale  structure  at 
the  unlocked  conditions  of  Fig.  3.38d.  As  regards  small  scale  motion,  the  12 Hz 
and  20  Hz  data  look  very  similar  to  those  in  the  unforced  case  but  rather 
different  from  that  in  the  locked  region. 

The  intensity  transformation  table,  which  displays  the  color  assignments  for 
Figure  3.36,  is  shown  in  Figure  3.37.  Note  that  low  product  concentrations  are 
labeled  blue,  medium  levels  are  green,  and  high  levels  are  red.  Comparing  the 
unforced  flow,  Figure  3.38a,  with  the  8  Hz  forced  case.  Figure  3.36b,  indicates 
that  the  forced  flow  contains  more  low  concentration  areas  than  the  natural 
flow  (i.e.,  Figure  3  36b  has  more  blue  and  green  area  than  Figure  3.36a).  It  is 
obvious,  in  the  natural  case,  that  the  composition  varies  from  vortex  to  vortex. 
One  structure  in  Figure  3.36a  contains  a  considerable  amount  of  red  while 
another,  at  the  extreme  right,  has  almost  no  red  and  little  yellow.  The  8  Hz  case 
does  not  show  as  strong  a  variation;  each  vo  tex  looks  more  or  less  like  any 
other  The  20  Hz  image.  Figure  3.36d,  also  shows  little  variation  from  structure 
to  structure  but  has  more  red  areas  than  the  8  Hz  case.  Notice  the  thin  braids 
in  Figure  3.36d  compared  to  those  in  Figure  3.36b. 

Monocolor  flow  images,  in  which  grey  level  is  proportional  to  concentration, 
for  the  low-Reynolds-n umber  data  are  shown  in  Figure  3.38.  Comparing  these 
images  with  the  high-Reynolds-number  images,  it  is  seen  how  mixing  takes  place 
on  a  contorted  interface  in  the  present  case.  Comparing  the  2.5  Hz  and  1.5  Hz 
data  (Figure  3.38c  and  3.38d)  with  the  unforced  flow  (Figure  3.38a)  suggests  that 
the  increased  mixing  is  due  to  an  increase  in  the  amount  of  interface.  Both 
these  flows  apparently  benefit,  as  regards  mixing,  from  the  enhanced  vortex 
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activity  upstream  of  the  measuring  station  which  is  responsible  for  more 
stretching  of  the  mixing  interface  In  Figure  3.38c  the  Wygnanski-Oster  parame¬ 
ter  is  Xw0  =  1.8,  i.e.,  the  data  were  acquired  in  the  middle  of  the  frequency- 
locked  region.  For  1.5  Hz  data  (Figure  3.38d)  X„0  =  0.9  and  a  pairing  can  be 
observed.  At  5  Hz  forcing  (Figure  3.38b)  the  flow  has  a  very  unusual  structure. 
The  Wygnanski-Oster  parameter  for  this  case  is  3.1  so  we  expect  these  data  to  be 
m  the  recovery  region.  The  vortices  in  this  case  contain  hollow  cores,  as  indi¬ 
cated  previously,  with  no  mixed  fluid  in  them  while  the  braids  are  greatly  accen¬ 
tuated  and  contain  most  of  the  product. 

Concentration  ftistagrams.  The  concentration  histogram  is  another  method 
of  displaying  the  information  from  the  experiments.  It  is  a  plot  of  probability 
versus  concentration  and  cross-stream  coordinate,  y\  i.e.,  the  concentration 
histogram  is  a  surface  in  a  three-dimensional  space  in  which  height  is  probabil¬ 
ity,  width  is  concentration,  and  depth  is  y  location.  Such  a  surface  is  shown  for 
the  unforced,  high-Reynolds-number  case  in  Figure  3.39  and  for  the  BHz  forced 
case  in  Figure  3.40.  The  unforced  concentration  histogram  has  little  variation 
in  the  y  direction.  The  forced  case,  Figure  3.40,  also  shows  little  variation  in  the 
y  direction  in  the  central  part  of  the  layer  where  the  vortex  cores  exist.  The 
edges,  on  the  other  hand,  do  have  y  variation  which  is  consistent  with  the  ear¬ 
lier  observation  that  the  braids  contain  a  larger  amount  of  low  concentration 
fluid  than  does  the  unforced  layer. 

A  less  complicated  form  of  presentation  is  obtained  by  integrating  in  the  y 
direction  to  produce  a  plot  of  probabilty  versus  concentration.  Figure  3.41 
shows  the  y  -averaged  concentration  histogram  for  each  of  the  high-Reynolds- 
number  data  sets.  Averaging  in  the  y  direction  produces  a  “total  histogram". 
The  probability  is  expressed  in  thousands  of  occurences  (25  on  the  plot  means 
25,000  realizations)  and  concentration  has  been  normalized  so  that  pure  high 


speed  fluid  is  zero  and  pure  low  speed  fluid  is  one.  i.e.,  the  concentration  axis  is 


— 2 — .  Comparison  of  Figure  3.41a,  the  unforced  case,  with  Figure  3.41b,  the 

8  Hz  forced  case,  shows  a  shift  to  the  left  in  the  forced  flow.  The  probability 
peak  in  the  natural  case  is  much  less  prominent  in  the  8  Hz  data.  This  forced 
flow,  Figure  3.41b,  shows  a  definite  increase  in  the  amount  of  fluid  with  normal¬ 
ized  concentration  between  0.0  and  0.15.  This  is  consistent  with  the  increased 
amount  of  low  concentration  fluid  in  the  braids  for  the  forced  flow  as  discussed 
earlier.  The  histograms  for  the  other  forced  cases,  Figure  3.41c  and  Figure 
3.41  d,  differ  only  slightly  from  the  natural  concentration  histogram.  The  peak 
in  Figure  3.41c  is  not  a9  strong  as  in  the  natural  case  while  the  one  in  Figure 
3  41d,  and  the  entire  histogram  for  that  matter,  appears  very  similar  to  that  for 
the  natural  flow. 

The  y  -averaged  concentration  histograms  for  the  low-Reynolds-number  flows 
are  shown  in  Figure  3.42.  The  peaks  observed  at  high  Reynolds  numbers  do  not 
occur  for  these  data.  It  is  interesting  that  all  four  concentration  histograms  in 
Figure  3.42  look  almost  the  same  even  though  the  structure,  evident  in  the  flow 
images,  is  so  different. 

3.6.  Summary  of  Forcing  Effects  on  Mixing 

The  measurements  of  product  thickness  obtained  by  the  absorption  method 
(section  3.3)  and  from  integration  of  the  laser  induced  fluorescence  concentra¬ 
tion  profiles  (section  3.5)  are  collected  in  what  may  be  called  “enhancement 
plots”  in  Figures  3.43a  and  3.44a.  The  effects  of  forcing  on  reaction  product  are 
demonstrated  by  plotting  the  ratio  of  forced  product  thickness  to  natural  pro¬ 
duct  thickness  as  a  function  of  dimensionless  downstream  distance,  Xwo.  Simi¬ 
lar  plots  for  the  ratio  of  product  growth  rates  are  given  in  Figures  3.43b  and 
3.44b  In  these  plots  a  value  greater  than  one  indicates  that  forcing  causes 


enhancement  in  either  product  thickness  or  product  growth  rate. 

The  high-Reynolds-number  enhancement  plots  are  shown  in  Figures  3  43a 
and  3  43b  for  product  and  product  growth  rate  respectively.  Recall  that  points 
which  fall  above  the  dotted  line  in  Figures  3.43  represent  flows  in  which  product 
was  enhanced  by  forcing.  It  is  clear  in  Figure  3.43  that  product  is  increased  in 
the  region  below  Xwa  =  1.0  ti.e.,  the  enhanced  growth  region).  The  forced- 
growth  plot  of  Figure  1.2  suggests  that  the  shear  layer  growth  rate  is  twice  as 
large  in  this  region  so  it  might  be  expected  that  twice  as  much  product  would 
result.  The  data  show  values  even  greater  than  two.  This  is  due  in  part  to  the 
fact  that  the  enhanced  growth  region  in  this  particular  flow  is  undergoing  the 
mixing  transition  so  the  natural  flow  has  a  small  product  thickness.  The  low- 
frequency,  high-Reynolds-number  data  are  suspect  because  the  visual  growth 
rate  was  smaller  than  what  it  should  have  been  indicating  that  the  top  and  bot¬ 
tom  of  the  test  section  may  have  been  limiting  growth.  Figure  3.43b  shows  that 
the  production  rate  is  enhanced  in  the  enhanced  growth  region  with  values 
greater  than  two  near  the  splitter  plate.  More  important,  however,  is  the  result 
that  forcing  stops  production  entirely  in  the  frequency-locked  region 

Broadly  speaking,  the  results  show  that  the  frequency-locking  phenomenon  in 
forced  shear  layers  serves  to  inhibit  mixing  in  high  Reynolds-number  flows. 

The  enhancement  plots  for  the  low-Reynolds  number  data  are  shown  in  Fig¬ 
ures  3.44a  and  3.44b  for  product  and  product  growth  rate  respectively.  The 
trend  in  these  plots  is  quite  different  from  the  high-Reynolds-number  case. 
Both  product  and  product  growth  rate  are  enhanced  in  the  frequency-locked 
region  The  reason  for  this  is  not  entirely  clear  but,  since  enhancement  by  the 
the  primary  structures  is  suppressed,  it  must  be  connected  with  the 
modification  of  the  secondary  structure.  The  secondary  structure  is  just  begin- 
ing  to  form  in  these  flows  at  the  downstream  location  where  the  frequency- 
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locking  occurs.  Photographs  of  the  flow  (Figure  3.8)  show  that  the  secondary 
structure  is  influenced  by  the  forcing.  However,  the  connection  between  details 
of  the  secondary  structure  and  improved  mixing  are  not  well  understood. 

To  summarize,  in  contrast  to  the  high-Reynolds-number  case,  the  frequency¬ 
locking  phenomenon  produced  more  mixing  in  the  low-Reynolds  number  flows. 
The  results  indicate  that  this  is  caused  by  the  effect  of  forcing  on  the  secondary 
structure  rather  than  the  primary,  two-dimensional  structure. 

3.7.  Discussion 

The  dramatic  effect  of  periodic  forcing  on  the  downstream  development  of 
the  thickness  of  the  shear  layer  is  demonstrated  by  the  plot  of  Browand  and  Ho, 
shown  in  Figure  1.2.  In  this  plot,  both  growth  rate  and  downstream  distance  are 
made  dimensionless  with  the  forcing  frequency,  i.e.,  eF/Uc  is  taken  to  be  a 
function  of  xaF/Uz  a  Xn  .  That  plot  is  influenced  heavily  by  the  original 
results  of  Oster  and  Wygnanski,  obtained  at  high  values  of  unit  Reynolds 

p 

number,  i.e.,  for  —  «  1  and  the  mixing  transition  very  close  to  the  splitter 

plate.  In  general,  however,  the  Reynolds  number  (through  its  effect  on  initial 

F 

thickness),  as  well  as  the  dimensionless  group  — .  are  parameters  of  the  flow 

r  o 

as  discussed  earlier,  c.f.  section  1.3.  Thus  one  can  expect  that  in  the 
( 9F/UC  ,  xcF/Ue )  plane  there  will  be  a  family  of  curves,  as  in  our  Figure  3.14. 
The  asymptotic  condition  represented  by  the  Browand-Ho  plot  provides  a  useful 
reference.  In  our  experiments,  the  experiment  which  most  closely  approaches 
the  asymptotic,  high  Reynolds  number  state  and  yet  displays  a  sufficiently  large 

p 

range  of  Xwo  is  the  one  summarized  in  Figure  3.17,  with  —  =  0.13.  It  is  com- 

F o 

pared  with  the  Browand-Ho  reference  curve  in  Figure  3.14,  where  the  ordinates 

9F  SvimF 

have  been  rescaled  from  —  to  — — ,  assuming  6^  =  8  e.  Our  results  are  simi- 
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lar  to  but  lie  somewhat  above  the  Browand-Ho  curve.  Part  of  the  discrepancy 
may  be  connected  with  the  rescaling  of  the  thickness. 

One  of  the  main  goals  of  this  research  was  to  describe  the  effects  of  forcing 
on  the  distribution  of  reaction  product  along  the  layer  and  to  compare  that  with 
the  unforced  case.  It  was  found  that  these  effects  of  forcing  could  not  be 
described  simply,  due  to  the  role  of  the  additional  parameters  discussed  above. 
In  particular,  the  upper  Reynolds  number  limit  of  our  facility  made  it  difficult  to 
isolate  the  mixing  transition  effects  from  the  forcing  effects.  The  broad,  qualita¬ 
tive  results,  however,  may  be  summarized  as  follows:  for  post-mixing-transition 
shear  layers,  in  which  the  length  of  the  pre-mixing-transition  region  is  small,  ini¬ 
tial  part  of  the  flow  forcing  reduces  the  amount  of  mixing  and  of  reaction  pro¬ 
duct  formed:  for  shear  layers  which  are  predominantly  pre-mixing-transitional, 
forcing  appreciably  increases  the  amount  of  product.  In  the  post-mixing- 
transition  case,  the  reduction  of  product  must  be  due  mainly  to  reduced 
entrainment  in  the  locked-in  region  and  beyond.  In  the  pre-mixing-transition 
case,  the  increased  product  may  be  partly  due  to  the  enhanced  growth  rate  of 
the  vorticity  thickness  and  whatever  enhanced  entrainment  that  produces,  but 
this  cannot  account  for  the  increases  observed  and  we  conclude  that  there  are 
important  effects  from  changes  in  the  developing  secondary  structure. 

Coming  back  to  the  simpler,  post-mixing-transition  case,  the  relation  of 
enhanced  growth  rate  and  entrainment  to  the  product  can  be  illustrated  as  fol¬ 
lows.  Dimotakis  (19B3)  has  shown  that  the  particular  mixture  ratio  inside  the 
shear  layer,  which  depends  on  the  ratio  of  entrainment  from  the  two  sides,  is  a 
function  of  growth  rate.  In  an  unforced,  linearly  growing,  self-similar  layer  this 
mixture  ratio  will  therefore  be  invariant  with  x.  The  non-linear  growth  rate 
along  a  driven  layer  implies  varying  entrainment  ratio,  according  to  Dimotakis’ 
model,  and  thus  varying  mixture  ratio  However,  the  mixture  ratio  indicated  by 
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the  amount  of  product  at  any  downstream  station  Is  an  integrated  measure 
dependent  on  all  the  events  up  to  that  point.  Comparing  the  forced  case  of  Fig¬ 
ure  3.17.  at  F  =  BHz.  with  the  unforced  flow  at  the  same  Reynolds  number,  we 
note  that  the  interaction  region  I.  which  has  twice  the  unforced  growth  rate,  has 
an  entrainment  ratio  of  1.73,  compared  to  1.31  for  the  unforced  layer,  accord¬ 
ing  to  Dimotakis'  theory.  Further  downstream  in  the  locked-in  region  (II)  the 
growth  rate  is  zero  and  so  the  entrainment  ratio  is  1,  according  to  Dimotakis’ 
model,  but  the  actual  level  of  entrainment  is  zero.  Weighing  the  contributions 
of  these  two  regions  (which  are  of  equal  length)  to  the  product  at  the  end  of  the 
locked-in  region,  predicts  34%  less  product  than  in  the  unforced  layer.  This 
agrees  well  with  the  measured  32%  reduction. 

A  further  demonstration  of  the  effect  of  forcing  on  composition  mixture  is 
given  in  Figure  3.45;  a  photograph  from  the  "disappearing  vortex”  experiment. 
The  strength  of  acid  in  the  low  speed  stream  was  increased  such  that  a  larger 
amount  of  high  speed  fluid  was  required  before  the  mixed  fluid  appeared  red; 
i.e.,  before  the  mixed  fluid  reached  the  pH  at  which  phenolphthalein  changes 
from  transparent  to  red.  Figure  3.45  shows  photographs  of  natural  flow  and  of 
the  8  Hz  flow  for  condition  in  which  the  acid  was  strong  enough  to  require  five 
parts  high  speed  fluid  to  one  part  low  speed  fluid  for  reaction  to  occur.  Thus, 
dark  areas  in  the  photographs  are  mixed  regions  where  there  is  at  least  five 
times  as  much  high  speed  fluid  as  low  speed  fluid.  The  unforced  flow  shows  vir¬ 
tually  no  product,  indicating  that  the  composition  of  mixed  fluid  never  exceeds 
flve-to-one  in  favor  of  high-speed  fluid.  In  the  forced  flow,  on  the  other  hand, 
the  cores  are  reacted  in  the  early  part  of  the  frequency-locked  region  but  disap¬ 
pear  as  they  move  downstream.  This  experiment  was  therefore  named  the 
"disappearing  vortex"  experiment.  It  shows  that  the  vortices  start  out  with  at 
least  five  times  more  high-speed  fluid  than  low-speed  fluid  in  the  mixed  cores 
and.  furthermore,  that  the  composition  ratio  of  the  mixed  fluid  reduces  as  the 
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vortices  move  downstream. 
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Effects  of  Forcing  on  Wake  Structure  and  Hiring 


4.1.  Structure  and  Miring  in  the  Unforced  Wake 

A  wake,  in  the  context  of  these  experiments,  is  created  by  setting  the  flow- 
speeds  of  the  two  streams  equal  to  each  other.  The  shear  flow  in  the  test  sec¬ 
tion  then  reflects  the  fact  that  the  splitter  plate  has  drag  and  it  is,  therefore,  a 
wake  in  the  truest  sense.  A  Reynolds  number  for  the  wake  can  be  defined  as 


Re 


v 


(4.1) 


where 

Um-Ux=  Uz  (4.2) 


and 


Gw  —  Gj  +  0%. 


(4.3) 


ev  is  the  sum  of  the  momentum  thicknesses  of  the  boundary  layers  on  the  two 
sides  of  the  splitter  plate. 

At  very  low  speeds,  the  flow  is  steady  with  mixing  only  on  a  thin  interface 
between  the  two  fluids.  Photographs  of  the  phenolphthalein  reaction  for 
Um  =  2  cm /sec  (implying  a*  =  3.2  mm  and  Re*w  =  63)  are  indistinguishable  from 
the  picture  in  Figure  3.1a  of  the  steady  laminar  shear  layer.  A  mean  velocity 
profile  was  not  obtained  for  this  flow  but  it  is  clear  that  it  would  show  that  the 
wake  thickness  is  considerably  larger  than  the  product  thickness  as  was  the 
case  for  the  laminar  shear  layer. 


I 


An  example  of  an  unforced  wake  at  a  higher  Reynolds  number,  where  periodic 
oscillations  have  appeared  in  the  flow,  is  shown  in  Figure  4.1  for 
Um  -  12  rm/sec,  implying  ew  -  1.3mm  and  Re*w  =  180.  The  photograph  shows 
that  this  wake  is  strongly  affected  by  the  side-walls  of  the  test  section  (Figure 
4.1a).  A  wedge  shaped  contamination  region  observed  on  both  side  walls  con¬ 
tains  most  of  the  mixed  fluid.  The  cross-section  view  (Figure  4.1b),  made  with 
the  laser  induced  fluorescence  technique,  shows  that  the  structure  along  the 
center  plane  of  this  flow,  away  from  the  side-wall  regions,  looks  like  a  two- 
dimensional  Karman  vortex  street  with  mixing  on  a  thin  interface  which  is  dis¬ 
torted  by  the  vortices.  The  cusped  wave  pattern  produced  in  the  interfacial 
boundary  by  the  action  of  the  vortices  is  similar  to  those  observed  experimen¬ 
tally  and  in  a  computer  simulation  for  a  mixing  layer  (see  Figure  3.3  and  Figures 
B.l  and  B.3  in  Appendix  B).  The  mixing  mechanism  is  similar  to  that  in  the  pre¬ 
mixing  transition  shear  layers  discussed  in  Chapter  3. 

At  sufficiently  high  speeds,  the  wake  becomes  "turbulent"  as  in  the  example 
shown  in  Figure  4.2,  for  which  Um  =  50  cm /sue ,  thus  6^  =  0.64  mm  and 
Re-  -  320  and  it  is  clear  that  a  transition  has  occurred.  In  contrast  to  the 
case  at  lower  Reynolds  number,  the  wake  structure  here  is  basically  three- 
dimensional  even  at  the  largest  scale,  as  evident  from  the  plan  view  in  Figure 
4.2.  This  is  consistent  with  results  of  Cimbala  (1984)  who  found  that  the  far 
wakes  of  circular  cylinders  appear  fundamentally  three-dimensional  in  plan  view 
at  sufficiently  high  Reynolds  number.  Also  relevant  are  the  experiments  of 
Breidenthal  (1980)  who  demonstrated  the  tendency  of  large  vortical  structures 
to  be  three-dimensional  in  wakes,  c.f.  two-dimensional  in  mixing  layers. 

For  the  flow  which  comes  off  a  splitter  plate,  as  in  Figure  4.2,  there  would 
appear  to  be  less  tendency  to  form  an  initially  two-dimensional  vortex  street  as 
in  wakes  behind  bluff  bodies.  Here  the  thr'  ’tmensional  structure  appears  to 
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develop  from  the  beginning. 

If  the  initial  wave  length  and  frequency  are  simply  scaled  with  the  rules 
L  -L  1. 

Re  ~  £/  5  ~  Uz  \  U  2  \  F  ~  U/6~  U2  ,  then  the  ratios  of  initial  wave 
lengths  and  frequencies  in  the  flows  in  Figures  4.2  and  4  1  are  0  5  and  8,  respec¬ 
tively.  Comparing  the  pho'os  one  can  discern  this  change  of  scale,  suggesting 
that,  as  in  the  examples  of  Breidenthal  and  Cimbala,  the  "vortex  street”  scaling 
laws  may  be  extended  to  higher  Reynolds  number  even  when  the  spanwise  insta¬ 
bility  that  contributes  to  the  three  dimensionality  of  the  primary  structure  is 
very  strong. 

The  unforced  wake  of  a  flat  plate,  then,  has  basically  three  Reynolds-number 
regimes:  steady  laminar  flow  at  low  Reynolds  number;  “turbulent”  flow  with 
three-dimensional  large  structure  at  high  Reynolds  number;  and  a  transitional, 
non-steady,  periodic,  flow  between. 

4.2.  Effects  of  Forcing  on  Structure  of  the  Wake  —  Visual  Results 

This  section  is  devoted  to  presentation  of  the  visual  results  from  forced  wake 
experiments.  The  changes  in  structure  brought  about  by  forcing  are  discussed 
as  are  effects  on  mixing  (though  photographs  allow  only  a  qualitative  assess¬ 
ment  of  mixing).  The  results  are  presented  in  order  of  increasing  Reynolds 
number. 

Forced.  Laminar  Wake.  Figure  4.3  shows  plan,  side,  and  oblique  views  of  the 
steady  laminar  wake  subjected  to  0.5  Hz  forcing.  Like  the  laminar  shear  layer, 
the  laminar  wake  responds  to  the  two-dimensional  forcing  with  a  very  three- 
dimensional  wavy  interface.  The  wake  does  not  respond  with  as  large  an  ampli¬ 
tude  as  the  shear  layer  indicating  that  the  shear  layer  is  more  unstable  than  the 
wake  at  these  low  Reynolds  numbers.  The  finite  width  of  the  test  section  and  the 
thick  side  wall  boundary  layers  at  low  speeds  are  probably  responsible  for  the 


-  62  - 


three-dimensional  response.  It  would  be  interesting  to  repeat  these  low  Rey¬ 
nolds  number  experiments  with  a  wider  test  section  to  see  if  the  response  ls 
three-dimensional  in  the  center  of  the  test  section,  away  from  side  wall 
influence. 

Structure  in  a  Forced  Wake  at  Re# =  160.  The  Re#^  =  160  wake  showed  a 
great  variety  of  interesting  changes  in  structure  when  forced.  These  flows  were 
first  investigated,  in  the  present  work,  with  the  Laser  induced  fluorescence  tech¬ 
nique.  Photographs  of  the  flow  field  for  several  forcing  frequencies  are  shown  in 
Figure  4.4.  The  acid-quenched  fluorescein  was  slightly  fluorescent  even  in  the 
"off"  state  so  pure  fluid  from  one  side  appears  black,  pure  fluid  from  the  other 
side  appears  grey,  and  the  mixed  fluid  is  white  in  a  black-and-white  picture.  The 
laser  sheet  was  aligned  along  the  center  plane  of  the  test  section  normal  to  the 
shear  layer,  for  these  photographs. 

The  photographs  in  Figure  4.4  are  all  for  the  same  flow  speed.  Only  the  fre¬ 
quency  varies  from  one  case  to  another.  As  can  be  seen,  forcing  drastically 
changes  structure  and  mixing  in  this  flow.  The  unforced  flow  is  shown  in  Figure 
4.4a.  The  basic  structure  is  a  somewhat  tenuous  Karman  vortex  street  with  mix¬ 
ing  occurring  on  a  thin  interface.  When  the  flow  is  forced  at  a  frequency  equal 
to  the  naturally  occuring  frequency  (Fq  *  5  Hz  )  a  more  organized  version  of  the 
natural  flow  results  as  seen  in  Figure  4.4b.  Apparently,  the  forcing  creates 
stronger  vortices  since  the  interface  is  wrapped  up  more  in  the  driven  case. 
Forcing  at  less  than  the  natural  frequency,  F  =  1.7 Hz,  produces  a  peculiar 
"double  Karman  vortex  street”  as  seen  in  Figure  4.4c.  This  corresponds  to  forc¬ 
ing  at  a  subharmonic  in  the  shear  layer,  The  flow  forms  structures  at  or  near 
the  natural  frequency  with  a  "modulation”  at  the  forcing  frequency.  In  the 
shear  layer,  the  modulation  leads  to  an  amalgamation  but  the  wake  adopts  the 
double  vortex  street  configuration  without  any  tendency  toward  amalgamation 
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within  the  length  of  the  test  section. 

When  the  wake  is  forced  at  frequencies  greater  than  the  natural  frequency 
the  initial  Karman  vortex  street  which  is  formed  is  frequency-locked  to  the 
input  disturbance.  However,  at  large  downstream  distances  a  breakdown  of  the 
initial  street  is  observed,  In  Figure  4.4d,  which  is  at  a  forcing  frequency  of  6  45 
Hz,  the  breakdown  or  rearrangement  can  be  seen.  Notice  the  small  mushroom¬ 
like  blobs  of  pure  fluid  from  one  side  that  jet  out  into  the  opposite  stream  on 
each  side.  At  a  forcing  frequency  of  9.8  Hz  (Fig.  4.4e)  the  breakdown  of  the  ini¬ 
tial  street  is  accompanied  by  a  dramatic  increase  in  mixing,  which  will  be  dis¬ 
cussed  later.  Increasing  the  frequency  further  to  iu.75  H2  produced  the  photo¬ 
graph  shown  in  Figure  4.4f,  The  emergence  of  a  new,  larger  scale,  Karman  vor¬ 
tex  street  following  the  breakdown  is  clearly  visible.  Notice  that  vortices  of  like- 
sign  contain  uniformly  mixed  cores  and  those  above  the  wake  centerline  have  a 
core  concentration  different  from  those  below  the  centerline. 

Close  inspection  of  the  photographs  in  Figure  4.4  reveals  that  there  is  a 
change  in  the  spacing  ratio  of  the  initial  vortex  street  with  change  of  forcing  fre¬ 
quency.  Figure  4  5  shows  how  the  lengths  a  and  b  are  defined  for  a  Karman  vor¬ 
tex  street.  The  wavelength,  a,  was  measured  from  each  of  the  photographs  in 
Figure  4.4  and  compared  with  the  forcing  wavelength.  Um/F.  The  agreement 
was  very  good  in  every  case.  The  separation,  b ,  was  also  measured  from  each 

photograph  in  Figure  4.4  and  values  for  the  spacing  ratio,  were  calculated. 

Figure  4.5  shows  how  the  spacing  ratio  varies  with  forcing  frequency  (the 
abscissa  is  forcing  frequency  made  dimensionless  by  the  natural  frequency  in 
the  unforced  case,  F0  =  5  Hz).  The  spacing  ratio  is  close  to  0.5  for  all  forcing 
frequencies  except  the  frequency  that  causes  the  "highly  rruxed  wake"  of  Figure 
4  4e  Intrigumgly,  the  spacing  ratio  for  the  highly  mixed  wake  is  nearly  equal  to 
the  theoretical  Karman  spacing  ratio  of  0.S81. 


Simultaneous  side  and  plan  views  of  the  same  flow  using  reacting  phenol- 
phthalein  visualization  axe  shown  in  Figure  4.8  for  several  of  the  frequencies 
included  in  Figure  4.4.  It  is  obvious  in  these  pictures  that  the  side  wall  influence 
is  strong  at  this  Reynolds  number.  The  reasons  for  this  will  be  discussed  in 
detail  latter  in  connection  with  a  mode)  for  the  effect.  What  is  seen  in  Figure  4.6 
is  that  forcing  significantly  affects  the  mixing  in  this  particular  wake.  Notice 
that  the  low  frequency  forced  wakes  mix  primarily  along  the  side-walls  (which 
was  also  the  case  without  forcing  as  shown  in  Figure  4.1)  while  the  high  fre¬ 
quency  forcing  causes  mixing  primarily  in  the  center  of  the  test  section. 

Secondary  (Stream-wise )  Structure  in  Forced  Wakes.  Photographs  of  the 
phenolphthalein  reaction  are  shown  for  a  Re^  =  220  wake  in  Figure  4.7.  The 
sidewall  influence  ;s  not  as  strong  at  this  Reynolds  number  as  it  is  for 
Re*w  =  160.  The  photographs  show  that  this  wake  responds  to  forcing  by  assum¬ 
ing  a  two-dimensional  Karman  vortex  street  configuration  and  that  three- 
dimensional  secondary  instabilities  form  on  the  spanwise  vortices.  The  secon¬ 
dary  instabilities  lead  to  formation  of  streamwise  vortices  with  a  topology 
different  from  the  streamwise  vortices  in  the  shear  layer  (Bernal,  1981). 

It  was  found  that  larger  forcing  amplitude  in  a  similar  flow  (Re®  =  200) 
caused  the  formation  of  strong  streamwise  vortices.  This  wake,  when  forced  at 
an  amplitude  of  5%  of  the  free-stream  speed,  is  shown  in  Figure  4.8.  The 
unforced  flow  is  shown  for  reference  in  Figure  4.8a;  note  tht  side  wall  influence 
at  this  Reynolds  number  is  similar  to  the  Re„w  =  180  wake.  Figure  4.8b  shows 
the  flow  forced  near  the  natural  frequency  with  the  larger  amplitude.  It  is  seen, 
in  the  central  part  of  the  plan  view,  that  the  secondary  instabilities  form  strong, 
closely  spaced,  streamwise  vortices  in  this  flow.  The  secondary-to-prima-y 
wavelength  ratio  is  smaller  for  this  forced  wake  than  it  is  for  the  natural  shear 
layer,  e.g.  Bernal  (1981)  found  a  ratio  of  approximately  0.5  for  the  mixing  layer 


-  65  - 


while  the  secondary-to-pnmary  spacing  ratio  measured  from  Figure  4.5b  is  0.12, 
a  factor  of  four  smaller 

Structure  in  a  Forced  Turbulent  Wake.  Figure  4  9  shows  the  flow  of  Figure  4.2 
being  forced  at  an  amplitude  of  2%  at  various  frequencies  which  are  much  lower 
than  the  natural  initial  frequency  {Fa  =  40 Hz)  as  estimated  above.  The  wake 
reponds  strongly  to  these  low  forcing  frequencies;  in  each  case  the  forced  wave 
length  is  Um/F.  compared  to  the  natural  initial  wave  length  U^/AOHz.  The 
response  is  strongest  at  F  -  9  1  Hz  for  which,  in  the  plan  view,  the  wave  appears 
to  have  developed  into  a  vortex  street  at  the  forcing  wave  length.  The  initial 
structure  at  the  natural  frequency  F0  =  40  Hz  appears  to  be  simply  incor¬ 
porated  into  the  larger,  forced  structure,  in  a  process  analogous  to  the  "collec¬ 
tive  interaction"  described  by  Ho  for  mixing  layers  forced  at  low  frequencies. 
Comparing  with  the  forced-mixing-layer  regimes  of  Figures  3.11  to  3.14  .  the 
regimes  visible  in  these  pictures  would  appear  to  correspond  to  the  enhanced 
growth  region  rather  than  to  the  later,  locked-in  region  of  no  growth.  The 
amount  of  product  is  seen  to  be  larger  than  in  the  unforced  case  but,  unfor¬ 
tunately.  no  quantitative  measurements  of  it  were  made. 

4.3.  Side-Wall  Effects  at  Moderate  Reynolds  Number 

It  has  been  shown  that  the  side-walls  of  the  test  section  have  a  strong 
influence  on  the  Re9ifl  =  160  wake  flow  discussed  earlier  (see  Figure  4  4  and  Fig¬ 
ure  4  6).  This  flow  and  the  side-wall  effect  in  it  will  be  discussed  in  detail  in  this 
section. 

When  this  work  was  begun  it  was  hoped  that  we  would  find  we  could  increase 
mixing  by  forcing  This  was  the  case  for  the  flow  shown  in  Figure  4.4e  where  an 
obvious  increase  in  mixing  occurs.  The  same  flow,  when  viewed  with  the  phenol- 
phthalein  technique,  was  found  to  have  a  highly  three-dimensional  structure. 


Simultaneous  side  and  plan  views  of  this  case  with  the  side-wall  interaction 
causing  strong  mixing  are  shown  in  Jlgure  4.10.  The  highly  mixed  region  h,  con¬ 
centrated  in  the  center  of  the  test  section  and  the  side-wall  influence  in  the 
upstream  Karman  vortex  street  is  substantial.  Further  insight  into  the  struc¬ 
ture  of  this  flow  were  obtained  by  looking  at  sectional  views  of  the  type  Bernal 
(1901)  first  obtained.  The  laser  sheet  was  placed  perpendicular  to  the  free- 
stream  flow  direction  and  photographed  from  underneath  the  test  section  look¬ 
ing  upstream.  Ideally,  a  picture  of  a  cross-section  through  the  flow  at  a  fixed  x 
location  is  obtained.  The  Bernal  section  at  x  =  12  cm  for  the  highly  mixed  wake 
is  shown  in  Figure  4.11.  The  picture  shows  that  the  flow  at  this  downstream 
location  possesses  streamwise  vorticity  in  the  form  of  counter-rotating  vortex 
pairs. 

The  photographs  suggest  that  the  structure  of  the  highly  mixed  wake  is 
caused  by  the  interaction  of  the  Karman  vortex  street  with  the  side-walls  of  the 
test  section.  The  vortices  in  the  forced  street  convect  with  the  local  velocity  so 
their  ends  will  be  tipped  in  the  upstream  direction  by  the  velocity  defect  in  the 
side  wall  boundary  layers  as  shown  in  Figure  4.12.  The  effect  is  shown,  in  Figure 
4.12,  for  only  one  vortex.  This  tipping  of  course  also  occurs  for  shear  layers, 
where  vorticity  of  only  one  sign  exists.  But  in  the  wake,  which  has  vortices  of 
both  signs,  the  tipping  has  more  serious  effects.  It  will  set  up  counter-rotating 
vortex  pairs  in  the  streamwise  direction,  along  the  sides  of  the  test  section  walls, 
whose  mirror  images  and  mutual  induction  will  tend  to  make  them  lift  off  the 
walls.  The  vortices  are  observed  in  the  Bernal  section  (Figure  4.11).  As  the  vor¬ 
tices  propagate  toward  the  center  from  the  two  sides  they  tend  to  churn  up  the 
fluid  they  encounter.  At  the  same  time,  the  three-dimensional  instability  in  the 
central  part  of  the  flow  is  enhanced.  This  and  possibly  smaller-scale  structure 
greatly  enhances  mixing.  Note,  in  Figure  4.10,  how  the  vortices  propagate 
rapidly  at  first  but  slow  down  as  the  pair  from  the  other  side  approach  until 
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they  cease  to  propagate  any  more  in  the  z  direction. 

Figure  4.13  demonstrates  a  simplified  model  of  how  this  process  takes  place. 
(It  was  suggested  by  L.  Sigurdson  in  a  private  communication).  The  views  in  Fig¬ 
ure  4.13  are  from  downstream  looking  upstream  such  that  the  free  stream  flow 
is  perpendicular  to  and  out  of  the  plane  of  the  page,  i.e.  same  as  the  Bernal  sec¬ 
tion  view.  Figure  4.13a  represents  an  x  location  near  the  splitter  plate  where 
the  tipping  has  occured  and  the  streamwise  vortices  are  established.  The 
images  of  the  streamwise  vortices  in  the  side  wails  are  indicated.  Actually,  the 
image  system  extends  cyclically  to  infinity  in  both  directions  though  only  the 
first  set  of  images  is  shown.  The  induced  velocities  on  each  vortex  are  indicated. 
The  counter-rotating  pairs  lift  off  the  side  walls  and  propagate  toward  the  center 
of  the  test  section  until  the  equilibrium  condition  in  Figure  4.13c  is  reached.  At 
this  point,  the  configuration  is  steady.  Notice  that  the  spacing  between  the 
streamwise  vortices  in  the  equilibrium  configuration  is  twice  the  distance  from 
either  pair  to  the  side  wall  from  which  it  originated.  This  is  consistent  with  the 
observed  equilibrium  spacing  seen  in  Figure  4.10  and  in  the  Bernal  section  (Fig¬ 
ure  4. 1 1). 

An  experiment  was  conducted  to  test  the  effect  of  the  thickness  of  the  side- 
wall  boundary  layer  on  the  highly  mixed  wake.  A  false  side-wall  was  inserted  in 
the  test  section.  A  triangular  notch  was  cut  in  the  false  wall  into  which  the 
splitter  plate  fit  so  that  the  boundary  layer  on  the  false  wall  had  developed  over 
3  cm  at  the  splitter  plate  tip.  This  should  be  compared  to  the  boundary  layer  on 
the  real  side-wall  which  had  developed  over  the  entire  contraction  length  of 
approximately  00  cm.  The  resulting  flow,  forced  at  the  high  mixing  frequency 
and  viewed  with  the  phenolphthalein  reaction,  is  shown  in  Figure  4.14.  The  flow 
is  not  perfectly  symmetric  in  the  wider  passage  but  the  tipping  process  and 
inward  propagation  alluded  to  above  still  occurs.  Hence,  it  is  only  necessary 


that  the  no-slip  condition  be  present  on  the  side  walls  for  the  high-mixing  wake 
structure  to  occur.  The  thickness  of  the  side-wall  boundary  layers  is  of  secon¬ 
dary  importance. 

4.4.  Mixing  Measurements  in  the  Moderate  Reynolds  Number  Wake 

In  order  to  quantify  the  increased  mixing  in  the  highly  mixed  wake,  measure¬ 
ments  of  laser  induced  fluorescence  were  conducted.  The  data  were  taken  for  a 
non-reacting  case  in  which  fluorescein  was  used  to  mark  purr  fluid  from  one 
side.  This  is  referred  to  as  a  dilution  experiment.  Eight  different  frequencies 
were  studied  at  a  downstream  location  *  =  25  cm.  Flow  images  of  the  data  are 
shown  in  Figure  4.15.  The  color  assignment  for  the  images  labels  pure  fluid 
from  one  side  as  blue  and  pure  fluid  from  the  other  as  red.  Mixed  fluid  appears 
yellow  and  green.  Figure  4.15a  is  the  unforced  case  and  the  weak  Karman  vortex 
street  structure  is  observed,  with  mixing  only  on  a  thin  interface.  Figure  4.15b 
shows  the  subharmonic  forcing  with  its  double  Karman  street  structure.  It  is 
interesting  to  note  that  one  wavelength  in  this  flow  contains  five  vortices;  two 
sets  of  counter  rotating  pairs  and  one  extra.  The  reason  for  this  is  not  clear;  it 
is  possibly  connected  with  weak  harmonics  in  the  forcing.  The  remaining  images 
display  different  stages  of  the  frequency-locked  Karman  vortex  street  and  down¬ 
stream  breakdown  discussed  earlier.  Figure  4.15d  demonstrates  that  a  nearly 
perfect  vortex  street  can  be  produced  by  forcing  at  or  near  the  natural  fre¬ 
quency.  Figure  4.15e  shows  the  begining  of  breakdown  of  the  initial  street. 
Notice  how  "tongues”  of  pure  fluid  from  one  side  extend  into  the  opposite 
stream  and  also  take  note  of  the  jet-like  flow  between  the  vortex  rows.  Figure 
4.15f  demonstrates  the  reformation  process  in  which  small  vortices  from  the 
upstream  street  begin  to  rearrange  themselves  into  a  larger  scale  Karman  vor¬ 
tex  street.  Figure  4.15g  is  the  highly  mixed  wake  flow  discussed  earlier;  the 
increase  in  mixed  fluid  is  evident  in  the  flow  image,  Finally,  at  high  frequencies 
(see  Figure  4.15h),  the  emergence  of  a  new  street  from  the  upstream,  small- 
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scale  street  is  observed. 

The  increased  mixing  in  Figure  4.15g  was  measured  from  a  concentration  his¬ 
togram  of  mixed  fluid.  Histograms  for  the  natural  wake  and  for  the  highly 
mixed  wake  are  shown  in  Figure  4.16.  It  can  be  shown  that  the  product  thick¬ 
ness  is  the  first  moment  of  the  histogram  from  a  dilution  experiment  when  the 
free  stream  delta  functions  are  removed,  e  g.  Koochesfahani  (1984).  The  narrow 
spike  at  the  left  in  Figure  4.16a  corresponds  to  pure  fluid  from  one  side  and  the 
broader  peak  on  the  right  corresponds  to  pure  fluid  from  the  other  side.  Figure 
4.16a  shows,  therefore,  exactly  what  is  obvious  in  the  image;  namely,  there  is 
very  little  measureable  mixed  fluid  in  the  unforced  case.  Figure  4.16b,  the  histo¬ 
gram  for  the  highly  mixed  wake,  has  the  same  peaks  corresponding  to  unmixed 
fluid  but  now  has  another,  broad  peak  in  the  middle  corresponding  to  mixed 
fluid.  Not  surprisingly,  most  of  the  mixed  fluid  occurs  at  a  50-50  mixture  of 
fluid  from  the  two  sides.  Computation  of  the  first  moment  of  these  histograms 
shows  that  thei  3  is  more  than  an  order-of-magnitude  increase  in  mixing  on  the 
centerline  due  to  forcing. 

4.5.  Discussion 

First,  as  regards  the  change  of  flow  regime  with  Reynolds  number  for  the 
unforced  wakes,  it  is  interesting  to  compare  with  corresponding  developments  in 
the  wakes  of  bluff  bodies,  especially  circular  cylinders,  for  which  considerable 
information  is  available.  For  bluff  bodies,  the  reference  length  is  the  cylinder 
diameter  d  ,  which  for  the  flat  plate  is,  ideally,  zero.  Momentum  thickness  is 
related  to  d  by  CDd  -  2  0w  where  CD  is  the  drag  coefficient  of  the  cylinder. 
Thus,  if  we  base  the  cylinder  wake  Reynolds  number  on  we  might  expect 

to  observe  a  rough  correspondence  between  critical  Reynolds  numbers  in  the 


two  flows. 


For  a  circular  cylinder,  it  is  known  that  the  wake  is  steady  up  to  about 
Re*  =  40,  i.e.,  Re^  =  34  (calculated  with  =  l-7)  while  for  the  splitter  plate 
the  flow  was  still  steady  at  Re^  =  63.  For  the  circular  cylinder,  the  nonsteady 
wake  is  pure  periodic  from  Re*  =  40  to  160  (Roshko,  1954),  i.e.,  Re^  =  34  to 
100  .  then  passes  an  “irregular"  or  transitional  range,  and  settles  down  to  regu¬ 
lar  but  "turbulent"  vortex  shedding  at  about  Re*  =  300  or  Re*  =  100.  A  simi- 
lar  course  of  development  is  described  by  Sato  and  Kuriki  (1961)  for  the  wake  of 
a  thin  flat  plate  (but  they  did  not  cite  specific  Reynolds  numbers).  The  splitter- 
plate  wake  in  our  experiments  seems  to  have  settled  down  to  a  "turbulent” 
regime  (Figure  4.2)  by  Re*^  =  320.  Thus  the  Reynolds-number  range  in  which 
the  flat-plate  wake  makes  its  transition  from  steady  to  turbulent,  periodic  flow 
(63  <  Re*w  <  320)  differs  somewhat  from  that  for  the  circular  cylinder  which  is 
not  surprising,  considering  the  considerable  differences  in  initial  vorticity  distri¬ 
bution.  For  the  flat-plate  wake.  Breidenthal  determined  the  beginning  of  the 
mining  transition,  where  product  rate  begins  to  increase  dramatically,  at  about 
Re*w  =  100,  i.e.,  in  the  middle  of  the  transition  range;  it  probably  corresponds  to 

the  onset  of  "irregular”  fluctuation  in  the  cylinder  wakes.  The  end  of  the  mixing 
transition  observed  by  Breidenthal  for  the  flat-plate  wake  was  at  about 
Re*w  a  800;  thus  the  flow  in  Figure  4.2  may  still  be  undergoing  the  mixing  transi¬ 
tion  although  it  appears  to  have  a  well  developed  three  dimensional  structure. 

The  wake-flow  regimes,  then,  can  also  be  labelled  as  "steady",  "pre-mixing- 
transition"  and  "post-mixing-transition",  as  for  the  mixing-layer  flows.  In  the 
mixing  layer  the  changes  from  one  regime  to  the  next  develop  in  downstream 
distance,  which  is  also  the  direction  of  increasing  Reynolds  number.  The  wake 
Reynolds  number,  on  the  other  hand,  is  invariant  with  downstream  distance,  to 
first  approximation,  and  the  regime  changes  occur  over  the  whole  wake,  so  to 
speak,  at  the  appropriate  changes  of  the  (uniform)  wake  Reynolds  number. 


Although  local  Reynolds  number  does  not  change  with  downstream  distance 
in  the  wake,  a  rescaling  of  the  large  vortex  structure  does  occur,  as  seen  in  the 
pictures  in  Figure  4.4.  For  the  unforced  case  3hown  in  Figure  4.4a,  the  "break¬ 
down’'  or  rescaling  occurs  at  z/9w  »  300,  which  is  the  same  as  the  values 
found  by  Cimbala  for  natural,  circular-cylinder  wakes.  When  forced  at  frequen¬ 
cies  higher  than  natural  the  breakdown  occurs  earlier,  for  example,  at 
x  /Qyu  a  70  in  Figure  4.4f,  This  smaller  distance  corresponds  to  the  smaller  ini¬ 
tial  scale  induced  by  the  higher  forcing  frequency. 

To  close  this  chapter,  a  note  on  the  case  of  anomalously  high  mixing  (Figure 
4.4e)  is  in  order.  First,  it  is  clear  from  the  various  photographs  in  this  chapter 
that  the  wake  flows  are  more  strongly  influenced  by  side-wall  boundary-layer 
interference  than  are  the  shear  layers.  The  difference  is  connected  with  the  fact 
that  in  the  shear  layer  the  vortices  have  only  one  sign  while  in  the  wake  the 
presence  of  vortices  of  opposite  signs  make  possible  the  strong  side-wall  interac¬ 
tion  after  their  deflection  by  the  no-slip  condition.  This  difference  in  primary 
vortex  structure  is  also  thought  to  be  connected  with  the  greater  tendency  to 
form  three-dimensional  structure  in  the  wake  (Roshko,  1976;  Breidenthal,  1980; 
Robinson  and  Saffman,  1983).  While  the  details  of  the  mechanism  leading  to  the 
greatly  enhanced  mixing  in  the  central  part  of  the  flow,  away  from  tne  side-wall 
contamination  region,  are  not  clear,  it  seems  plausible  that  the  events  are  asso¬ 
ciated  with  the  susceptibility  of  the  wake  vortex  system  to  spanwise  instability. 
The  side  wall  effects  are  strongest  and  enhancement  of  mixing  is  greatest  for  the 
flow  in  which  the  forced  initial  vortex  spacing  has  its  smallest  value,  which  turns 
out  to  be  close  to  the  well  known  Karman  value.  Smaller  vortex  spacing  implies, 
from  the  model  in  Figure  4.13,  higher  propagation  speed  away  from  the  wall. 
Whether  this  is  the  only  effect  in  the  enhanced  mixing  or  whether  there  are  addi¬ 
tional  ones,  e  g.,  strengthening  of  the  vortices,  is  not  clear. 


From  an  engineering  stand  point,  the  highly  mixed  wake  may  be  very  useful 
in  a  real  combustion  situation.  Not  only  does  the  amount  of  mixing  increase  as 
a  result  of  forcing,  but  the  mixing  occurs  in  the  center  of  the  test  section  rather 
than  along  the  side-walls  as  it  does  in  the  same  flow  without  forcing.  This  means 
that  heat  generated  from  combustion  can  be  localized  in  the  center  of  the 
combustor  rather  than  along  the  sides  were  heat  transfer  to  the  walls  may 
degrade  performance. 
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Chapter  5 

Conclusions 

The  effects  of  a  periodic  disturbance,  applied  to  one  of  the  free  streams,  on 
large-scale  structure  and  mixing  processes  in  chemically  reacting  shear  layers 
and  wakes  were  investigated.  A  wide  range  of  Reynolds  numbers,  encompassing 
three  different  flow  regimes  was  examined.  Two  different  methods  were 
employed  to  measure  the  amount  of  chemical  product  and  thus  the  extent  of 
molecular  scale  mixing,  Absorption  by  reacted  phenolphthalein  provided  cross¬ 
stream  average  product  thickness  and  laser  induced  fluorescence  intensity  pro¬ 
vided  the  product  concentration  distribution.  These  methods,  in  addition,  pro¬ 
vided  very  effective  flow  visualization  of  the  large-scale  structures  and  of  their 
response  to  the  periodic  forcing.  Following  is  a  summary  of  the  results. 


1.  The  effects  of  periodic  forcing  on  structure  of  the  shear  layer  may  be 

f 

summarized  as  follows.  For  —  <<  1  and  Re*  large  the  forced  shear 

F  o 

layer  possesses  a  frequency-locked  region  (II)  in  which  vortices  were 
spaced  at  the  forced  wavelength  and  amalgamations  do  not  occur. 
Upstream,  a  region  (I)  of  enhanced  interaction  exists  and  downstream, 


p 

a  region  (III)  of  recovery  to  the  unforced  state  occurs.  For  ~=-~  ~  1 

F 0 

but  not  small  compared  to  1,  the  interaction  region  (I)  consists  of  con¬ 
trolled  amalgamations,  i.e.  pairings,  triplings,  etc.  ,  of  vortices  formed 
at  the  splitter  plate.  The  present  measurements  of  the  distribution  of 
thickness  in  forced  shear  layers  agree  with  results  obtained  by  previ¬ 
ous  investigators.  The  results  confirm  and,  to  a  some  degree,  unify  the 


previous  results  of  Wygnanski  and  Oster  (1980)  and  Ho  and  Huang 
(1982). 

2.  The  effects  of  periodic  forcing  on  the  distribution  of  mixing  product  in 
shear  layers  is  different  at  high  and  low  Reynolds  numbers.  At  high 
Reynolds  number,  in  post-mixing-transition  shear  layers,  the  most 
prominent  result  is  that  mixing  is  greatly  reduced  or  even  completely 
stopped  in  the  frequency-locked  region  (II).  Upstream,  in  region  I,  the 
mixing  is  enhanced.  The  net  result  is  a  reduction  of  product  at  the 
end  of  region  H.  The  results  at  low  Reynolds  number,  in  pre-mixing- 

JT 

transition  shear  layers,  are  more  varied,  depending  on  and  Re*  . 

Fa  1 

The  most  prominent  result  is  that  for  some  cases  there  is  a  large 
enhancement  of  mixing  in  the  frequency-locked  region. 

3.  More  limited  measurements  were  made  on  wake  flows,  which  were 
obtained  by  setting  equal  flow  velocities  on  the  two  sides  of  the  splitter 
plate.  An  important  difference  between  the  shear-layer  and  wake  flow 
was  discovered,  namely  the  much  greater  susceptibility  of  the  wakes  to 
test-section  side-wall  influence.  The  difference  is  attributed  to  the  fact 
that  the  wake  has  vortices  of  both  signs  which,  when  swept  backward  in 
the  sidewall  boundary  layers,  interact  strongly.  A  model  based  on  this 
picture  accounts  fairly  well  for  the  strong  sidewall  effects  observed  at 
particular  forcing  frequencies.  Less  clear  is  why  under  these  condi¬ 
tions,  the  mixing  is  greatly  increased,  by  almost  an  order  of  magni¬ 
tude,  in  the  central  part  of  the  wake,  away  from  the  side-wa.  contami¬ 
nation  region;  it  almost  certainly  is  connected  with  effects  on  the 
secondary  structure. 
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The  most  surprising  result  was  the  measured  effect  of  forcing  on  mixing 
in  the  frequency-locked  region  of  the  shear  layer  at  high  Reynolds  number. 
For  a  particular  frequency  it  was  observed  that  mixing  in  that  region  was 
complete  stopped.  Oster  and  Wygnanski  had  previously  demonstrated  that 
the  growth  of  the  shear  layer  is  stopped  and,  even  more,  that  Reynolds 
stresses  are  reversed  in  this  region,  but  it  was  not  anticipated  that  mixing 
could  be  reduced  to  zero.  While  large-scale  entrainment  is  stopped,  the 
secondary  and  possibly  smaller  turbulent  structure  still  exists  at  these  Rey¬ 
nolds  numbers  and  it  is  not  clear  why  it  makes  no  contribution  to  chemical 
mixing  product. 

In  contrast,  at  low  Reynolds  numbers,  under  what  would  be  pre-mixing- 
transition  conditions  in  the  unforced  layer,  forcing  enhances  mixing  in  the 
frequency  locked  region.  In  this  case  the  secondary  structure  must  evi¬ 
dently  be  prominent  in  contributing  to  the  mixing. 

An  equally  surprising  finding  was  the  anomalous  structure  of  forced 
wakes  under  some  conditions.  The  strong  side-wall  effects  can  be  explained 
but  the  reasons  for  the  remarkable  increase  of  mixing  away  from  the 
sidewalls  is  not  yet  clear 


a> 
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Appendix  A 

Notes  on  Laser  Absorption 

A,  1.  Span  wise  Average  of  Absorption  in  a  Sheet 

Breidenthal  (1978)  used  absorption  of  light  along  a  beam  to  measure  product 
in  a  turbulent  shear  layer.  The  measurement  was  an  average  in  the  cross¬ 
stream  direction  (i.e,,  the  y  direction)  and,  when  referenced  to  a  free-stream 
concentration,  was  expressed  as  a  length  called  the  product  thickness. 

In  the  present  investigation,  a  laser  beam  was  spread  into  a  sheet  and  passed 
through  the  test  section  perpendicular  to  the  plane  of  the  shear  layer  A  bar 
photodiode  (Quantrad  model  LL  -  4  diffused  junction  silicon  photodiode)  was 
used  as  a  detector  to  monitor  the  transmitted  light  intensity.  The  device  is  a 
solid  piece  of  semiconductor  with  a  light  sensitive  area  0.1  inches  wide  and  4 
inches  long  Ideally,  the  current  from  the  diode  is  proportional  to  the  total 
power  incident  on  the  sensitive  area.  As  such,  the  measurement  is  a  spanwise 
(i.e.,  2  direction)  as  well  as  cross-stream  average  The  relationship,  derived 
below,  between  product  thickness  and  -  hotodiode  current  is  the  basis  for  the 
measurement  method 

A2.  Relation  Between  Photodiode  Current  and  Product  Thickness 

The  exponential  relationship  between  incident  and  transmitted  power  along  a 
ray  is  the  basis  for  calculating  the  product  thickness  from  absorption.  The 
input  power  is  a  function  of  3panwise  coordinate  z  as  is  the  output  power.  Their 


ratio  is 


where  a  is  the  absorption  coefficient  and  L(z )  is  the  local  value  of  the  product 
thickness.  The  photodiode  current  is  converted  to  a  voltage  with  a  standard 
operational  amplifier  circuit.  The  output  voltage,  V,  is  proportional,  therefore, 
to  the  integral  of  (2)  across  the  sheet. 

V  =  2fjr f  P*u{z)dz  (A-2) 

where  W  is  the  width  of  the  sheet,  0  is  the  conversion  constant  for  the  current 
to  voltage  converter  with  units  of  volts/ amp,  and  fx  is  the  sensitivity  of  the  pho¬ 
todiode  in  amps/watt.  Assuming  /*<»»(*)  is  a  constant,  Equations  (A.1)  and  (A.2) 
imply 

y  s  r  a-*H,)dz  (A. 3) 

"  0 

Now,  the  spanwise  average  product  thickness,  L,  is  the  average  of  /.(*); 

r 

Z=±fL(x)dz.  (A.  4) 

A  relationship  is  desired  between  the  measured  quantity,  V,  and  the  spanwise 
average  product  thickness  L.  In  order  to  achieve  this  a  linearization  in  (A. 3)  is 
clearly  necessary.  Introduce  Ljm  as  the  largest  value  L{z)  ever  achieves  and 
non-dimensionalize  in  the  following  way: 


£•(*•)=  iifil. 

mav 


■  Z  • 

*  =  y  ■  a  =  clLt 


(A.  5) 


(A.6) 


Then  has  a  maximum  value  of  1  and  (A. 3)  becomes 

i 

V  =  0/zP*  /  . 

a 


(A.7) 


Assume  0  <  o#  «  1  and  expand  the  integrand  in  (A.7). 


V 


(A.  8) 


or  in  terms  of  L 


V 


=  l  -0*7^—+  ©(o'2) 


(A.9) 


ftfiPjn  is  the  voltage  when  no  product  is  present;  i.e. 


Vm  =  foPm 


(A.  10) 


Finally, 


(A.  11) 


L  =  j-  l-f^ta**) 

«  Yi n 

Since 

1  -  tt-=  0(a*)  (A.  12) 

i 

the  error  in  Z  Is  of  order  a*  if  the  higher  order  terms  in  (A.l  1)  are  neglected.  In 
the  experiment,  a  molar  concentration  of  3  x  10"*  gave  an  absorption  coefficient 
of  0.1  cm'1.  One  centimeter  is  a  reasonable  estimate  for  Lou  .  so  a*  would  be 
0.1  and  the  error  would  be  102S. 

A.  3.  Calibration  of  Phenolphthalein  as  an  Absorbing  Medium 

Experiments  were  conducted  to  determine  the  absorption  coefficient,  a,  for 
phenolphthalein.  Mixtures  were  prepared  according  to  the  recipe  given  by 
Breidenthal  (1976);  alcohol  and  soft  water  plus  phenolphthalein  for  one  reac¬ 
tant  and  soft  water  plus  sodium  hydroxide  for  the  other.  The  base  was  mixed  so 
as  to  produce  a  pH  of  11.7  and  the  phenolphthalein  was  mixed  to  concentrations 
of  1.0  x  10"*  .  3.33  x  10"*  ,  and  1.0  x  10~®  molar  to  test  the  effect  of  concentra¬ 
tion  on  a.  A  beam  of  laser  light  was  passed  through  a  test  cell  consisting  of 
three  1  cm  cuvette  cells.  First,  one  of  the  cells  was  Oiled  with  a  pH  11.3  mixture 
of  the  two  reactants  and  the  other  two  were  Oiled  with  clear  water.  The  output 
voltage  of  the  detector  was  recorded.  Next,  two  cells  were  Oiled  with  reacted 
fluid  and  one  with  water  to  obtain  absorption  for  a  2  cm  path  length.  Finally,  all 
three  cells  were  Oiled  with  reacted  fluid.  The  results  are  shown  in  Figure  A.l  for 
the  three  different  phenolphthalein  concentrations  tested.  Equation  (A.l)  is 
applicable  since  calibration  was  done  with  a  beam  and,  appropriately.  Figure  A.  1 
shows  the  natural  logarithm  of  the  ratio  of  output  to  input  power  as  a  function 
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of  path  length.  Apparently,  the  absorption  coefficient  is  a  linear  function  of  con¬ 
centration.  The  results  in  Figure  (A.l)  imply  the  empirical  relation 

a  =  (constant)  x  concentration 


where 

( constant )  =  3.0  x  104  cm-1 

and  it  is  understood  that  concentration  refers  to  the  molar  concentration  of 

phenolphthalein. 

• 

Another  set  of  experiments  were  conducted  to  determine  the  effect  of  base 
concentration  on  the  absorption  coefficient.  Figure  A.2  shows  the  absorption 
curves  for  three  different  base  concentrations.  The  absorption  coefficient  is 
weakly  dependent  on  base  pH  for  the  range  of  values  shown  in  this  plot.  The 
weak  dependence  demonstrated  in  Figure  A.2  was  not  a  factor  in  the  experi¬ 
ments  since  accurate  measurement  of  the  volume  of  concentrated  base  that 
was  added  to  the  soft  water  was  not  difficult.  The  volume  was  of  the  order  of  a 
liter  for  a  pH  11.7  solution.  The  range  of  pH  values  shown  in  Figure  A.2 
corresponds  to  a  factor  of  25  variation  in  the  volume  of  concentrated  base,  i.e.. 
from  1/5  liter  to  5  liters  of  concentrate.  The  pH  of  the  base  in  all  the  experi¬ 
ments  varied  only  from  11.8  to  11.8.  Each  "batch’’  of  fluids  was  calibrated  at 
run  time,  in  the  way  described  above,  and  the  resulting  value  of  a  for  that 
"batch"  was  used  to  determine  the  product  thickness. 

A.  4.  Data  Reduction 

The  first  time  the  laser  sheet  power  was  measured  with  the  bar  photodiode,  a 
substantial  60  Hz  component  was  discovered  which  was  due,  most  likely,  to  the 
laser  power  supply.  In  any  case,  it  was  removed  by  recording  both  the  input  and 
output  power  as  functions  of  time.  Since  the  product  thickness  is  proportionaL 


to  the  difference  between  input  and  output  power,  within  the  assumptions  given 
above,  the  common  60  Hz  components  cancelled  each  other. 

Data  were  digitized  with  an  8-bit  computer.  In  order  to  take  advantage  of  the 
full  dynamic  range  of  the  computer,  the  D.C.  levels  were  subtracted  out  and  gain 
was  applied  to  the  remainder.  The  gains  and  D.C.  levels  were  recorded  for  each 
run  as  was  the  absorption  coefficient  for  the  fluids  used.  Later,  with  the  com¬ 
puter,  the  original  signals  were  recovered  by  applying  the  proper  (inverse)  gain 
and  restoring  the  D.C.  levels.  The  spanwise  average  product  thickness  as  a  func¬ 
tion  of  time  was  computed  on  a  point  by  point  basis  using  equation  (A.  11)  with 
higher  order  terms  neglected. 
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Appendix  B 

Roll-Up  of  a  Passive  Interface  by  Rows  of  Point  Vortices 

B.l.  Introduction 

The  motivation  for  the  model  presented  here  was,  for  the  most  part,  the  work 
of  Jimenez  (1980).  He  presented  models  for  the  visual  growth  of  the  turbulent 
shear  layer  which  bare  a  striking  resemblence  to  photographs  of  pre-mixing- 
transition  shear  layers  taken  during  the  course  of  the  present  study,  e.g.,  the 
photo  of  the  forced  flow  in  Figure  3.3  . 

One  of  Jimenez's  models  concerns  the  roll-up  of  a  thin  vortex  sheet.  An  ini¬ 
tially  flat  vortex  sheet  is  disturbed  and  allowed  to  roll-up  under  the  action  of  its 
own  induced  velocity.  The  resulting  shape  of  the  vortex  sheet  is  qualitatively 
similar  to  photographs  of  a  chemically  reacting  liquid  flow;  specifically,  a  spiral 
type  roll-up  occurs. 

Another  of  Jimenez's  models  concerns  the  roll*up  of  a  passive  interface  by  a 
row  of  point  vortices.  In  this  model  all  of  the  vorticity,  within  a  given  wavelength 
in  the  initial  sheet,  is  concentrated  into  a  single  point  vortex.  The  rest  of  the 
sheet  contains  no  vorticity  but  still  represents  a  material  surface  separating 
fluid  above  from  fluid  below.  The  interface  moves  with  the  fluid  and,  as  such, 
will  roll-up  under  the  action  of  a  row  of  point  vortices  with  the  given  spacing. 

Surprisingly,  the  two  models  produce  similar  interface  even  though  physically 
the  models  represent  opposite  ends  of  the  spectrum,  so  to  speak.  All  that 
matters,  apparently,  is  the  net  vorticity,  not  how  it  is  distributed;  at  least  as  far 
as  the  distortion  of  a  material  surface  between  high  and  low  speed  fluid  is  con¬ 
cerned.  The  advantage  of  the  point  vortex  model  is  that  it  is  a  significantly 


simpler  problem  to  solve.  For  this  reason,  we  shall  investigate  Jimenez’s  point 
vortex  model  more  deeply. 

When  examined  in  detail,  several  features  of  the  passive  interface  solution  do 
not  agree  with  experimentally  observed  flow  structures.  There  is  difficulty  with 
an  interface  that  intersects  a  point  vortex  since  the  velocity  singularity  at  the 
vortex  center  causes  infinite  distortion.  Experimentally,  it  is  known  that  the 
centers  of  the  large-scale  structures  contain  pure  high  speed  fluid  for  pre- 
mixing-transition  flow  and  that  there  exists  a  decided  asymmetry  in  the  spiral 
roll-up.  Figure  3.3  shows  both  features  quite  clearly;  the  centers  of  the  stuc- 
tures  are  pure  black  and  the  roll-up  is  not  symmetric.  Jimenez’s  point  vortex 
model  is  entirely  symmetric  with  equal  amounts  of  high-speed  and  low-speed 
fluid  within  the  spiral. 

To  alleviate  these  inconsistencies,  it  is  presently  proposed  that  a  passive 
interface  which  does  not  intersect  the  vortices,  but  rather,  one  that  is  displaced 
intially  toward  the  low  speed  side  will  better  agree  with  Figure  3.3  .  The  infinite 
distortion  is.  thus,  obviously  eliminated  and  it  is  likely  that  cores  of  pure  high 
speed  fluid  will  result.  The  problem  is  with  justification  of  such  a  model.  It  is 
hard  enough  to  accept  concentration  of  all  the  vorticity  in  a  given  wavelength 
into  a  single  point  on  the  interface  but  concentrating  it  to  a  point  off  the  inter¬ 
face  seems  even  more  removed  from  reality.  A  simple  justification  is  that  the 
vorticity  of  the  "correct”  sign  comes  from  the  high  speed  boundary  layer.  It 
seems  locigal,  therefore,  that  the  emerging  vortices  should  exist  above  the 
material  line  originating  from  the  splitter  plate. 

Koochesfahani  (1984)  made  an  observation,  based  on  linear  stability  theory 
for  the  spatially  growing  layer,  about  the  roll-up  process.  He  found  that  the  vor¬ 
ticity  contour  that  is  spaced  20%  of  the  boundary  layer  thickness  toward  the  low 
speed  side  resembles  the  experimentally  observed  interface  quite  closely.  This 


-  84- 


confirms  the  supposition  that  the  center  of  vorticity  is  above  the  material  line 
separating  the  two  fluids.  From  a  physical  standpoint,  linear  stability  theory 
preserves  more  of  the  physics  of  the  problem  than  our  model  but  with  added 
compuational  complexity.  Linear  stability  theory  is  only  valid  for  small  distur¬ 
bances  to  a  basically  parallel  flow  so  it  is  not  capable  of  predicting  what  happens 
at  downstream  distances  where  substanial  vorticity  concentration  has  taken 
place.  For  that  regime,  some  simplifications  are  necessary.  The  present  model 
retains  the  basic  elements  suggested  by  Koochesfahani's  observation  without 
the  computational  difficulty. 

Cantwell  (1981)  found  a  difference  between  vorticity  maxima  and  particle 
path  spiral  singularities.  His  study  of  transition  in  the  axisymmetric  jet  showed 
that  an  off-axis  spiral  point  develops  even  though  the  vorticity  is  a  maximum  on 
the  axis  of  symmetry. 


B.2.  Statement  of  the  Problem 

The  fluid  motion  is  assumed  two-dimensional,  inviscid,  and  incompressible 
and  as  such  satisfies  the  Euler  equations: 


du 

3* 


+ 


(B.l) 


?HL  +  u  OIL  a.  v  !*L  =  _  L  3£_ 
dt  dx  dy  p  dx 


(B2) 


dv  .  dv  ,  dv  1  dp 

dt  dx  dy  p  dy 


(B3) 


The  x  and  y  components  of  the  vorticity  vector  are  identically  zero  and  the  z 
component,  a,  is  given  by 
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9v  _  du 
dx  dry 


(B.4) 


When  the  vorticity  is  concentrated  in  isolated  singularities  the  complex  potential 
can  be  used  to  obtain  solutions  to  the  system  (B.l)  (B.2)  (B.3)  in  the  usual  way. 
Given  the  velocity  field,  the  interface  location  is  found  by  solving  an  initial  value 
problem.  Denote  the  coordinates  of  a  point  on  the  interface  by  (£(£).*?(* ))■  This 
point  moves  with  the  fluid  and  is  therefore  also  a  fluid  particle  whose  motion  is 
given  by 

=  (B.5) 

^-  =  *(*.17.0.  (B.6) 

Initial  conditions  for  (B.5)  (B.6)  are 

C(*o)  =  $0  ■  rj{t0)  =  T70  .  (B.7) 


Two  different  types  of  interface  patterns  can  be  generated  from  solutions  of 
the  system  (B.5)  (B.8)  (B.7).  The  two  will  be  referred  to  as  temporal  and  spatial 
interfaces.  The  temporal  interface  is  the  locus  of  solution  points  for  all  possible 
values  of  •  when  t0  and  rj0  are  fixed.  Physically,  this  represents  the  distortion 
of  an  interface  which  at  t  *  t0  is  a  straight  line  y  -  ij0.  The  spatial  interface,  on 
the  other  hand,  is  the  locus  of  solution  points  for  all  values  of  t0  ,  when  (0  and  ffo 
are  fixed.  As  such,  it  is  the  streakline  originating  from  the  point 
(*.  y)  *  (Co.  tjo). 
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The  goal  in  the  present  approach  is  to  take  simple,  known  solutions  to  the 
Euler  equations  which  contain  vorticity  and  integrate  equations  (B.5)  and  (B.6) 
to  determine  the  interface  distortion. 

B.3.  Analytical  Solution  for  Stuart's  Velocity  Field 

Introduce  the  streamfunction,  if/,  such  that 

u(x,y,t)  =  gf ,  v{x,y  ,t)  =  -  .  (B.0) 


The  Euler  equations  can  then  be  written  as 

du  +  du_  _  d^_  du_  _  q 
dt  dy  dx  dx  dy 


(B.9) 


dzj/  dzj/ 
dx 2  dy 2 


+  u  =  0  . 


(BIO) 


Stuart  (1967)  found  a  solution  to  (B.9)  and  (B.10)  which  can  be  used  to  model 
the  coherent  structures  of  the  turbulent  shear  layer.  His  exact  solution  is 


V'(*  <y 't )  =  Ucy  +  — In 


cosh 


+■  A  cos 


X  —  Urn  t 


(B.U) 


&U  is  the  velocity  difference  across  the  layer  and  h  is  proportional  to  the  vortex 
spacing  The  vortices  move  with  speed  Uc  in  the  positive  x  direction.  The  non- 
dimensional  parameter,  A,  is  a  vorticity  concentration  parameter.  4=0 
corresponds  to  a  parallel  flow  with  hyperbolic  tangent  velocity  profile  and  4  =  1 
corresponds  to  the  flow  due  to  an  infinite  row  of  potential  vortices.  For  further 
properties  of  the  solution  (B.ll),  the  reader  is  referred  to  Stuart’s  paper. 


Using  (B.ll)  in  (B.8)  the  velocity  field  is  found 


sinh 

HU 

u. 

h 

2 

cosh 

v. _ ^ 

+  A  C09 

x-Uet)' 

h 

H  U 


A  Sin 


x  —  Uct 


2 

cosh 

1 L 

+■  A  cos 

N 

1 

.P 

h 

l  *  i 

Inserting  this  velocity  field  into  the  particle  path  equations  of  the  previa 
tion  yields  the  following  problem  for  the  interface  solution: 


dt 


Uc  + 


HU 

2 


sinh 

2L 

h 

cosh 

> 

H 

+  A  cos 

\-uet 

ft 

1 

h 

A  sin 


t-U't 


dt  2 

cosh 

2. 

+  A  cos 

t-U't 

h 

h 

£(*o)  =  fo  .  *?(*o)  ~  Vo  ■ 

Introduce  non-dimensional  variables  In  the  following  way: 


(B  12) 


(B.13) 


sec- 


(B.14) 


(B.15) 


(B.16) 
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x 


$-Uct 

h 


a 

.  y  = 


x0 


(o  -  C/„ 

h 


Equations  (B.14),(B.  15),  and  (B16)  then  become 

dx  _  sinhy 
dt  coshy  +■  A  cosx 


dy  _  A  sinz 
dt  coshy  +  A  cosx 


x(t  0)  =  z0  .  y(?0)  =  yo 


Dividing  equation  (B.19)  by  (B.18)  yields 


dy  _  A  sinz 
dx  sinhy 


which  implies 


sinhydy  -  A  sinxcfx  =  0 


(B.  17a) 


(B.17b) 


(B.18) 


(B.19) 


(B  20) 


(B21) 


(B.22) 


Integration  of  (B22)  gives 
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coshy  +■  A  cos*  =  constant 


(B.23) 


Initial  conditions,  (B.20),  determine  the  constant  of  integration.  Introduce  $ 
such  that 


e*  =  coshy  +  A  cos2  =  coshy0  +  A  cosS0 


(B.24) 


This  simply  shows  that  the  interface  point  moves  on  the  streamline  i>  deter¬ 
mined  by  initial  conditions.  The  function  is  independent  of  time.  From  (B.24) 
and  (B.18)  it  is  seen  that 


dx 

— ar  =  e^sinhy  . 
d  t 


(B.25) 


Since  sinhy  =  VCoshzy  -  1  ,  Equation  (B.24)  can  be  used  to  write  (B.25)  as 


=  ((1  -  A  s  cos  x)z  -  e  “**]  * 


(B.28) 


[(1  -  A  e  ~*cos  z)z  -  e  _8*j  * 


(B.27) 


Integrating  (B.27)  produces 


?-w 


-a  cos z )z  -  6: 


(B.20) 


where 
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a  -  Ae~*.  b=e~*  (B.29) 

o  and  b  are  constants.  The  solution  is  found  by  determining  x(t)  from  (B.28) 
and  then  y(t)  from  (B.24). 

Physically,  it  is  known  that  there  is  a  critical  value  of  the  non-dimensional 
streamf unction,  ip.  corresponding  to  the  boundary  of  the  so  called  ‘‘cat’s  eye”. 
This  value  is  easily  found, 

$  =  ln(l  +A)  .  (B.30) 

Outside  the  cat’s  eye  boundary  e*  >  1+4  and  the  quantity  inside  the  square 
root  in  (B.28)  is  positive  for  all  real  values  of  z.  Inside  the  cat’s  eye  e*  <  1+4 
and  the  denominator  of  the  integrand  in  (B.28)  possesses  zeros  for  certain  real 
z.  A  change  in  integration  variable  will  put  (B.28)  into  a  form  where  the  above 
behavior  is  more  apparent.  Let 

<■  =  tan  2- ,  z  =  2  arctan  (B.31) 


(B.33a 


<■  =  tan  — .  Co  =  tan  — 


a 


1-g  +6  * 

1+0+6 


P  = 


o  +6  -1 
o  -6  +1 


(B.33b' 


Inside  the  cat’s  eye 


1  -4  <  e^  <  1+4 


which  implies 


(B.34] 


1  -g  <  6  <  1  +  g  . 


(B.35) 


Hence  a  and  fi  are  positive  real  numbers.  Outside  the  cat’s  eye 

e*  >  1+4 


which  implies 


6  <  1  —  o  . 


(B.30) 


(B.37) 


Once  again  a  is  positive  real  but  now  (I  Is  positive  and  purely  imaginary. 

The  integral  in  (B.32)  is  in  a  form  that  allows  the  solution  to  be  expressed 
easily  in  terms  of  Jacobian  elliptic  functions.  The  fundamental  relation  defining 
these  functions  is  the  following: 


r»  w  *  ~  w  ~  - 
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/=/ 


(B.38) 


y ,  considered  as  a  function  of  /  and  k ,  is  expressed  as 


y  =  sn  (f  ,k ) 


(B.39) 


Two  other  functions  are  defined  in  terms  of  sn 


cn(f,k)  =  >/l  -y‘ 


dn(f.k)  =  vT^ 


Reciprocals  of  these  three  functions  are  denoted  as  follows: 

ns(f.k)  =  k)  ,  nc(f,k )  =  ^  k)  .  nd(f.k)  =  ; 


Byrd  and  Friedman  (1971)  give  a  complete  development  of  the  Jacobian  elliptic 
functions  and  of  their  properties.  The  similarity  between  (B.38)  and  the  integral 
in  (B.32)  establishes  the  connection  between  the  desired  solution  and  Jacobian 
elliptic  functions. 
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B.4.  Exact  Solution  Inside  the  Cat's  Eye  Boundary 

Since  the  streamlines  inside  the  cat’s  eye  are  closed  orbits  we  expect  2  to  be 
a  periodic  function  of  ?  in  (B.28).  The  period  of  one  revolution  around  the  vor¬ 
tex  center  can  be  shown  to  be 


_ 8  K[s) _ 

V(i  +  b  +  6)(1+b  -6)(az-t-/?e) 


(B.40) 


where 

*  = 


(B.41) 


and  /f(s)  is  the  complete  elliptic  integral  of  the  first  kind. 

Since  t  and  f  0  are  known,  determination  of  the  number  of  complete  revolu¬ 
tions,  which  have  occured,  is  possible.  Let  this  quantity  be  denoted  by  N. 


N-INT 


t  -t0 
T 


(B.43) 


The  symbol  INT(x )  was  borrowed  from  Fortran.  It  simply  means  to  truncate  x 
after  the  decimal  point.  For  example 

INT( 2.123)  =  2  , 

the  result  being  a  simple  Integer. 

The  function  x(t)  is  bounded  above  and  below  in  a  predictable  way.  Define  a 
quantity  x  *  such  that 


.h. 
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cosx* 


1-6 

a 


(B.43) 


Then 


x*  <  x(T)  <  2rr  - x* . 


(B.44) 


With  the  above  information,  the  solution  needs  to  be  separated  into  four 
different  cases  corresponding  to  initial  locations  in  one  of  four  quadrants  about 
the  vortex  center.  The  four  cases  are  as  follows: 


1)  Vo  <  0  ,  x  *  <  x0  <  rr 

(B.45a) 

2)  y0  >  0  .  x  *  <  x0  <  n 

(B.45b) 

3)  y0  >  0  ,  n  <  x0  <  2n  -  x  * 

(B.45c) 

4)  y0  <  0  ,  n  <  x0  <  2rr  -  x  *  . 

(B.45d) 

For  each  case,  it  will  turn  out.  five  sub-cases  are  needed,  so  twenty  separate  for¬ 
mulas  represent  the  complete  solution.  More  development  will  be  given  for  Case 
1  in  (B.45a).  The  other  cases  follow  in  a  similar  manner. 

Consider  Case  1.  Define  the  time  t  *  as 

*. 

f  =  fg{S)dt  (B.46) 

* 

where 

7(f)  =  [{1  -«  cosf)z  -  62]  *  .  (B.47) 


Given  the  number  of  periods,  N  from  (B.42),  we  can  define  five  separate  possibil' 
ities.  They  are 
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a)  0  <  (?-?,)-  AT  <  f* 

(B.48a) 

b)  r  <  (T  -?0)  -AT <  + 

(B.48b) 

c)  *'+  l-r<  (?  -To)  —  nt <  **  +  i-r 

(B.48c) 

d)  **+  ^-r<  (?  -?o)  -  AT  <  t#+  |-r 

(B.40d) 

«)  t'  +  2-r  <  (?  -?0)  -at<  r . 

4 

(B.40e) 

Take  Case  la  for  example.  Equation  (B.2S)  implies 

A 

t  -t0-NT  =  fg($)dt  (B.49) 

with  restrictions 

0  <  (? -T0)  -  NT  <  f  ,  y<0,  x*  <  *  <  x0  .  (B.50) 


Equation  (B.49)  implies 

■o  £ 

? -t0- NT  =  f  g{t)df  -  f  g{t)dt . 

Using  (B.46) 

•a 

t  -tz-NT  =  r  -  f  g{$)d$  . 


Using  (B.47),  (B.31).  (B.33),  and  (B.43)  the  above  can  be  written  as 
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t -t0-NT  - 


*< 

V(^  +  az)«*-/3z) 


(B  51) 


where 


o  —  (1  +  a  +  b)(l  +a  —  b)  . 


Manipulation  aimed  at  obtaining  a  form  similar  to  equation  (B.38)  will  eventually 
produce 


4"  tam 
0 


|-Va (a2 +  0*)  (*'-?  +  ?„  + AT)  . 


a 


(B.52) 


which  defines  the  solution  z(t)  in  terms  of  the  Jacobian  elliptic  function 
nc(f,k )  where  /  and  k  are  the  complicated  expressions  in  square  brackets 
above.  The  twenty  formulas  mentioned  earlier,  of  which  (B.52)  is  one,  are  col¬ 
lected  together  in  section  B.6. 

B.5.  Exact  Solution  Outside  the  Cat's  Eye  Boundary 

A  fluid  particle  that  starts  outside  the  cat’s  eye  boundary  will  remain  outside 
and  simply  move  over,  or  under,  the  vortex  centers  depending  on  whether  the 
particle  starts  above,  or  belo.w  the  boundary.  Equation  (B.2B)  may  be  written 
then,  for  the  outside  flow,  as 

5 

above  ?  -  ?0  =  f  g($)d$  ,  x  >  x0  ,  p  >  0 

*o 


(B.53a) 
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bglow:  t  -tQ  -  f  g($)dt ,  x  <  x0  ,  y  <  0  (B.53b) 


where  g(()  is  given  in  equation  (B.47). 

The  characteristic  time  for  the  outside  flow  is 

2w 

t  =  fg{S)dt  (B.54) 

o 


which  is  analogous  to  the  period,  T,  for  the  inside  flow.  Physically,  it  is  the  time 
it  takes  a  particle  to  travel  from  a  location  directly  above  one  stagnation  point 
to  a  location  directly  above  the  next  stagnation  point.  The  methods  used  previ¬ 
ously  allow  r  to  be  found. 


*nr) 

aV^l  4-a+bHl+a—  6) 


(B.55) 


where 


r . 

a 


(B.50a) 


1  — tx  +6 

X 

i 

1  —a  —b 

l+a+6  ( 

■  7  = 

l+o  -6 

(B.56b) 


Following  a  strategy  similar  to  that  used  for  the  inside  problem,  the  number  of 
vortex  centers  traversed  can  be  defined  as 
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M  =  /NT 


(B.57) 


Now  define  new  coordinates  with  the  effect  of  non-zero  M  subtracted  out. 
Denote  these  by  a  prime. 


t'  =  t  - 10  -Mr,  0  <  f'<  r 


(B.58a) 


Za  =  x-xa-2  nM  ,  0  <  x,  <  2  n  ,  y  >  0 


(B.58b) 


xi  =  z0  -  x  -  2n  M  ,  0  <  x*  <  2rr  ,  y  <  0 


(B.50c) 


The  subscripts,  a  and  6 ,  refer  to  above  and  below  the  boundary;  a  for  y  >  0 
and  b  for  y  <  0.  It  is  useful  to  define  an  origin  xOT.  Above  the  boundary.  x00  is 
taken  to  be  the  x  coordinate  of  the  first  stagnation  point  to  the  left  of  x0.  Below, 
it  is  the  first  stagnation  point  to  the  right. 


y  >  0.  x00  =  2  n  Af0  ,  M0  = 


xq  >  0 


-1+/7VT  ,  x0  <  0 

C77 


(B.59a) 


y  <  0,  Xqq  =  2  tt  Mq  ,  A/0  = 


i+zm'  .  x0>  0 


^  <  0 


(B.59b) 


vr*. 


:i 


r 


a 


L 


■  7  j"  ■  r;  i-j  tf'. 1 7-;  r  WJ  ?*7  ',WJ 


■  *  "  i  ■  »  ■  i  ■  »  ' 


r  1 ,'  W  .  *  .  ■  ^ 
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then 


y  >  0,  0  <  x0  -  x00  <  2* 


y  <  0,  0  <  x00  -  xc  <  2rr 


(B.60a) 

(B.SOb) 


Consider  y  >  0;  y  <  0  will  be  very  similar.  Equation  (B.53)  can  be  written,  using 
(B.50)  and  (B.59),  as 

A  .  .  «  A  ,  ,  A  A 

*■♦*0  *.+*0-*00 

*'=  f  =  f  g(S)dt 

•o  *o -*oo 


*«+*o~*oo  *o~*oo 

t'  ~  f  ~  f  0(()d( 

o  o 


(B.61) 


The  inequalities  in  (B.50)  and  (B.60)  imply 


0  <  x,  +  x0  -  x00  <  4rr 


0  <  x0  -  x00  <  2rr 


(B.62a) 

(B.02b) 


Define  a  time,  tq  ,  as  the  time  for  a  particle  to  travel  from  x oo  to  xq  ■  Then 


*o'«oo 

Tos  f 


Putting  r0  into  the  elliptic  function  form,  using  (B.31),  produces 


(B.63) 
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4)  ~t  <  ?'  +  t0  <  2t,  3n  <  Xg+£0-xoo  <  4n 

Gt 


(B.67d) 


Consider  Case  1.  Straightforward  manipulation  of  (B.66)  gives 


t '  +  r0  = 


2F{<p\.r) 

aVM +a+6)(l+a— i 


(B.68) 


where 


.  _  1  .  *;  +  x0-5oo 

tan  59  a  =  —  tan  - - - 


(B.69) 


Using  properties  of  elliptic  functions.  (B.68)  and  (B.89)  Imply 


[*a  +  Xo— *oo  c(f')  [  1  ]» 


7  [  1  -ca(?') 


(B.70) 


where 


: (?’)  =  cn  S-vli+a-t-b^l+a-b)  (?’  +r0)  ,  r 


(B.71) 


Equation  (B.70)  defines  z*  and  hence  z  as  a  function  of  f ’  and  hence  of  t .  x(t) 
is  the  desired  solution.  The  other  three  cases  follow  in  similar  fashion  as  does 
the  solution  below  the  cat's  eye  boundary. 
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B.6.  Complete  Solution 

For  the  record,  the  complete  solution  to  the  problem  (B.18).(B  19)  (B.20)  is 
given  below  x0,  y0 ,  t0,  and  t  are  given  and  x  and  y  are  desired.  The  formulas 
for  x{t )  will  be  given  realizing  that  those  for  y(t)  are  obtained  from  equation 

(B.24). 

Inside  Solution: 

Case  1:  x*  <  x0  <  n  ,  y0  <  0 

la)  0  <  t  -t0-NT  <  t  *  ,  x'  <  x  <  x0  ,  y  <  0 

-r- tan  =  nc  +•  ^ .  (f  *  -?  +  ta  +  NT)  ,  s 

P  2  * 

16)  f  <  7  -t0-NT  <  t*  +  -~T  ,  x'  <  x  <  rr  ,  y  >  0 


1  „ 

—  tan 
0 


^  ^  (f  —  <q  -  NT  —  f  *) 

!C 


1  c)  f  *  +  f-r  <  ? -to -NT  <  f  *  +  — r  ,  n  <  x  <  2rr  -  x  *  ,  y  >  0 
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2c)  I -T  -t*  <  t  -T0  -NT  <  |-r  -t#  .  rr  <  *  <  Zn  -z*  ,  y  <  0 

**■ 


jj*  tan.  32^*.  =  nc  j  _S_±i_(?_?0._^r+r-i-r),  s 


2d)  ? -tQ- NT  <  T  -t* ,  r'  <  x  <  7T  ,  y  <  0 


4- tan  f-  =nc  a  +P  -+(T-t*-t  +t0  +  NT)  ,  s 

p  2  « 


2a)  T  -t‘  <  t-ta-NT  <T,x*<x<xo,y>0 


^  tan  f -  =nc  (t  -t0-NT -T  + 1*)  ,  s 

p  2  « 


Case  3:  tt  <  x0  <  2tt  -  x  ,  yQ  >  0 


3a)  0  <  t  -?0  -  ATT  <  t  *  ,  x0  <  x  <  2rr  - x*  .  y  >  0 


j *“  ^ 


nc  +?0  +  NT)  ,  s 

fC 


36)  t * 


3c)  f 


3d)  t * 


3e)  f* 


>  •*  -v-  .  -  v  •_  *v 
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<  t  -tQ -NT  <  t‘  +  jT  , 


<  x  <  2rr  -z* ,  y  <  o 


7tan 


2rt  -x 


-  TIC 


^S-u-^-NT-n 


*■  j-T  <  t  -t0-NT  <  «*  + 
4 


x*  <  x  <  tr  ,  y  <  0 


Vg2*/?8 

K 


(f*  +  jT-1+1q  +  NT)  , 


•—7’  <  ?-?0-7V7’<  f  *  +  t-7\ 
&  4 


x  <  x  <  rr  ,  y  >  0 


K 


(t-to-NT-t* 


h  —T<  t  —tg  -NT  < 


T,  rr  <  x  <  x0  ,  y  >  0 


•  tan 


2rr  -  x 


=  nc 


Vg2  +  ff2 


(i'  +  r-f+ro  +  ivr) , 


Cass  4:  n  <  x0  <  2rr  -  x*  ,  y0  <  0 
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4a)  0  <  t  -?0  -NT  <  ^T-t*  ,  n  <  x  <  x0  ,  y  <  0 

i-tan[^]=nc  j^2.(f-?0-JVr+O.  * 

46)  ~-T  -  t‘  <  t  -t0-NT  <  ,  z'  <  x  <  n  ,  y  <  0 


Jtan  f-  =nc  [  ~  ^  (jT -f -t +t0  + NT)  ,  s 


4c)  f -ta-NT  <  %-T  -t* ,  *#  <  *  <  n  .  y>  0 

jtan  =nc  — ^ t*- |-r)  ,  s 
4d)  | -T-t'  <  t  -f  0  -  AT  <  f  -  r  ,  rr  <  x  <  2rr  -x*  ,  y  > 


tan.  277  2~X-  =nc  (f -f' -f  +  ?0  +  Wr)  , 


4e)  f  -  t*  <  t-t0-NT  <  7\x0<x<2tt-x*,v<0 
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cot 


4jt  +200  -zo-r0 


=  -F 

7 


cn 


f-(2r-r-r0) 

iv 


Below  Outside  Solution:  y0  <  0  ,  y  <  0 


Case  1:0  <  f '  +r0  <  ^-r ,  0  <  xi  +xqo -z0  <  rr 


cot 

X(J  +  Xoo  —  Xo 

=  — F 

cn 

2 

7 

|(?'+r0).  r 


Case  2:  — r  <  f '  <  r,  n  <  xj  -f  x0o  -x0  <  2n 

<w 


cot 


2tT  -f  X0  -  Xft  -  Xqo 
2 


=  -  F 

7 


cn 


fir-t'-Tv) 


Case  3:  t  <  f‘  +r0  <  ~r  .  2n  <  xj  +x00  -x0  <  3tt 


cot 


A  •  .  A  A  a _ 

x#+xoo-*o~2tt 


-F 

7 


cn 


~  (t'  +  T0 -r) 


Case  4:  tj-t  <  t’+r0  <  2t  ,  3n  <  x*+x00-x0  <  4tt 


B.7.  Numerical  Solution  and  Interface  Patterns 


B.7.1.  Numerical  Mat  hods.  Equations  (B.5)(B.8)  and  (B.7)  can  be  solved 
directly  numerically  by  Euler  integration.  The  differential  equations  are  approx¬ 
imated  as  difference  equations  for  a  small  time  step,  At.  Starting  at  t0.  the  solu¬ 
tion  (£.77)  is  found  at  fo  +  Af .  Using  this  solution  as  initial  conditions  the  solu¬ 
tion  at  t0  +  2At  is  found  and  so  on  until  the  final  time  t  is  reached.  This  method 
will  work,  with  no  additional  complication,  for  any  velocity  field  even  if  exact 
integration  is  not  possible.  The  disadvantage  is  that  it  is  computationally  very 
slow  since  high  accuracy  requires  a  small  time  step.  A  PDP  11/44  computer 
took  as  long  as  5  hours  to  compute  a  typical  interface  of  interest. 

The  exact  solution  of  the  previous  section  can  be  evaluated  numerically  as  an 
alternative  to  Euler  Integration.  The  Jacobian  elliptic  functions  are  approxi¬ 
mated  by  taking  the  first  several  terms  of  their  Fourier  series  expansions,  is  The 
interface  pattern  computed  in  5  hours  by  Euler  Integration  takes  only  6  minutes 
to  compute  when  the  exact  solution  is  evaluated. 


B.7.2.  Results  of  Numerical  Evaluation  of  Exact  Solution..  Results  for  the 
spatial  interface  problem  introduced  earlier  will  be  presented  first.  Recall  that 
the  interface  in  this  case  is  the  same  as  the  streakline  that  originates  from  a 
point  (x0.  yo)-  From  the  discussion  in  Section  B.l,  our  interest  is  in  the  case 
when  x0  is  the  x  coordinate  of  the  splitter  plate  tip  and  yo  is  small  and  negative. 
For  demonstration  purposes,  suppose  U\  =  10  cm/sec,  Uz  =  5  cm/sec,  and  take 
the  forced  wavelength  to  be  n  cm.  The  interface  pattern  for  A  =  1.0  and 
y0  -  -0.55  cm  is  shown  in  Figure  B.l  a.  The  similarity  between  this  computed 
interface  and  the  observed  interface  for  turbulent  pre-mixing  transition  shear 
layers  is  quite  remarkable;  c.f.,  Figure  3.3.  As  expected,  an  origin  below  the 
splitter  plate  produces  cores  of  pure  high  speed  fluid.  A  cusp  forms  which  rolls 
around  the  vortex  center  as  It  moves  downstream.  Inside  the  cat’s  eye 
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boundary,  the  interface  forms  a  double  spiral,  the  two  spiral  arms  terminating 
at  the  cusp.  A  second  cusp  forms  along  the  low  speed  edge.  This  cusp  occurs  in 
the  outside  flow,  using  the  terminology  of  the  previous  sections.  The  formation 
of  these  cusps  can  be  understood  by  looking  at  the  temporal  roll-up  for  a 
moment.  The  streamlines,  which  are  fixed  for  all  time,  are  shown  in  Figure  B.2 
along  with  the  initial  interface.  Pomt  (a)  moves  around  on  its  streamline  while 
points  (b)  and  (c)  move  on  their  streamline.  The  period  for  (a)  to  orbit  the  vor¬ 
tex  center  is  less  than  that  for  points  (b)  and  (c)  because  (a)  is  closer  to  the 
vortex.  Point  (a),  being  at  the  mid-point  between  vortex  centers,  will  orbit  in  a 
shorter  time  than  any  other  point,  between  (d)  and  (e),  on  the  interface  and 
thus  will  be  the  cusp  point.  In  the  outside  flow,  point  (f)  takes  less  time  to  move 
from  one  vortex  to  another  than  any  other  between  (d)  and  (g).  Point  (f)  is, 
therefore,  the  cusp  point  in  the  outside  flow.  These  cusps  are  observed  in  the 
real  flow,  Figure  3.3. 

Figure  B.l  demonstrates  the  effect  of  vorticity  concentration  on  the  interface 
pattern.  Figure  B.la  is  the  solution  for  A  =  1.0,  the  potential  flow  solution,  while 
Figures  B.lb.  B.lc,  B.ld  are  the  solutions  for  smaller  values  of  A  and  hence,  for 
less  concentrated  vortices.  The  release  point  is  fixed  for  all  four  figures  at 
y0  =  -0.55  cm  Recall  that  the  streamline  pattern  does  not  change  significantly 
with  A.  The  width  of  the  cat's  eye  at  its  widest  point  decreases  with  decreasing 
A.  If  we  take 


f{x  ,y )  =  ln(coshy  +  A  cosx  ) 


(B.72) 


then  the  width  of  the  cat's  eye  divided  by  the  spacing  of  the  vortices  is 
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(4 


r 


JfL-  —  in  [(1  +  2A)  +  • 

2n  7T  1  1 


( B.73 ) 


This  function  is  given  in  the  table  below  for  various  values  of  A. 

TABLE  B.l 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

0.6807 

0.5123 

0.4537 

0.3797 

0.2480 

0.0 

The  width  decreases  monotonically  with  decreasing  A.  This  effect  is  visible  in 
Figure  B.l  by  comparing  the  width  of  the  last  vortex  on  the  right  in  Figures  B.l  a, 
B.lb,  and  B.lc.  In  Figure  B.ld  the  cat’s  eyes  are  so  narrow  that  the  y0  »  -0.55 
release  point  places  the  interface  entirely  in  the  outside  flow.  The  conclusion 
from  Figure  B.l  is  that  A  does  not  effect  the  qualitative  properties  of  the  inter¬ 
face.  As  long  as  part  of  the  interface  falls  inside  the  cat’s  eye  boundary,  the 
cusp  and  spiral  arms  will  form.  Qualitatively  speaking  then,  there  is  no 
difference  between  a  concentrated  point  vortex  and  a  vortex  with  non-zero  core 
size  as  far  as  the  shape  of  the  interface  is  concerned. 

The  effect  of  release  point  for  fixed  vorticity  concentration  is  demonstrated 
in  Figure  B.3.  When  the  release  point  is  near  the  z-axis,  as  in  Figure  B.3a,  most 
of  the  interface  is  inside  the  cat's  eye  so  the  spiral  arms  are  longer  and  wrap 
around  the  vortex  center  more  times  than  would  be  the  case  for  y0  further  from 
the  z-axis.  Notice  also  how  the  cusp  in  the  outside  flow  changes  with  y0.  I a  Fig¬ 
ure  B.l  a  it  is  much  sharper  than  in  B.ld.  The  time  for  the  cusp  to  form,  how¬ 
ever.  is  not  a  strong  function  of  release  point. 


B.7.3.  Numerical  Solution  by  Euler  Integration.  The  lack  of  qualitative 


changes  in  the  interface  pattern  with  vorticity  concentration  in  the  foregoing 
discussion  suggests  that  it  is  sufficient  to  limit  further  consideration  to  point 
vortices  where  the  analysis  is  simplified  by  complex  variable  theory.  In  fact,  this 
simplification  allows  an  extension  of  the  model  to  include  effects  of  a  vortex 
pairing  on  the  interface.  Jimenez  (1980)  presented  an  exact  solution  for  orbit¬ 
ing  rows  of  point  vortices  which  will  be  used  to  model  the  pairing  process.  This 
solution  is  briefly  reviewed  in  what  follows. 

The  complex  potential  for  two  infinite  rows  of  point  vortices  with  spacing  L  is 


■w(z)  -  In  sin  y-(z  -zt)  -t-  7^- In  sin  y-(z  -z2). 


(B?4) 


zx  and  z2  label  the  position  of  the  two  rows  while  7,  and  yz  are  the  strengths  of 
the  vortices  in  row  1  and  row  2  respectively.  zx  and  z2  are  functions  of  time 
determined  by  the  condition  that  each  vortex  moves  with  the  induced  velocity  of 
all  others.  This  implies 


dzx  iyz  n,  , 

—  =2Tcot  T*  ' 


(B  75) 


d*a  *7i  .  rr,  , 


(B.76) 


If  we  let 


f3*!-** 


(B.77) 


7  =  7i  +  7a 


(B.78) 
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±L=  12. 

dt  2 1 


cot 


(B.79) 


Separating  (B.79)  into  real  and  imaginary  parts  will  produce  a  system  of  equa¬ 
tions  very  similar  to  equations  (B.  18)(B.19)(B.20)  of  Section  B.3.  Without  loss  in 

generality,  choose  l  =  and  7  ~  TT  and  let  f  ■  Assuming  =  jz  ,  some 

manipulation  leads  to  the  following  system  of  equations. 

=  1 _ sinh  (y  -y)  +  1 _ sinh  ( y  +  n) 

dt  2  cosh  (y  -rj)  -  cos  (x  —  £)  2  cosh  (y  +rj)  -  cos  (x  +£)  k 

-  _2 _ sin  (x  -Q _ 1 _ sin_(x+_g)  8 

dt  2  cosh  (y  -rj)  -  cos  (x  -£)  2  cosh  (y  +rj)  -  cos  (x  +£)  '  ’ 


<LLS  L-  sinh  2 y 
dt  ”  2  cosh  Ztj  -  cos  2£ 


(B.B2) 


dq  _  _  1 _ sin  2g _ 

dt  2  cosh  Ztj  —  cos  2£ 


(B.B3) 


with  initial  conditions 

y(0)  =  x0,  y  (0)  =  Vq,  e(0)  =  e0.  t?(0)  =  0  .  (B.84) 


Physically  the  system  (B.80)  through  (B.84)  represents  the  following  problem. 
Start  with  an  infinite  row  of  point  vortices  with  spacing  n  and  strength  At  t 

=  0,  select  pairs  and  displace  them  from  their  equilibrium  location  toward  each 
other  by  an  amount  fo  The  vortices  will  move  since  they  are  no  longer  at  their 
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equilibrium  location.  The  developing  interface  is  described  by  the  system  (B.80) 
through  (B.84)  with  (£,77)  reflecting  the  vortex  locations  and  ( x  ,y )  the  interface 
location.  The  system  of  equations  is  easily  solved  by  Euler  integration. 

As  a  check,  £0  =  0  should  produce  the  same  interface  pattern  as  was  seen  ear¬ 
lier  since  no  displacement  initially  wilL  leave  the  vortices  at  their  equilibrium 
positions  for  all  time.  Figure  B.4  shows  the  developing  temporal  interface  for 
such  a  case.  The  cusp  and  spiral  arm  behavior  discussed  earlier  is  apparent  and 
the  result  checks  with  the  £0  =  0  limit. 

Jimenez  (1980)  describes  the  interface  we  speak  of  as  being  like  a  rubber 
band  that  is  attached  to  fluid  particles  at  time  equals  zero  and  forced  to  stretch 
or  compress  as  dictated  by  the  prescribed  fluid  motion.  Even  for  steady  flow  the 
rubber  band  shape  will  continuously  change  with  time.  It  is  interesting  to  see 
how  much  the  rubber  band  gets  stretched  and  where  it  gets  stretched  the  most. 
The  impression  one  gets  from  Figure  B.4  is  that  a  great  deal  of  stretching  takes 
place.  Strain  is  an  important  consideration  if  this  model  is  taken  a  step  further 
by  treating  the  interface  as  a  strained  diffusion  flame.  Such  an  extension  would 
allow  computation  of  the  product  of  a  chemical  reaction  between  two  species 
separated  by  the  interface.  This  problem  is  not  solved  in  this  work.  However, 
the  strain  on  the  interface  is  computationally  simple  to  find  Consider  two 
points  on  the  interface;  points  (a)  and  (b)  in  Figure  B.2  for  example.  Let  the  dis¬ 
tance  between  (a)  and  (b)  be  L  at  t  =  0 .  At  some  later  time,  (a)  and  (b)  will  be 
some  other  distance  apart,  say  Lx .  The  strain  is  the  new  length  divided  by  the 
original  length.  Now  consider  many  points  all  with  the  same  L.  At  a  later  time 
the  spacing  between  the  initially  equidistant  points  is  proportional  to  the  local 
strain.  This  strain  is  displayed  in  Figure  B.5  for  the  interface  roll-up  in  Figure 
B.4. 
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The  data  In  Figure  B.5  are  computed  starting  with  300  equally  spaced  points 
extending  over  four  wavelengths.  This  means  the  points  are  0.04189  apart  at 
t  =  0  .  Each  plot  in  Figure  B.5  shows  the  spacing  between  points  as. a  function  of 
point  number  for  a  certain  elapse  tune.  Figure  B.5a,  at  the  upper  left,  is  the 
first  time  step  shown  and  time  increases  to  the  right  and  downward.  Note  the 
frequent  rescaling  of  the  plots.  Figure  B.5a  shows  that  at  first  there  are  almost 
equal  amounts  of  stretching  and  compressing,  The  minimum  in  the  strain 
corresponds  to  the  cusp  point  inside  the  cat's  eye.  As  time  transpires  (i.e.,  mov¬ 
ing  from  plot  to  plot  in  Figure  B.5)  you  see  that  a  net  stretching  begins  to  occur. 
Two  local  maxima  occur  corresponding  to  locations  on  the  two  spiral  arms  on 
opposite  sides  of  the  inside  cusp  point.  The  apparent  “turbulence"  in  these 
plots  is  most  likely  due  to  the  numerics  since  nothing  in  the  flow  could  cause 
this  kind  of  randomness.  The  true  strain  would  probably  look  like  a  smoothed 
out  version  of  the  curves  in  Figure  B.5. 

At  some  definite  time,  a  change  in  the  strain  takes  place.  In  Figure  B.5g 
(third  row,  first  column)  a  peak  is  observed  at  the  extreme  left.  Subsequent 
frames  show  this  peak  increases  rapidly  until,  for  the  last  plot,  it  dominates  all 
others.  This  strain  maximum  occurs  at  the  point  where  the  interface  intersects 
the  cat’s  eye  boundary;  point  (d)  in  Figure  B.2.  Consider  two  points  very  near 
point  (d);  one  inside  the  cat’s  eye  and  one  outside.  Point  (d)  moves  toward  the 
stagnation  point  with  time  and  eventually  ends  up  at  the  stagnation  point  in 
infinite  time.  The  point  inside  the  boundary  also  moves  toward  the  stagnation 
point  at  first  but  will  make  a  sharp  turn  to  the  right  at  some  time.  Similarly, 
the  point  outside  the  boundary  will  move  toward  the  stagnation  point  at  first 
turning  sharply  to  the  left  at  some  later  time.  Thus,  the  two  points  which  were 
initially  very  close  together  remain  close  to  each  other  for  a  while  and  then 
separate  by  large  distances  at  a  time  associated  with  point  (d)  approaching  the 
stagnation  point.  The  outside  point  will  leap-frog  from  vortex-to-vortex  to  che 


left  while  the  inside  point  will  orbit  around  the  vortex  above  point  (a)  in  Figure 
Q.Z.  For  infinite  time,  the  points  will  become  infinitely  far  apart.  Jimenez  (1900) 
refers  to  this  process  as  a  production  of  interface  at  the  stagnation  points.  In 
his  solutions,  yQ  =  0,  so  the  stretching  occurs  instantly.  In  our  model,  y0  <  0 
and  it  takes  a  finite  amount  of  time  for  the  strong  stretching  to  occur. 

B.7.4.  Results  for  the  Case  of  Pairing  Rows.  Consider  the  system  (B.00) 
through  (B  04)  when  the  vortices  are  displaced  initially  from  their  equilibrium 
position;  i.e.,  £ 0  *  0  •  The  developing  temporal  interface  is  shown  in  Figure  B.6. 
Once  again,  the  interface  pattern  looks  very  much  like  the  pictures  of  the  actual 
flow;  c  f.,  pairing  vortices  at  the  downstream  end  of  Figure  3.3.  The  cusp 
behavior,  outlined  earlier  for  £o  =  0,  is  similar  for  the  pairing  vortices.  A  cusp 
forms  in  the  inside  flow  for  both  vortices  involved  in  the  pairing  as  do  the  spiral 
arm  pairs.  The  outside  flow,  below  the  vortices,  develops  cusps  also  with  the 
right  vortex  possessing  a  more  pronounced  one  than  the  left  vortex  Something 
new  shows  up  in  the  last  two  frames  of  Figure  B.6.  A  cusp  forms  on  the  other 
side  of  the  vortex  which  started  on  the  left.  For  the  steady  row,  £0  =  0,  this  can 
never  occur.  In  that  problem,  the  portion  of  the  initially  straight  interface  that 
is  outside  the  cat's  eye  is  either  entirely  above  or  entirely  below.  In  other  words, 
the  outside  problem  was  separable  into  two  problems;  an  above  problem  and  a 
below  problem.  Once  a  particle  started  above,  it  remained  above  and  likewise 
for  the  below  problem.  When  the  flow  is  unsteady,  such  as  for  £o  *  0 ,  this  is  not 
the  case.  Figure  B.6  is  an  example  of  the  effect  of  unsteady  flow.  The  interface, 
outside  the  instantaneous  cat's  eye  at  t  =  0,  is  entirely  below.  The  new  cusp, 
seen  in  the  last  two  frames,  suggests  that  part  of  the  interface  has  crossed  the 
boundary  to  become  part  of  the  above  flow  since  that  is  the  only  place  such  a 
cusp  could  form. 


B.S.  Summary 

The  roll-up  of  a  passive  interface  by  rows  of  point  vortices  has  been  studied 
for  the  case  when  the  interface  initially  starts  below  the  vortex  centers.  Many 
features  of  the  real,  forced,  pre-mixing  transition  shear  layer,  shown  in  Figure 
3.3,  appear  in  the  solutions.  Moreover,  these  features  exist  only  for  the  off-axis 
interface  and  not  for  the  symmetric  case  of  Jimenez  (i960). 

Cores  of  pure  high  speed  fluid  were  found  as  wens  the  cusp  and  spired  arm 
features  evident  in  Figure  3.3.  In  addition,  cusp  development  in  the  braids  also 
occured  which  matched  well  with  the  picture.  Study  of  the  strain  on  the  inter¬ 
face  showed  that  it  took  a  finite  length  of  time  for  the  large  strain  associated 
with  the  stagnation  point  flow  to  exert  am  influence  on  the  interface.  In 
Jimenez’s  solution,  the  influence  was  instantaneous  because  the  interface  inter¬ 
sected  the  stagnation  point  at  t  =  0 .  Study  of  orbiting  rows  of  point  vortices 
showed  interface  patterns  similar  to  those  observed  for  pairing  vortices  in  the 
real  shear  layer.  Specifically,  a  cusp  forms  on  the  high-speed-side  like  those 
visible  in  Figure  3.3. 
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Figure  1.1  :  Breidenthal 's  Product  Measurements  (unforced  flow) 
P  -  Product  Thickness  .  8  -  Vorticity  Thickness 


Growth  Law  for  the  Forced  Shear  Layer  (from  Browand) 
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Figure  2.1  :  Forcing  Meachanism  Rotor  and  Stator 

(a)  Maximum  Open  Area 

(b)  Minimum  Open  Area 
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Figure  2.3:  Installed  Forcing  Mechanism 
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To  High  Speed 
Stream 


Figure  2.4  :  SohCMtlo  of  Forcing  MochsnlM 


Figure  2.5:  Photograph  of  the  Laboratory 


(a)  U4  =  2.0  cm/sec  .  U2  =  1.0  cm/sec  ,  Reej  -  31 


(b)  -  14.8  cm/sec  .  U2  -  8.6  cm/sec  ,  Re9i  «  86 


Figure  3.1  :  Three  Flow  Regimes  of  the  Unforced  Shear  Layer 
(a)  Laminar  ,  (b)  Transitional  .  Cc)  Turbulent 
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Figure  3.2  :  Mean  Velocity  Profile  in  tne  Laminar 
Shear  Layer  at  x  -  20  cm. 
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Figure  3.3:  Non-reacting  L.I.F.  Photograph  of  a  Forced  Transitional 
Shear  Layer.  U.  «B  cm/sec  ,  U,«5  cm/sec 


(b)  forced  pairing 


Figure  3.4  :  Controlled  Amalgamation  dy  SuDharmonic  Excitation 


Figure  3.S  :  Effects  of  F  -  0.5  Hz  forcing  on  a  Laminar  Shear  Layer 
11^2  cm/sec  ,  U8-l  om/sec  .  (a)  Simultaneous  Plan  and 

Side  Views  (b)  Oblique  View  looking  down  and  upstream 


Figure  3.7  :  Side  Wall  Interaction  for  a  Transitional  Shear  Layer 
Uj  ■=  18  cm/sec  .  U2  ■>  13  cm/sec  .  r*=0.72.  Refl  -  90 


Effects  of  Forcirtfl  on  Shear  Layer  Structure  for 
Ut  -  17.4  cm/sec  and  U8 -  8 . 3  cm/sec  (part  i) 


Figure  3.B:  Effects  of  Forcing  on  Shear  Layer  Structure  for 
Uj  -  17.4  cm/aec  and  UE- 8.3  cm/sec  (part  2) 


Figure  3.8:  Effects  of  Forcing  on  Shear  Layer  Structure  for 
U1  -  17.4  cm/sec  and  U8  ■  8. 3  cm/sec  (part  3) 


Figure  3.9:  Simultaneous  Plan  and  Side  Views  of  the  Unforced  Transitional 
Shear  Layer,  U, =  17  cm/sec  .  U- ■=  8  cm/sec 


gure  3.10:  Simultaneous  Plan  and  Side  Views  of  a  Transitional 
Shear  Layer  Forced  at  2.4  Hz,  Uj  «=  17  cm/sec  , 

Ua  =  0  cm/sec 


(a)  unforced  flow 


(t>)  F  -  12  Hz  ,  F/F0-0.19P  X  =  4  cm 


a 


0  (d)  F  =  4  Hz  .  F/F0-0.06.  X  =  12  cm 


L  Figure  3.12  :  Time  exposures  of  the  Phenolphthalein  Reaction 

I  •  Ji  “  70  cm/sec  ,  U3  =  27  cm/sec 
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20  25 


r  -  0 . 38 
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Figure  3. iB  :  Product  Thickness  and  Growth  Rate  from  Absorption 
Measurement  —  High  Reynolds  Number  ,  F  ■  3.9  Hz 
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Figure  3.18  :  Product  ThicKness  and  Growth  Rate  from  Absorption 
Measurement  —  Low  Reynolds  Number  .  Unforced 


Uj  -  17.4  cm/sec 


r  -  0 . 48 


=  0 .54  mm 


Figure  3.21  :  Product  Thickness  and  Growth  Hate  from  Absorption 
Measurement  —  Low  Reynolds  Number  .  F  ■  4.5  Hz 
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Figure  3.22  :  Product  Thickness  and  Growth  Rate  from  Absorption 
Measurement  —  Low  Reynolds  Number  .  F  ■  5.8  Hz 
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Figure  3.24  :  Product  Thickness  Fluctustion  from  Absorption  Product 
Messurement  —  High  Reynolds  Number  .  F  -  3.8  Hz 
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Figure  3.25  :  Product  Thickness  Fluctustion  from  Absorption  Product 
Measurement  —  High  Reynolds  Number  ,  F  ■  0.0  Hz 
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Figure  3.27  :  Product  Thickness  Fluctuation  from  Absorption  Product 
Measurement  —  Low  Reynolds  Number  ,  F  -  1.8  Hz 
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Fiflure  3.29  :  Product  Thickness  Fluctuation  from  Absorption  Product 
Measurement  —  Low  Reynolds  Number  ,  F  -  4.5  Hz 
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Figure  3.32  :  Mean  Product  Concentration  Profiles  for  a  Transitional  Layer 

U.  -  17 . 4  cm/sec  .  U2  -  8.3  cm/sec 


TIME  (seconds) 


Figure  3.33  :  Time  Series  of  Product  Concentration  in  the  Unforced  Flow 

U,  -  70  cm/sec  ,  Ua  -  27  cm/sec 
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Figure  3.34  :  Time  Series  of  Product  Concentration  in  a  Forced  Layer 
F  -  0  Hz  .  Uj  -  70  cm/sec  ,  U2  -  27  cm/sec 
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Figure  3.36  :  Time  Traces  of  Product  Concentration  in  a  Transitional  Layer 
Data  shown  are  for  the  y  location  of  maximum  mean  product. 

U1  -  17.4  cm/sec  ,  U2  -  8.3  cm/sec 


REPRODUCED  AT  GOVERNMENT  FVPFNSE 


(c)  F  =  12  Hz  Xwo  “  2.8 


300 


Figure  3.37  :  Intensity  Transformation  Table  used  for  Figure  3.36 


(a)  unforced  flow 


Figure  3.38  :  Flow  Images  from  L.I.F.  Measurements  in  a  Transitional  Layer 

Uj  ■»  17.4  cm/sec  ,  U2  «*  8.3  cm/sec 


Flour*  3.41:  Mean  Hlstogr 


Figur*  9.42:  Hsan  Histogram.  U.  -  17.4  os/sao,  Ua  -  B .  3  c«/aec, 


Figure  3.43a:  Ratio  of  Forced  Product  Thickness  to  Natural  Product  Thickness 
as  a  Function  of  Oimensionless  Downstream  Distance.  Xm. 

High  Reynolds  Number  Results. 
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High  Reynolds  Number  Results. 


Figure  3.45  :  The  'Disappearing  Vortex'  Experiment.  Uj - 70  cm/sec . 

Ua  -  27  cm/sec .  Reacted  fluid  is  black  only  when  the 
ratio  of  high-speed  to  low-speed  fluid  in  the  mixed 
region  is  greater  than  five. 


(a)  Side  and  Plan  Views  of  Phenolphthalein  Reaction 


(b)  Reacting  L.I.F.  Cross-section.  Scale  different  than  (a)  . 


Figure  4.1  :  Unforced  Wake  at  a  Reynolds  Number  of  160 


Figure  4.2:  Simultaneous  Plan  and  Side  Views  of  a  Turbulent  Wake 

Ub"50  cm/sec,  Unforced 


(a)  Plan  and  Side  Views 


Effects  of  Forcirg  on  a  Laminar  Wake  (F  ■=  0 
(a)  Simultaneous  Plan  and  Side  Views  (t>)  00 
—  looking  down  and  slightly  upstream 
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v  v  - 
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(f)  F  -  10.8  Hz 


Figure  4.4  :  Effects  of  Forcing  on  a  Transitional  Wake 
U„«  12  cm/sec  ,  Re#  -  160 


Figure  4.5  :  Effect  of  Forcing  on  Spacing  Ratio  in  the  Initial  Karmen 
Vortex  Street  of  a  Transitional  Wake  (  U_  -  12  cm/sec  ) 
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Figure  4.B:  Reacting  Phenolphthalein  Pictures  of  the  Re#  -  1B0  Wake. 

U„-  12  cm/sec  .  1.3  mm 


unforced  f 
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Figure  4.9  :  Photographs  of  a  Forced  "Turbulent"  Wake  (Re, 


320)  . 
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Figure  4.10  :  Simultaneous  Plan  and  Side  Views  of  the  Highly  Mixed  Wake 


Figure  4.11  :  Cross-section  View  of  the  Highly  Mixed  Wake  using 
Non-Reacting  L.I.F.  Visualization 
x  -  12  cm 


Figure  4.12  :  'Tipping'  of  Spenttlse  Vorticity  6y  the  No-Slip  Condition 

on  the  Test  Section  Side-Malls. 


Figure  4.13  :  Model  for  Motion  of  the  Streemwlse  Vortices  In  the  Highly  Mixed 
Make  .  The  view  point  Is  from  dOMnstreem  looking  upstream. 

(a)  near  splitter  plate  (b)  Intermediate  stage 
(c)  steady  equilibrium  conf lguretion 


Figure  4.14  :  Effect  of  a  ‘False  Wall”  on  the  Highly  Mixed  Wake 
Simultaneous  Plan  and  Side  Views 


ing  L.I.F.  Measureme 
cm/sec,  x  ■=  25  cm.  (a 
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Figure  4.16  :  Mean  PDF  from  Non-Raactlng  L.I.F.  Maaaurementa 
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Figure  A. 2  :  Effect  of  Base  Concentration  on  the 
Absorption  Coefficient.  « 
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St re an line  Pattern  and  Initial  Interface  for  the  Temporal  Caae 


Figure  B.4:  Development  of  the  Temporal  Interface  for  a  Rom  of  Point  Vortlcea 
Dimensionless  elapse  time  la  Indicated  for  each  plot. 


Flgurs  B.S  :  Development  of  Strsin  in  the  Temporal  Interface  Solution 

shown  In  Figure  B.4  .  Abscissa  for  each  plot  Is  Point  Label 
Number  (Sosle  Is  0  to  300  for  each  plot) .  Ordinate  for  each 
plot  la  Dimensionless  Spacing  between  Initially  equidistant 
points.  Ordinate  scales  range  from  0  to  the  following  maxima 
for  each  plot  :  (a)  0.2  (b)  0.2  (o)  0.9  (d)  0.9  (e)  0.9 

(f)  0.9  (g)  1.0  (h)  1.0  (1)  1.0  ())  2.0  00  2.0  (1)  3.0 

Dimensionless  elapse  times  are  as  indicated. 


Figure  B.B:  Development  of  the  Temporal  Interface  for  a  Pair  of  Orbiting 
Rows  of  Point  Vortices.  Dimensionless  elapse 
time  is  indicated  for  each  plot. 


